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The  formulation  of  a  rigorous  statistical  mechanical  solution 
theory  relates  thermodynamic  properties  in  a  fluid  mixture  to 
volume  integrals  of  molecular  distribution  functions.   These 
distribution  functions  are  complex  functions  of  molecular  interactions 
and  the  macroscopic  state  of  the  system  including  the  composition 
although  for  many  purposes  the  composition  dependence  of  thermodynamic 
properties  can  be  represented  by  functional  expansions  about 
limiting  compositions  in  the  mixture. 

Molecular  distribution  functions  of  different  orders  are  related 
through  density  derivatives  of  the  distribution  functions.   A  first- 
order  perturbation  theory  has  been  developed  to  relate  the  distribution 
functions  of  a  real  fluid  to  the  similar  functions  in  a  hard  sphere 
fluid  through  differences  in  the  intermolecular  potentials.   The 
calculated  radial  distribution  function  for  a  Lennard-Jones  6-12  model 
agrees  well  with  molecular  dynamics  results.  A  one-parameter  correspond- 
ing states  theory  for  the  compressibility  integral  gives  qualitative 
predictions  for  thermodynamic  properties  in  subcritical  systems. 
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Corresponding  states  correlations  involving  only  a  single 
characteristic  parameter  for  each  substance  express  the  isothermal 
compressibility  of  pure  liquids  and  the  partial  molar  volumes  of  gases 
at  infinite  dilution  as  universal  functions  of  the  reduced  solvent 
density.   A  generalized  isothermal  equation  of  state  reproduces 
pressures  and  volumes  in  compression  experiments  of  pure  liquids  up 
to  high  pressures. 

The  isothermal  compressibility  of  liquid  mixtures  is  described 
by  a  one-fluid  theory  with  composition  dependent  parameters.   A 
postulated  linear  variation  of  direct  correlation  function  integrals 
between  composition  extremes  gives  accurate  predictions  of  partial 
molar  volumes  and  activity  coefficients  at  intermediate  compositions. 
The  experimental  triplet  correlation  functions  in  a  representative 
subcritical  system  did  not  agree  with  the  Percus-Yevick  or  HNC 
approximations . 

An  experimental  determination  of  the  volumetric  and  vapor- 
liquid  equilibrium  properties  of  gas-solvent  mixtures  is  reported. 
Volumetric  properties  in  the  liquid  phase  are  expressed  by  functional 
expansions  in  the  solute  mole  fraction,  with  a  one-fluid  model  for  the 
solvent  mixture.   Generalized  correlations  are  developed  from  experimental 
data  to  relate  the  activity  coefficient  parameters  of  H„  and  CH,  in 
liquid  solvents  to  the  reduced  solvent  density.   Several  expressions 
relating  the  Henry's  constant  of  a  gas  in  a  mixed  solvent  to  those  in 
the  pure  solvents  are  analyzed  with  experimental  data.  Vapor  liquid 
equilibrium  calculations  show  that  some  of  these  expressions  with  the 
generalized  activity  coefficients  yield  K-factors  and  bubble  pressures 
in  good  agreement  with  experimental  data. 
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CHAPTER  1 
THERMODYNAI-IICS  OF  GAS-LIQUID  SYSTEMS 

Our  everyday  experience  presents  a  panorama  of  processes  in 
gas-liquid  systems.   They  vary  in  nature  from  the  complex  transfer 
of  oxygen  into  the  blood  stream  of  animals,  to  the  simple  dissolution 
of  carbon  dioxide  in  the  ubiquitous  coke  bottle.   A  major  effort  of 
chemical  and  petroleum  industries  is  contacting,  treating  and  separating 
gas  and  liquid  streams. 

In  all  of  these  instances,  the  equilibrium  distributions  of 
material  in  the  two  phases  follow  simple  thermodynamic  relations.   This 
is  not  to  imply  that  equilibrium  conditions  obtain  in  these  processes, 
for  then  there  would  be  no  transfer  of  materials.   The  equilibrium 
distributions  reflect  the  thermodynamic  constraints  on  the  system. 
In  the  engineering  domain,  these  equilibrium  distributions  are  a  key 
factor  in  the  selection  of  processing  conditions  and  the  design  of 
process  equipment. 

Thermodynamic  properties  of  matter  in  bulk  are  the  manifestation 
of  the  molecular  interactions  between  the  different  species  present 
in  the  system.   Statistical  mechanics  provides  the  formal  link  between 
the  molecular  characteristics  of  a  physical  system  and  its  bulk 
thermodynamic  properties.   This  work  describes  the  use  of  a  statistical 
mechanical  solution  theory  based  on  molecular  distribution  functions 
for  correlating  the  thermodynamic  properties  of  gas  solvent  systems. 

1.1  Phase  Equilibrium  in  Gas-Liquid  Systems 

Generally,  for  thermodynamic  equilibrium  between  the  phases  of  a 


multicomponent  system,  the  chemical  potentials  of  each  species  are 
equal  in  all  phases,   i.e., 

yj"'   is  the  chemical  potential,  or  partial  molar  free  energy  of 

species  i  in  phase  j,  and  depends  on  temperature,  pressure  and  composition 

of  the  phase. 

For  any  phase,  solid,  liquid  or  gas,  the  chemical  potential  of 

the  individual  species  can  be  expressed  in  terms  of  the  fugacities,  f., 

2 

as 

y^  -  uj  -  RTin  f^/f°  (1.2) 

y.  and  f ,  are  the  chemical  potential  and  fugacity  in  the  standard 
state.   Of  these  two  quantities,  one  can  take  on  an  arbitrary  value, 
but  they  are  not  both  independent.   The  fugacity,  f .,  also  depends  on 
temperature,  pressure,  and  composition.   The  fundamental  equilibrium 
relation  (1.1)  reduces  to 

,(1)  .  ,(2)  _  ,(3)  (,.3) 

The  central  problem  of  phase  equilibrium  thermodynamics  is  the 
description  of  the  component  fugacities  in  terms  of  temperatures, 
pressures,  and  compositions.   For  vapor-liquid  equilibrium, where  the 

species  are  liquids  at  the  temperature  of  the  system,  a  suitable 

2 
description  of  equation  (1.3)  is 

P   -L 
.v  «       ,„refv    ^o^     I     i 


^^y^P  =  x^Y.  (P"^)  .  f°^  exp  J    ^  dP  (1.4) 

,ref 


P 


which  relates  the  vapor  and  liquid  phase  fugacities.   $  is  the  vapor 

phase  fugacity  coefficient  and  represents  the  deviation  from  ideality 

of  the  individual  components  in  the  mixture,  and  the  nonideality  of 

mixing.   For  an  ideal  mixture  of  ideal  gases,  i.    is  unity.   f°  is  the 

1  1 

fugacity  of  the  liquid  in  the  standard  state;  this  is  the  pure 
component  at  the  temperature  of  the  system,  and  under  its  own  saturation 
pressure  (p'^^  =  P^^^) .   y^(P^^^)  is  the  activity  coefficient  in  the 
liquid  phase;  it  represents  the  composition  dependence  of  the  chemical 
potential  at  the  reference  pressure.   (The  importance  of  this  quantity 
in  the  thermodynamics  of  liquicfe  cannot  be  over  emphasized.)   For  an 
ideal  solution,  y^  =  1  at  all  x^.   In  a  nonideal  mixture,  y.(?^^^)    is 
normalized  such  that 

Y^  -►  1   as   x^  -^  1  (1.5) 

The  activity  coefficient  is  unity  when  the  component  is  in  its  standard 
state.  The  exponential  term  corrects  for  the  change  in  fugacity  of  the 
liquid  from  the  reference  pressure  to  the  pressure  of  the  system. 

The  choice  of  the  pure  component  as  the  standard  state  for 
fugacities  serves  very  well  in  the  description  of  liquid  mixtures 
wherein  each  species  can  exist  in  the  pure  state  at  the  system 
temperature.   However,  this  definition  is  unsuitable  in  considering 
gas-liquid  systems,  for  which  the  lighter  components  are  super- 
critical at  the  mixture  temperature.   Here,  the  solute  at  infinite 
dilution  in  the  solvent  can  be  chosen  as  the  standard  state;  the 
solvent  may  be  a  pure,  or  mixed,  solvent.   The  liquid  phase  fugacity 
for  the  noncondensable  components  becomes 


r.L  *,  ref.    /r,ref. 

fi  =  x.Yj^(P   )H^(P   )  exp 


P   -L 

RfdP  (1.6) 


pref 

where  the  reference  fugacity  H.   ,  commonly  called  the  Henry's 
constant,  is 

x,-»-0  1    ^^ 

*  ref 

The  activity  coefficient  y.(?        )  is  now  normalized  such  that 

Y*(P"^)  -^1   as   x^  ->  0  (1.8) 

ref 
where  P    is  usually  the  bubble  pressure  of  the  solvent  at  the 

temperature  of  the  system.   The  product  x.y.CP'^^  )  H.  (P^^  )  describes 

i  i       i^ 

the  variation  of  the  liquid  phase  fugacity  with  composition. 

An  important  relation  exists  between  the  Henry's  constant 
and  activity  coefficient  of  a  gas  dissolved  in  a  mixed  solvent  and  the 

corresponding  quantities  for  the  gas  dissolved  individually  in  each  pure 

3 
solvent.   It  is 

Jin  H    -I     x.£nH,.   =«  T  x ,  (lim  2.nY..)  -  lim  Iny.  (1.9) 

j=l  "^        j=l    ^i"*"^        ^i 

where  the  multicomponent  solvent  mixture  consists  of  N  solvents  with 
mole  fraction  x . ,  j=l,...,N.   H,   and  H.  .  are  the  Henry's  constants 

of  the  gas  in  the  multicomponent  solvent  and  the  jth  solvent, 

*       * 

respectively;  y       and  y .      are  the  respective  activity  coefficients. 

Then  if  the  solute  interaction  with  the  individual  solvents  is  known, 
equation  (1.9)  can  be  used  to  calculate  the  solute  nonideality  in  a 


multicomponent  solvent  mixture  from  the  Henry's  constant,  which  is 
an  ideal  solution  property.   Conversely,  if  the  solute  nonideality  in 
the  mulricompcnent  solvent  can  be  approximated  by  some  average  of  the 
nonideality  in  the  individual  solvents,  the  solubility  in  the  multi- 
component  solvent  can  be  determined  from  those  in  the  individual 
solvents.   However,  one  must  recognize,  the  important  and  unavoidable 
assumption,  made  in  equation  (1.9),  that  the  composition  dependence  of 
the  activity  coefficient  as  determined  from  dilute  solution  (or  low 
concentration)  behavior  can  be  extrapolated  to  the  physically 
unrealizable  limit  of  x  =  1. 

1.2   Isothermal  Pressure  Dependence  of  Thermodynamic  Properties 

The  fugacity  of  a  component  in  a  mixture  depends  on  the 
temperature,  pressure  and  composition.   Its  pressure  dependence  is 
given  by  the  partial  molar  volume 


SZnf . 

1 

9P 


T,x 


1 
RT 


(1.10) 


The  derivative  is  taken  at  constant  temperature  and  composition;  the 
volume  changes  to  satisfy  the  general  equation  of  state  of  the  system. 
(The  variance  of  the  intensive  quantities  is,  of  course,  determined 
by  the  phase  rule.)   For  a  finite  pressure  change. 


£n  f, 


-  «.n  f, 


2  - 


(1.11) 


1  P, 


In  most  liquid  systems,  the  partial  molar  volume  does  not  vary 
appreciably  with  pressure  and  (1.11)  is  simplified  to 


£n  f. 


-  Jin  f. 
2         ^1 


RT 


(1.12) 


When  the  vapor  and  liquid  phases  are  assumed  ideal,  and  the  partial 
molar  volume  is  that  at  infinite  dilution,  we  have 


^         H^(P"^  =  0)     v°P 


y,P 


exp 


RT 


(1.13) 


which  gives  the  pressure  dependence  of  the  solubility  in  an  ideal 
liquid  solution. 

The  effect  of  pressure  changes  on  a  fluid  mixture  is  given  by 
the  isothermal  compressibility,  <,  which  is  defined  as 


8P 


T,n 


(1.14) 


where  n  denotes  the  constant  composition  of  the  system.   In  integrated 
form,  equation  (1.14)  is 


AP  -  P2  -  P^ 


P< 


(1.15) 


Given  k  as  a  function  of  density,  this  equation  can  be  used  to  relate 
corresponding  changes  in  pressure  and  density.  Equation  (1.15)  is  an 
isothermal  equation  of  state. 
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CHAPTER  2 
DISTRIBUTION  FUNCTION  SOLUTION  THEORY 

A  distribution  function  solution  theory  for  equilibrium 
properties  of  multicomponent  fluid  systems  was  first  proposed  by 
Kirkwood  and  Buff.   By  rigorous  statistical  mechanical  methods, 
they  related  thermodynamic  properties  of  a  mixture  to  volume 
integrals  of  the  radial  distribution  functions.   The  shape  of  the 
radial  distribution  functions  in  a  fluid  mixture  result  from 
intermolecular  attractions  and  depend  in  a  complex,  and,  as  yet, 
incompletely  determined,  way  on  the  macroscopic  variables  of  the 
system  and  on  the  intermolecular  potentials  of  the  different  species 

present.   Currently  available  theories  for  predicting  the  distribution 

2 

functions,   albeit  approximate,  can  provide  an  adequate  basis  for 

determining  thermodynamic  properties  through  the  use  of  statistical 
mechanical  relations  to  be  presented  here.   The  properties 
considered  are  the  isothermal  compressibility,  the  partial  molar 
volume  of  each  component,  and  the  composition  derivative  of  the 
chemical  potentials. 

2.1  Radial  Distribution  Function  Solution  Theory 

The  framework  of  the  grand  canonical  ensemble  is  used  in 


the  development  of  the  solution  theory.   Relations  are  obtained  between 
density  fluctuations  and  integrals  of  the  radial  distribution  functions, 
and  again  between  density  fluctuations  and  the  composition  dependence 
of  the  chemical  potentials.   Elimination  of  the  density  fluctuations 
leaves  the  desired  relationships. 

The  system  of  interest,  with  volume  V,  is  considered  to  be  in 
a  large  heat  bath  and  open  with  respect  to  molecules  in  the  system. 
The  bath  provides  a  reservoir  of  heat  at  temperature  T  and  of  molecules 

at  chemical  potentials  u^,   u^ u„.   The  system  is  then  characterized 

by  the  thermodynamic  variables  T.  V.  u,.  i  =  1.  M.   (M  is  the  nu.ber 
of  chemical  species  present.)  There  are  fluctuations  in  the  number  of 
molecules.  N^.  of  each  species  within  the  system.   The  average  singlet 
and  pair  densities  and  the  fluctuations  in  the  density  of  each  species 
are  related  through  equations  of  the  form 


! 


P^   (rJdr,  =•  <N  > 
a   —1  —1     a  av 


(2.1) 


«o   J.  z   I  —z  a  6  av  \^   a  av  \ 


av 


^^1% 


<N  N  >   -  <N  >   <N  > 
a  3  av     Qc  av  S  av 


(2.2) 


-  6   <N  > 
0(6  a   av 


(2.3) 


Here 


(2) 


as  -l'-2-^  ^^  ^^^   average  pair  density  for  mole 


(1) 


cules  of  species 


type  a, 


a  and  e  and  p^   (r^)  u   the  singlet  density  for  molecules  of 

<N^>  is  the  average  number  of  molecules  of  species  a  within  the  volume  V. 


10 


(2) 
related  to  the  average  pair  and  singlet  densities: 


The  radial  distribution  function   g''^'' (r^  ,r^)  or  g  „  (r) ,  is  simply 

dp     ^1   i.  CIO 


P^f  (r^.r^)  =  P^'^r^)P3^'^r2)g^f  (r^.r^)  (2.4) 


c  =  p^"^\r,)  =  <N  >   /V  (2.5) 

a    a   —1      a  av 


Ir^-r^I  (2.6) 


c  is  the  bulk  concentration  of  species  a;  r  is  the  scalar  difference 
a 

(2) 
of  r,  and  r„.   Equations  (2.4)  and  (2.6)  assume  that  g  „  (r)  depends 
—1     — Z  ap 

only  on  the  relative  distance  between  r^^  and  r^^ ,  and  not  on  their 

(2) 
exact  locations  within  the  system,   g   (r)  is  generally  a  function  of 

ap 

density,  temperature,  composition  and  the  intermolecular  forces  between 
molecules  of  species  a  and  6. 

Equations  (2.4)  and  (2.3)  together  yield  the  desired  relation 
between  density  fluctuations  in  the  grand  ensemble  and  the  integrals 
of  the  radial  distribution  functions 


a  av  3  av  a 


where  the  integral  extends  over  all  sets  of  relative  coordinates  of  the 

(2) 

pair  a  and  B.   It  is  instructive  to  note  that  g   (r)  as  written 

represents  systems  of  monatomic  fluids  only,  where  the  molecules  have 
only  translational  degrees  of  freedom.   If  the  system  under  study 

consists  of  molecules  with  structure,  these  must  be  accounted  for  in  the 

(2) 

expression  for  g        (r)    and   the  integration  extended  over  all  possible 

ap 
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(2) 

configurations.   In  the  present  discussion,  g   (r)  is  considered  to 

ap 

include  the  averaged  effect  of  all  degrees  of  freedom  of  the  molecule. 

In  the  grand  canonical  ensemble  the  relation  between  density 
fluctuations  and  the  composition  dependence  of  the  chemical  potentials 


<N  N„>   -  <N  >   <N  > 

a  B  av     a  av  B  av 


\\,\/W 


(2.8) 


where 


ae 


kT 


3N, 


(2.9) 


T.V.N 


and  |A  I  is  the  cof actor  of  A  ^  in  the  determinant  det/A/;  \i      is  the 
chemical  potential  per  molecule  of  species  a.   Then,  eliminating  the 
density  fluctuation  between  (2.8)  and  (2.7),  we  get 


where 


1^ 
kT 


3^ 

3N, 


T.V.N 


\al/"l'l 


Bo-c6„  +  c  c„G   „ 
a6    a  ag    a  S  a6 


=a6"l^«a 


(2) 

6 


(r)  -  l}dr 


(2.10) 


(2.11) 
(2.12) 


Equation  (2.10)  deals  with  fluid  mixtures  at  constant  volume  and 
temperature.   To  obtain  relations  at  constant  pressure,  we  use  the 
Gibbs  Duhem  equation 

M    c  3y 


K  "4  ^ 

a-1   ^   g  J 


(2.13) 


T.P.N 


and  a  mathematical  relation 
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a 


T,V,N 


c 

9N 


6 


T,P,N 


g  B 


(2.14) 


Y  '  '  Y 

where  v  is  the  partial  molar  volume,  per  molecule,  of  species  a,  and 
K   is  the  isothermal  compressibility  of  the  mixture.   Use  of  the  additive 
property  of  the  partial  molar  volumes 


M 

I 
a=l 


N  V     =  V 

a  a 


(2.15) 


results   in 


KkT  =  |B|/y  y  c  cJb   I 


(2.16) 


v^  =  {  V  c  |B  I  }/y  y  c  cJb  I 

a  a  8 


3u 

c 

3N. 


T,P,N 


kT 

vTb] 


3  D 


a6' 


I  I   <^  S|S 

^    i.        n    ft  ' 


6  a 


a  S'    aB 


(2.17) 


(2.18) 


Equations  (2.16),  (2.17)  and  (2.18)  are  the  key  equations  of  the  solution 

theory.   For  ease  of  manipulation,  they  are  conveniently  expressed  in 

.  .  3 

matrix  form.   Redefining 


N 


ctB 


(g^f(r)  -Ddr 


and  further 


aB   kT 


r  9y 

c 

9N, 


T.P.N 


=  IL 

^aB   kT 


3p 

c 

3N, 


T.V.N 


equations  (2.16),  (2.17)  and  (2.18)  are  written  as 
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v/<RT  =  i'^B"-'-  X  i  (2.19) 


"    1   t  -1 
pv.  =   y  B.,/±     B  X  i 


(2.20) 


where 


X  D  =  I  -  ixi^'h^   i  i^)/±^   b'^x  i  (2.21) 


B  =  I  +  X  G  (2.22) 


A,  D  and  G  are  matrices  with  elements  a   ,  d   ,  and  G      ,   respectively, 
X  is  a  diagonal  matrix  of  mole  fractions,  x   -  N  /N,  i_  is  a  unit  column 
vector.   The  basic  equation  (2.10)  relating  chemical  potentials  to  the 
radial  distribution  function  Integrals  is 

A  -  (B  x)"-"-  (2.23) 

2.2  Direct  Correlation  Functions 

O'Connell  has  shown  that  integrals  of  the  direct  correlation 
function  may  be  used  in  the  distribution  function  solution  theory.   The 

direct  correlation  function  is  related  to  the  radial  distribution 

2 

function  by  the  Ornstein-Zernike  integral  equation 


g(r)  -  1  =  c(r)  +  p 


c(r,r'){g(0,r')  -  l}dr'         (2.24) 


Like  g(r),  c(r)  is  also  a  function  of  the  intermolecular  potential  and 
the  macroscopic  variables  of  the  system. 

There  is  no  conceptual  advantage  in  relating  thermodynamic 
properties  to  integrals  of  the  direct  correlation  functions  rather  than 
to  those  of  the  radial  distribution  functions.   However,  on  the 
molecular  scale,  as  r,  the  distance  between  molecular  centers  increases, 


lA 


the  function  c(r)  goes  to  zero  much  more  rapidly  than  the  function 


g(r) 


2,6 


Hence,  it  may  be  possible  that  the  Integral 


c(r)dr 


can  be  calculated  more  accurately  from  a  molecular  theory  than  the 
corresponding  integral  {g(r)  -  l}dr_.   Anticipating  these  developments, 
the  further  results  of  the  distribution  function  solution  theory  are 
presented  in  terms  of  integrals  of  direct  correlation  functions. 

The  spatially  integrated  form  of  the  multicomponent  Ornstein- 
Zernike  integral  equation  is  (Appendix  A) 


G  -  C  +  G  X  C 


where 


"^ij    V 


clf(r,.r2)dr2 


(2.25) 
(2.26) 


Then  from  (2.23)  and  (2.25)  we  have 


xA  =  B~  =I-xC 


(2.27) 


and  the  equations  for  thermodynamic  properties  become 


M 


pv  =  {1  -  I      X  C   }/{l  -I        I      C       X  x^} 
j=l  ^   ^       j=l  k=l  J^  J  ^ 


(2.28) 


and 


l/p<kT  =  {1  -  y   Y   x.x,  C,  } 
j=l  k=l  J^  J^ 


X  D  = 


x  F  i  i^  ' 


I  - 


i  I  X  i 


(2.29) 
(2.30) 


where 


F  =  I  -  X  C 


(2.31) 


Equation  (2.30)  has  elements  of  the  form 
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kT 


3U, 


X  {l+C        -  I  X   (C     +C      )+  llx  X    (C.    C,    -C..C      } 
1  IJ        ^     n     In     jn       ^  ^  n  m'   m  jm     ij   nm 


r _n n  tn 


•J  V  r 

TPN  l~Zz.xxC 

^  n  m 

(2,32) 

In  principle,  one  can  use  currently  available  theories  for  the 

distribution  functions  in  the  above  expressions  to  calculate 

thermodynamic  properties.   However,  the  volume  integrals  of  the  direct 

correlation  functions  C    (and  the  radial  distribution  functions  G  ), 

•J  ij 

depend  on  the  macroscopic  variables  of  the  system.   Their  composition 
dependence  must  be  known  before  equations  (2.32)  or  (2.18)  can  be 
integrated  to  obtain  the  Gibbs  free  energy  of  the  mixture. 

However,  this  restriction  is  by  no  means  absolute.   As  in 
other  branches  of  mathematical  physics,  one  avenue  of  progress  is  via 
functional  expansions  about  those  limiting  compositions  in  the  mixture 
at  which  the  distribution  functions  can  be  determined  with  some  degree 
of  certainty,  e.g.,  the  pure  components.   Then,  within  reasonably  small 
ranges  about  these  compositions,  the  thermodynamic  properties  can  be 
evaluated  by  treating  the  distribution  functions  as  independent  of 
composition. 

2jJ — Taylor's  Series  Expansions  for  Thermodynamic  Properties 

The  integrals  C^^ ,  equations  (2.25)  -  (2.32),  are  generally 
dependent  on  the  system  temperature,  density,  and  the  composition.   The 
limiting  composition  dependence  is  determined  from  the  behavior  of  the 
pure  fluid  (solvent  component)  and  the  infinitely  dilute  solution 
(solute  component).   To  relate  the  discussions  to  the  thermodynamics 
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of  gases  in  mixed  solvents,  expansions  are  derived  for  the  following 
physical  systems: 

i)  Pure  fluid,  i.e,  pure  solvent, 
ii)   Solute  gas  in  single  solvent — binary  system, 
iii)   Solute  gas  in  binary  solvent — ternary  system. 
Component  1  is  the  solute,  2  and  3  are  solvents,  superscripts 
and   refer  to  solute  at  infinite  dilution  in  single  and  mixed 
solvent,  respectively.   The  general  expressions  will  be  derived  first, 
i.e.,  case  (iii),  and  the  results  for  (ii)  and  (i)  inferred  as  simplifica- 
tions of  the  general  case. 

2.3.1  Solute  Gas  in  Binary  Solvent 


The  general  form  of  functional  expansion  is 

2 


W). 


V^l 


W)     n  +  5^1 

X  =0    1 


3x, 


■) 


X,  r 


3  2w 


x,=0 


3x 


1  X  =0 


(2.33) 


where  W  is  the  thermodynamic  property  and  x^  is  the  solute  mole 

2,       ^ 


fraction.   The  derivatives 


9W 
3x 


1  i 


aw 

3x? 


are  obtained  from  relations 


(2.28),  (2.29)  and  (2.32),  respectively.   The  integrals  C   are  functions 
of  X- ,  and  the  derivatives  include  integrals  of  higher  order  direct 
correlation  functions.   The  relations  between  integrals  of  different 
orders  are  presented  in  Appendix  B.   The  relation  used  here  for  the 
composition  derivatives  is  equation  (B.9).   The  final  results  are: 


itXYj  = 


^+/o      2.    ^+,  3  2    3v  . 
-  f2^^^i  ~  ^1^  ~  3^  Y  ^1  ~  ''l^    ' 


(2.34) 
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where 


«t- 


-a-ct,,.{^ 


(2.35) 


34  = 


-    (1  -   C 


11 


+  -+ 

1  +  C     {1  +  u  -  pv^d 

<11>   ^  " 


-  w)} 


(1  -  D) 


_+ 


_+, 


-+. 


+  pv^[2(l   -   2t)(l  -   pv^)    -   3r  +  pv^{3u  +  D  -  pv^(D  +w)}] 

(2.36) 


pv^  =  pv^  +  x^ 


^              _+       {pv^[2u+t+D-pv^(EHw)]-C^^-r}  " 
2pv^(l  -  pv^)   +  


(1  -  D) 


+    ... 


(2.37) 


_+ 


PkRT 


i-  =    (1  -   D)   +  X    [2(D   -    t)    -  u  +  wpv^  -  D(l  -   pv^)]   +    ... 


where        D  ?   ^^022  +  ^3^33  ■*"  ^^2^3^23 


(2.38) 
(2.39) 


t  E  X2C^2  "^  ^"3^13 


u  =-   x2c^22  ^  2X2X30^23  "■  ^3^133 

w  E   x2c^22  +  3x^x3^23  "^  ^^2^3^33  "-  4^333 


r  E  x^c[^^  +  X3C^^3 


From  equation    (2.38)   with  x^  =  0 


PkRT 


(1  -  D) 


X  =0 


(2.40) 


re 


lating  (1  -  D)  to  the  compressibility  of  the  solvent  mixture.   Again, 


from  (2.37)  with  x  =■  0 
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^^^1^x^=0  ■  (P^l^""  = 


1  -  t 
1  -  D 


(2.A1) 


Thus  (1  -  t)  is  related  to  the  partial  molar  volume  of  the  solute  at 
infinite  dilution  in  the  solvent  mixture. 

2.3.2  Pure  Component 


Se 


tting  X  =  X  =  0  in  (2.38)  gives 


p<RT 


=  1  -  C 
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(2.42) 


relating  C^j   directly  to  the  isothermal  compressibility  of  the  pure 
solvent  2. 

2.3.3  Gas  in  Single  Solvent 

Setting  x^  =  0  in  (2.38)  and  (2.37)  yields 

^  =  (1  -  C°2)  +  x^[2(C°2  -  C°2)  -  C°22  ^  ^222  '\   "  ^^2^^  "  ^^1> 


(2.43) 


pv,  = 


pv°  +  x^ 


-C°^-C°,,.3pv°C°,.2C°,,PV°-pv°C°,-(pv°)2(C°2.C°,,) 


(1  -  C^2) 


where 


■  1  -  c 
-o       ^12 

p^i ': — :^ 


(2.44) 


(2.45) 


1  -  C 


22 


thus  C^2  is  related  to  the  partial  molar  volume  of  solute  1  at  infinite 
dilution  in  solvent  2.   Then,  with  x  =  0  in  (2.35)  and  (2.36) 


2f°  =  -  (1  -  C°  )  +  i2- 

^^    (1  -  C°  ) 


(2.46) 
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3f"  =  -  (1  -  C 
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^111^  -^ ^^ 


C°   )} 


1  -  C 


22 


+  pv°[2(l  -  2C°2)(1  -  pv°)  -  3C°^2  +  ^^l^^^m  "^  ^22  "  P^1^^22  "*"  ^222^^^ 


From  (2.34) 


3x, 


-2£) 


(2.47) 


(2.48) 


x^-0 


C   is  contained  in  f„,  and  is  thus  related  to  the  limiting  slope  of 
the  composition  dependence  of  the  solute  activity  coefficient. 

Equations  (2.42),  (2.45)  and  (2.48)  dictate  a  hierarchy  for  the 

o    o 
determination  of  the  direct  correlation  function  integrals  C  „.  C-2 

and  C^ ^  from  experimental  data.   It  is: 


pure  solvent  isothermal  compressibility  gives  C„„. 


C--  and  p.m.v.a.id.  give  C^.. 
,0   „o 


^22'  "^12  ^^^   limiting  slope  of  activity  coefficient  give  C^^. 


A  similar  hierarchy  exists  in  the  more  general  case  of  the  solute  in 

*       -    1 
binary  solvent.   With  experimental  values  of  (9£nYi/3x-),  pv.  ,  — =-;■ 

XX       X   p<Kl 

one  can  calculate  D,  t  and  C^ ^  from  equations  (2.40),  (2.41)  and 
(2.35),  respectively. 


However,  no  such  identification  can  be  made  between  the  thermo- 
dynamic properties  and  the  three  body  integrals  C  ,  and  C,  ,  .   This  is 
a  natural  consequence  of  the  Taylor's  series  expansion.   The  triplet 
functions  result  from  the  composition  and  density  derivatives  of  two 


body  integrals  in  the  composition  dependent  expressions  for  thermodynamic 
properties.  Therefore,  the  C.  must  be  estimated  from  experimental  data 
by  truncating  the  expansions  after  a  certain  number  of  terms. 

It  is  appropriate  to  :r.ake   scire  general  remarks  about  the 
expressions  derived  thus  far.   Being  expansions  about  zero  composition, 
they  should  be  valid  only  in  the  range  of  small  solute  concentrations. 
To  increase  the  composition  range  one  must  include  terms  from  second 
and  higher  derivatives  of  the  direct  correlation  function  integrals; 
these  will  include  triplet  and  higher  order  correlation  function 
integrals. 

The  expressions  for  isothermal  compressibility  and  partial  molar 
volume  indicate  an  important  difference  betveen  these  properties  and 
the  chemical  potential.   At  small  concentrations  where  the  linear 
expressions  in  x  are  valid,  the  chemical  potential  depends  on  two- 
body  correlations  only,  i.e.,  f„,  whereas  the  compressibility  and 
partial  molar  volume  expressions  include  two  and  three-body  terms. 

Equation  (2.42)  has  another  important  characteristic;  C-.  ^  , 

which  is  related  to  the  solute  nonideality,  appears  in  the  linear 

term  for  the  partial  molar  volume.   It  can  be  determined  from  the 

experimental  partial  molar  volume  in  its  linear  range;  i.e.,  from  the 

slope  at  zero  concentration.   This  presupposes  a  reliable  method  for 

determining  the  three-body  integrals  C^^_  and  C-^o*   ^  generalized 

method  for  doing  so  is  described  in  a  later  chapter.   (One  may  be 

tempted  to  make  the  simple  approximation  of  setting  the  C. .,  =  0, 

ijk 

but  this  is  not  always  physically  realistic.) 
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2.4  Linear  Composition  Dependence  of  Direct  Correlation  Function  Integrals 

It  is  possible  to  relate  the  thermodynamic  properties  of  the 
solute  in  the  binary  mixture  to  those  of  the  solute  in  the  pure  solvents 
by  postulating  simple  rules  for  the  composition  dependence  of  the 
quantities  C , . .   The  solution  theory  itself  does  not  provide  any  such 
prescriptions.   These  must  come  from  a  molecular  theory  for  mixtures. 
The  simplest  such  rule  is  a  linear  variation  of  the  form: 

* 
where  the  superscript   refers  to  the  solvent  composition  X_  =  1, 

and  the  superscript   refers  to  the  solvent  composition  X„  =  1.   A 

possible  variation  of  the  quantities  C  .  according  to  the  rule  (2.49) 

Is  shown  in  Fig.  2-1.   Quantities  that  can  be  determined  from 

experimental  binary  data  alone  are  shown  circled.   For  example, 

,0   „o 


Isothermal  compressibility  of  pure  2  gives  C-_;  C„„  and  partial  molar 
volume  of  3  at  Infinite  dilution  In  2  gives  C   ;  C_„,  C._  and  limiting 

slope  of  activity  coefficient  of  3  in  2  gives  C  _.   The  quantities 

o       * 

C^_  and  C^2  ^^^   different  from  the  two-body  terms  encountered  thus  far. 

o 
Cj^-  Is  the  Integral  of  the  direct  correlation  function  between  one 

molecule  each  of  1  and  3  when  both  are  entirely  surrounded  by  a  medium 

of  type  2.   It  is  not  directly  related  to  any  thermodynamic  property. 

A  naive  approximation  for  C^-  is  a  mean  of  C^^   ^^'^  ^to'  ^•®- »  simply 

the  arithmetic  mean. 

^13  ■  I  <^°12  -^  ^32^  <2.50) 

A  more  appealing  approximation  is 
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(C*3) 


(C*,) 


(C*3) 


12 


(0*3) 


32   "-23 


11 


'22 


'33 


(C°,) 


(C°,) 


«12> 


'13 


(C°3) 


Fig.  2-1.   Linear  composition  dependence  of  the 
integrals  c1".  in  system  of  solute  in 

binary  solirent.   Circled  quantities 

can  be  determined  from  binary  data  alone. 
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^13  "  ^3  •  ^2/^3 


(2.51) 


and  is  arrived  at  by  the  following  argument.   C^_  is  the  1-3  interaction 
in  an  environment  of  type  3;  C^ .  results  when  the  environment  is  changed 
to  type  2  and  this  is  achieved  macros copically  by  changing  the  volume 
in  the  ratio  of  the  pure  component  volumes.   This  assumes  that  the 
change  in  environment  is  accurately  described  by  size  effects.   When 
the  solvents  2  and  3  are  chemically  similar  and  are  composed  of  nearly 
equal  sized  molecules,  a  realistic  assumption  is 


^12  -  ^3 


(2.52) 


From  relationships  of  the  type  (2.49)  we  get 


PkRT 


o  o 


^2'^2^^   P3<3 


RT 


*      o 

+  2x^x^(1  -  ^2^23  ~   ^2^3^ 


-  X2X3(X2C22  -  X2C22  -  x^C^^  +  ^3^33^ 


(pv^) 


-  (pkRT) 


1  -  X2X3(C*2  +  C°3)  -  x^ 


1  - 


-o 

^1 


<°RT  ' 


-*  2 

+      o       *    ^^^1^ 


-    X- 


(2.53) 


_* 

V, 


1  - 


K3RT  ^J 
(2.54) 


(2.55) 


All  quantities  on  the  right  side  of  these  equations  can  be  determined 
from  binary  experimental  data. 

The  above  expressions  describe  the  thermodynamic  properties  of 
the  solvent  mixture  over  its  entire  composition  range,  i.e.,  from  pure 
2  to  pure  3.  The  component  activity  coefficients  are  expressed  in  the 
symmetric  convention,  i.e., 
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^2  "  ^2  "^  RT£n  X2Y2 


(2.56) 


o     * 
The  superscripts   and   again  refer  to  pure  2  and  pure  3,  respectively. 

Then 


RT   3x„ 


N^.T,? 


^Iny. 


-'2     ^ 


N3.T.P 


(2.57) 


Using   the   relation 


ay. 


r  3y, 


N, 


9N, 


3N, 
/ 

9x. 


9y, 


N2+N3 


N, 


N, 


(2.58) 


in  equation  (2.32)  leads  to 


N, 


3    3 


^  [1  +  C22  -  2x^0^2  -  2-3^3  ^  ^3(^23^23  -  ^22^3^^ 


(2.59) 


where 


2       2 
1  -  X2C22  -  ^3^33  -  2x^X30^3 


(2.60) 


At  the  compositions  of  2  and  3  at  infinite  dilution  in  each  other,  we 
get 


3  Any, 
~3x7 


x^-O 


and 


-e. 


3«.nY, 

~3x7 


X2=0 


C°  +  C°  -  2C°   -  (C°   C°  -   C°   C°   ) 
S3   ^22   ^^23    ^^22^33   ^23^23^ 

(1  -  C°2) 


*     *      *      **     ** 

S2  ^  ^33  "  ^^3  "  ^^2^33  ~  ^3^23^ 

* 
(1  -  C33) 


(2.61) 


(2.62) 


=  -e. 
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From  equation    (2.28) 


23  ,   -    .o 


^-^2 


~       =■    (PV3)  (2.63) 


•*•   ^23  "  ^  "  v°/<°RT  E  B2  (2.64) 

Similarly 

"^23  "  ^  "  ^2^^3^^  -  ^3  <2.65) 


Finally  from  equation  (2.42) 

C33  =  1  -  I/P3K3RT  =  a^  (2.66) 


*  * 


C°2  =  1  -  1/P2<2^T  =  a^  (2.67) 


Then  for  a  linear  composition  dependence  of  the  C   we  get 


C22  =  ^3X2(1  -  ^2)  +  x^(c,^   -  a^e^  -  (1  -  83)^)  (2.68) 

C33  =  a2X3(l  -  a3)  +  X2(a2  -  0263  -  (1  -  82)^)  (2.69) 

C23  =  X2B2  +  X3B3  (2.70) 

And,  substituting  in  equations  (2.29),  (2.28)  and  (2.59),  we  get 

2 

a  "  —^  =  1  -  x^a  (1  -  a,)  -  x^x-(a.  -  a.e.  -  T^T:  )  -  x^a,(l  -  ctj 


p<RT   "   "2^3'"   ^2'  "  ''2''3^°'3  ~  ^3^2  "  -^  "  1^3  -»  "  X2a2< 

2 


3' 


-  V2^^^2  -  "^2^3  -  1  -  ^2  )  -  2X2X3 (X3S3  +  X262)  (2.71) 

PV2  -  [1  -  X2a3(l  -  a^)  -  X2X3[a3  -  a^e^  -   6~^)^]   -  x^(x^&^+K^s^)]/a 

(2.72) 


26 


and 


dZny, 


9x 


2      J 


[1  +   (1  -   2x-){a„x„(l   -   a„)   -5-  x   [a_   -   a_e- 

X-        ax_  232  2  3      3  3   1 


-    (1  -  BJ]}-  2x-(x_6,  +  X-B-)   +  x_{(x.B^  +  x  3,)      -    [a^x    (1  -  a   ) 


3'   3  3  2'^2'        -^3^^-3-3  2^2 


3   2 


+  X    [a     -  a^e^  -    (1  -  62)    ]  ]  }la2^3*'-'-  "   ^3^   "^  x^Ca^  -  a^e^ 


-  1  -  B2  )}]] 


(2.73) 


These  expressions  are  used  in  Chapter  5   to  study   the   thermodynamic 
properties  of   liquid-liquid  systems. 
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CHAPTER  3 

PERTURBATION  THEORY  FOR  MOLECULAR  DISTRIBUTION 
FUNCTIONS  OF  DENSE  FLUIDS 


The  molecular  distribution  functions  in  a  real  fluid  are  quite 
complex  functions  of  its  molecular  interactions  and  the  macroscopic 
variables  of  the  system.  No  complete  theory  presently  exists  which 
can  accurately  predict  the  radial  distribution  function  (RDF)  in  a 
real  fluid  in  terms  of  its  intermolecular  potential  and  the  temperature 
and  bulk  density  of  the  system.  This  situation  is  even  more  complex, 
and  hopeless  from  the  point  of  view  of  prediction,  in  a  fluid  with 
many  chemical  species,  and  hence  different  molecular  interactions. 

The  problem  of  describing  the  RDF  and  relating  it  to  the 
intermolecular  potential  has  traditionally  been  approached  in  two 
distinct  ways.   The  first  method  is  the  postulation  of  theories  relating 
the  direct  correlation  function  (DCF) ,  the  RDF  and  intermolecular 
potential.   These  theories  are  postulated  or  derived  by  techniques  such 
as'  the  diagrammatic  methods  (HNC  theory)  or  by  functional  differentiation 
(Percus  Yevick) .   They  provide  a  relation  between  g(r),  c(r)  and 
(^(r),  where  g(r)  and  c(r)  depend  on  the  temperature  and  density  of  the 
system,  whereas  <})(r)  is  assumed  independent  of  the  macroscopic 

conditions.   These  theories  are  moderately  successful  in  describing 

12  3 
the  dense  liquid  state.  '  ' 

The  other  approach  to  the  problem  at  hand  is  perturbation 

theory.   In  this  method  the  real  fluid  is  modelled  as  a  perturbed  state 

about  a  reference  state  '  whose  properties  can  be  well  determined;  the 

perturbation  is  usually  in  the  intermolecular  potential,  and  the  RDF  is 
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determined  as  a  function  of  the  perturbation  parameters.   These  theories 
have  been  quite  successful  in  determination  of  thermodynamic  properties; 
however,  for  close  agreement,  second  order  terms  in  the  perturbation 
parameters  must  be  included. 

In  the  present  study  we  have  derived  a  simple,  but  rigorous, 
first  order  perturbation  theory  for  a  real  fluid.   For  computational 
convenience,  the  base  fluid  is  taken  to  be  the  hard  sphere  fluid, 
though  this  is  not  a  requirement  of  the  theory.   A  first  order  density 
expansion  is  obtained  relating  the  RDF's  and  DCF's  in  the  real  and  hard 
sphere  fluids.   Combined  with  the  Omstein-Zemike  equation  for  the  real 
fluid,  the  perturbation  expression  provides  a  scheme  for  determining 
the  real  fluid  distribution  functions.   To  reduce  the  solution  of  an 
integral  equation  to  that  of  computing  a  numerical  integration,  an 
intuitive  approximation  is  made  about  the  nature  of  the  real  fluid 
distribution  function.   The  validity  of  this  approximation  is  borne 
out  by  the  close  agreement  of  the  calculated  RDF  to  that  resulting  from 
computer  experiments.   The  DCF  of  the  real  fluid  is  then  calculated 
from  the  known  intermolecular  potential  and  calculated  RDF.   It  is 
also  numerically  integrated  to  yield   c(r)dr^,  which  has  earlier  been 
shown  to  be  related  to  the  isothermal  compressibility  of  the  pure  fluid. 
The  integral  is  computed  as  a  function  of  reduced  temperature  and 
density  giving  a  table  of  values  of  the  reduced  compressibility 
integral. 

In  multicomponent  systems,  the  DCF  and  its  integral  depends  on 
the  composition,  and  no  theory  is  yet  available  to  describe  this 
dependence  accurately.   Empirically,  the  composition  dependence  of  the 
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DCF  integral  is  represented  as  being  through  the  macroscopic  volume 
of  the  system  which  depends  on  the  composition.   Then,  using  represen- 
tative molecular  parameters  for  the  different  intermolecular 
attractions  present,  the  DCF  integrals  in  the  '.?.ixture  can  be  evaluated 
at  any  liquid  density.   This  method  is  used  to  determine  the 
thermodynamic  properties  in  several  mixtures  of  subcritical  components 
for  the  infinitely  dilute  composition  states.   The  agreement  between 
calculated  and  experimental  values  is  good  enough  to  classify  this 
theory  as  acceptable,  for  the  purposes  for  which  it  is  used.   Further 
work  in  this  regard  is  clearly  warranted. 

3,1  First-Order  Perturbation  Theory 

The  perturbation  expression  is  derived  by  the  method  of 
functional  differentiation  described  by  Percus;   It  is  based  on  the 
recognition  that  the  molecular  distribution  functions  are  functional 
derivatives  of  the  molecular  density  at  any  point  in  the  fluid  with 
respect  to  a  change  in  potential  on  the  fluid.   This  provides  a  precise 
relation  (hierarchy)  between  molecular  distribution  functions  of 
different  orders. 

Details  of  the  derivation  are  presented  in  Appendix  C.   The 
resulting  integral  equation  is: 

rf.  .       g''^(r)exD^^*'  ^""^^        ,    hs,  ,    hs,  , 
g^^(r)  -  1  +  ^ — LiiSiS 1  _  g  (r)  -  c  (r) 

hs,  ,   [B4>  (r)] 
g  (r)exp 


-g  <^)^^P   +  1  -  g^^s)  +  c^^(s) 


hs,  V 
g  (s) 


rf 


[g"(r')  -  1]   dr'  (3.1) 


3i 


where   s  =  ll.  -  £_'  '  .  Now  making  the  approximation  that  for  terms 
within  the  integral  we  can  write 


rf,  V    hs,  . 
g  (r)  -  g  (r)  (3.2) 


we  get 


rf  ,  .    ,rf .  . 
y  (r)  -  d   (r)  -  1  =  p 


[f^^s)yh^s)  +  d^^s)][e-^*'''<'^')yh^r')  -  l]dr' 


+  pj    [d^^s)  +l-e^*''^^^>][y'^^r')  -  1]   dr'     (3.3) 


a 

1 


-  1 


where  we  are  using  Rowlinson's  notation 

y(r)  =  g(r)e^*^'^^  ;   c(r)  =  f(r)y(r)  +  d(r)  ;  f(r)  =  e"^*^'^^ 

(3.4) 

Given  an  intermolecular  potential  function,  and  a  value  of  the  reduced 

* 

hard  sphere  diameter,  a  -  a/a,  where  a   is  the  distance  of  zero 

potential,  the  real  fluid  distribution  functions  can  be  calculated  for 
any  thermodynamic  state. 

Calculations  for  the  Lennard-Jones  6-12  potential  have  been 
made  at  two  thermodynamic  states  in  the  dense  liquid  region.   The 
hard  sphere  diameter  was  determined  by  the  Barker-Henderson^  prescription, 
and  the  hard  sphere  distribution  function  by  the  Walts  and  Henderson'^ 
modification  of  the  Percus-Yevick  equation.   The  resulting  RDF's  are 
shown  in  Fig.  3-1  and  compared  with  molecular  dynamics  values^  of 
Verlet.   For  one  of  the  two  states,  we  also  show  the  RDF  determined 

from  the  usual  first  order  Percus-Yevick  theory.   As  anticipated,  the 

rf 
calculated  g   (r)  follows  the  molecular  dynamics  result  closely  out  to 
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Fig.  3-1.  Radial  distribution  function  for  Lennard- Jones 
6-12  fluid. 
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large  reduced  distances  in  both  magnitude  and  phase,  and  is  clearly 

an  improvement  over  the  Percus-Yevick  result.   The  calculated  first 

rf 
peak  of  g   (r),  which  is  most  inportant  in  the  determination  of 

thermodynamic  properties  is  also  much  closer  to  the  molecular  dynamics 

peak  than  is  the  Percus-Yevick  peak.   This  indicates  that  thermodynamic 

properties  determined  from  the  g   (r)  calculated  from  (3.3)  should 

agree  more  closely  with  experimental  (computer)  results  than  do  those 

determined  from  the  Percus-Yevick  theory. 

rf 
To  obtain  improved  solutions  of  g   (r) ,  equation  (3.3)  could 

be  used  in  an  iterative  solution  scheme.   It  is  expected  that  this 

procedure  will  yield  results  in  even  closer  agreement  with  molecular 

dynamics  values.   In  the  same  fashion,  varying  the  hard  sphere 

diameter  in  ways  other  than  that  of  the  Barker-Henderson  prescription 

could  lead  to  improved  results. 

3..2  Reduced  DCF  Integrals 

The  agreement  obtained  between  the  calculated  RDF  and  Monte 
Carlo  results  indicates  that  the  perturbation  theory  developed  here  is 
a  satisfactory,  if  not  completely  accurate,  model  of  the  molecular 
interactions  in  a  Lennard-Jones  6-12  fluid.  We  now  postulate  that 
within  the  limits  of  accuracy,  the  Lennard-Jones  6-12  potential 
represents  the  interactions  in  a  real  fluid. 

In  spherically  symmetric  intermolecular  potentials,  the 
compressibility  equation  of  state,  or  equation  (2.42)  is 


-  r,-,.  =  1  -  ^TT 

pKcRT 


c(r)r^dr  (3.5) 
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In  terms  of  the  reduced  variables 


P  =  pa\  T  =  ^  ,   d  =f  ,  r  =f  (3.6) 


we  have 


e    (1  -  p  C)"-^ 

^     m     tr    — —  =  -^ ~ 

3     -   3  ~ 
a      p  d  T 


(3.7) 


C  =  AiTd^ 


c(f)  r^dr  (3.8) 


The  integral  C  has  been  calculated  for  state  points  covering  the  entire 
liquid  region,  using  the  perturbation  theory  developed  earlier. 
Fig.  3-2  shows  the  variation  of  this  quantity  with  reduced  density  at 
several  values  of  the  reduced  temperature  including  the  critical 
temperature  and  the  triple  point.   At  a  fixed  reduced  density,  this 
integral  is  not  a  strong  function  of  temperature,  and  in  the  dense 
liquid  region,  it  may  be  taken  as  approximately  constant.   It  does 
vary  strongly  with  reduced  density  along  an  isotherm.   This  is  not 
surprising  since  the  direct  correlation  function  c(r),  especially  at 
small  r,  is  a  strong  function  of  density. 

It  is  now  a  simple  matter  to  calculate  the  isothermal 
compressibility  of  the  fluid  at  a  given  thermodynamic  state.   All 
that  is  required  are  appropriate  values  of  the  potential  parameters 
a   and  e/k.   However,  before  establishing  a  general  rule  for  determining 
these  from  some  known  macroscopic  property  such  as  a  second  virial 
coefficient  or  shear  viscosity,  it  is  advantageous  to  examine  further 
macroscopic  implications  of  the  graphs  sho^.m  in  Fig.  3-2. 
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Fig.  3-2.   Reduced  direct  correlation  function  integral 

for  Lennard- Jones  6-12  fluid  from  perturbation 
theory. 
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The  reduced  compressibility  <    can  be  re-expressed  as 

,  ,  1 =  1  +  -  -  ^  (3.9) 

p<RT       p 

Since  we  know  from  macroscopic  evidence  that  <  depends  on  p  and  T, 

we  expect  that  1  +  —  -  ^  also  depends  on  p  and  T.   The  calculated 

p  <  T 

quantity  1  +  -^  -  -,  where  <  is  obtained  from  (3.7),  is  shown  in 
p  <  T 

Fig.  3-3.   Over  the  entire  liquid  range  it  does  not  vary  with  tempera- 
ture by  more  than  10%  at  a  fixed  reduced  density,  and  may  be 
considered  constant.   This  assumption  makes  isochoric  lines  on  a  P-T 
plot  linear. 

The  constancy  of  the  reduced  compressibility  over  the 
temperature  range  provides  the  microscopic  verification  for  a  macroscopic 
correlation  of  the  isothermal  compressibility  of  pure  fluids.   The 
development  of  the  macroscopic  correlation  will  be  discussed  in  detail 
later.   The  general  result  for  pure  liquids  is  that  (1  +  —  „„)  is  a 
universal  function  of  a  reduced  density  and  is  independent  of  temperature. 
The  reducing  density  is  close  to  the  critical  density  for  nonpolar 
liquids,  and  larger  than  the  critical  for  polar  liquids.   These 
correlating  densities  are  used  as  the  source  of  molecular  parameters 
in  the  microscopic  theory.   Existence  of  the  correlation  asserts  that 
macroscopic  corresponding  states  is  valid  for  this  property,  and  the 
representation  of  all  molecular  interactions  is  equivalent  to  postulating 
corresponding  states  on  the  microscopic  level.   Then  taking  Argon  as 
the  prototype  or  reference  substance,  i.e.,  as  satisfying  both 
macroscopic  and  microscopic  corresponding  states,  the  molecular 
parameters  for  the  other  fluids  are  determined  from  the  relation 
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Fig.  3-3.  Microscopic  correlation  of  isothermal  compressibility 
for  Lennard-Jones  6-12  fluid. 
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(3.10) 


Ar 


* 
where  v  is  the  reciprocal  reducing  density  for  the  macroscopic 

correlation.   Values  of  c  for  some  substances  are  shown  in  Table  3-1. 


TABLE  3-1 
INTERMOLECULAR  PARAMETERS  FOR  SIMPLE  FLUIDS 


Substance  -»- 

Ar       CH^ 

^2 

O2 

CO 

S\ 

* 

V  cc/mole  -'■ 

74.57    99.5 

90.1 

73.4 

93.1 

127.3 

a  A       -*■ 

3.405^   3.749 

3.627 

3.387 

3.666 

4.069 

H2O      CCl^ 

^6»6 

SO2 

S«8 

^^4^10 

46.4    276.0     255.0     115.0     200.0     255.0 
2.907    5.267     5.130     3.934     4.752     5.130 


CS2 

CH  OH 

NH^ 

Cyclohexane 

165.0 

101.5 

65.18 

311.0 

4.437 

3.774 

3.756 

5.481 

+. 


From  second  virial  coefficient  data 


The  value  of  a,   is  taken  from  second  virial  coefficient  data. 
Ar 

For  other  substances,  the  value  calculated  from  (3.10)  is  generally 
smaller  than  that  obtained  from  second  virial  coefficients.  The 
highly  polar  substances,  water  and  methanol,  have  relatively  much 
smaller  values  of  a. 


The  function  C  is  assumed  terp-ri"  r:i  Independent  it  a  fixed 

reduced  density  and  no  temperature  related  -noaecular  par;iae.;er  is 

necessary  to  characterize  the  microscopic  s-ace  of  the  fluid.  C  at 
a  given  p  is  determined  from  Fig.  3-1. 

2^2    Multicomponent  Systems 

In  a  mixture  the  quantities  C. .  depend  on  the  molecular 
interactions  and  the  macroscopic  state  of  the  system.   No  present  day 
molecular  theory  can  describe  the  composition  dependence  accurately, 
(The  multicomponent  perturbation  theory  of  Barker-Henderson-Smith 
is  a  first  attempt  at  describing  the  composition  dependence  of  the 
multicomponent  radial  distribution  function.) 

An  intuitive  method  is  proposed  here  to  represent  the  dependence 
of  C   on  the  system  variables.  It  is  based  on  two  postulates,  which 
are: 

1.  All  molecular  interactions  in  the  mixture,  <i> ,  , ,    are  of 
the  same  type;  i.e.,  they  are  representable  by  a  unique  reduced 
potential  function,  the  Lennard- Jones  6-12  potential. 

2.  The  composition  dependence  of  the  C  ,  is  completely 
described  by  the  macroscopic  density  of  the  system,  which  is 
composition  dependent.  This  assumption  neglects  the  effects  on  C  . 

of  composition  changes  at  constant  density  but  emphasizes  that  the  most 
significant  parameter  for  determining  C  .  is  the  macroscopic  density. 

On  the  basis  of  these  two  postulates,  the  general  reduced 
integral  C. .  may  be  written  as 
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where  a.,  is  the  size  parameter  for  the  i-j  ir.tarajcicn.   Here,  it  is 
calculated  from  the  arithmetic  mean  rule 

^i3  =  i  ^^ii  ^ '^jj)  ^'-^^^ 

These  equations  form  a  method  for  determin?lng  tha  integrals  C.^.  ,  C .  . 
and  C..  at  any  given  liquid  density  from  which  the  tner-.r.odynamic 
properties  can  be  calculated.   The  calculated  thermodynanilc  prctjercies 
for  some  binary  systems  are  shown  in  Table  3-2  and  cor^pared  with 
experimental  values  as  available.   The  composition  dependant  integrals 
C, .  have  been  calculated  only  at  the  state  of  infinite  dilution  of 
component  1  in  component  2.  For  all  systems  of  simple  molecules 
Ar,  0-  and  N„ ,  there  is  good  agreement  between  calculated  and 
experimental  isothermal  compressibilities  and  partial  molar  volumes  at 
infinite  dilution.  This  indicates  that  the  molecular  representation 
assumed  by  the  theory  is  quite  valid  for  these  systems.   For  the  other 
systems  of  more  complex  molecules,  the  agreement  is  only  qualitative, 
although  the  trends  of  all  calculated  quantities  are  in  the  direction 
of  the  experimental  values. 

Comparison  of  the  calculated  C. .  with  the  experimental  values, 
where  these  are  available,  indicates  that  the  calculated  C.  are 
not  negative  enough,  e.g.,  C^^  in  the  n-pentane  systems.   The  larger 
negative  values  necessary  are  obtained  only  by  calculating  the  quantities 
at  larger  reduced  densities  (i.e.,  the  a   values  are  too  sm.all) ,  or 

by  introducing  a  temperature  dependent  parameter  into  the  formulation 

o 
of  the  theory,  which  would  then  allow  somewhat  larger  negative  C  , . 

However,  it  is  more  likely  that  the  assumptions  of  conformal  solutions 
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(potentials)  and  pairwise  additivity  break  down  for  complex  molecules 
in  dense  systems. 

Nevertheless,  the  above  theory  provides  a  basis  for  che 
qualitative,  and  often  quantitative,  understanding  of  the  state 
dependence  of  the  direct  correlation  functions  and  their  integrals. 
Two  valid  conclusions  can  be  drawn: 

1.  The  calculated  integrals  C   ,  C^j  and  C  „  are  approximately 
in  the  correct  ratio  to  one  another. 

2.  When  the  calculated  integral  C^^   is  within  5-10%  of  the 
experimental  value,  the  calculated  C^ „  and  C^ .  are  also  within  this 
range  of  the  corresponding  experimental  quantity.   Wlien  the  theory  is 
inadequate  in  predicting  experimental  data,  it  is  due  to  the  poor 
modelling  of  the  pure  solvent  state. 
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CHAPTER  4 

MACROSCOPIC  STATE  DEPENDENCE  OF  DIRECT 
CORRELATION  FUNCTION  INTEGRALS 


In  the  preceding  chapter  we  studied  a  molecular  theory  for  the 
state  dependence  of  the  DCF  integrals.  Although  the  generalized 
microscopic  method  gave  some  adequate  results  for  thermodynamic 
properties,  there  are  still  areas  of  uncertainty  in  the  molecular 
approach,  and  no  real  clues  are  provided  for  improvements  to  the  molecular 
theory.   In  this  chapter  we  study  the  available  macroscopic  data  to 
determine  the  state  dependence  of  the  DCF  integrals. 

A  definite  hierarchy  exists  in  the  calculation  of  0-2 >  C  „ 
and  C^.  from  experimental  data,  and  this  order  is  followed  in  searching 
for  macroscopic  correlations.   Generalized  equations  for  the  temperature 
variation  of  the  pure  solvent  saturation  pressure  and  liquid  volumes 
are  combined  with  the  previously  derived  correlations  to  derive  the 
temperature  and  pressure  dependence  of  the  C  . 

4.1  Corresponding  States  Correlations  for  C„„  and  C^ „ 

The  DCF  integral  C-.  is  related  to  the  isothermal  compressibility, 
K,  o£  a  pure  fluid  by  the  equation 

It  is  experimentally  derived  from  P-V-T  measurements  of  compression 
experiments.   Commonly,  the  compression  data  are  filled  to  analytic  and  fitted 
expressions  such  as  the  Tait  Equation  (for  most  liquids)  and  the 
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Kuddleston  equation  (liquid  hydrocarbons).^  Experimentally,  <   for  a 
pure  liquid  is  found  to  be  a  strong  function  of  temperature,  and  very 
weakly  dependent  on  pressure.   As  the  temperature  approaches  the 
critical,  the  isothermal  compressibility  of  the  saturated  liquid 
becones  very  large. 

In  Fig.  4-1,  the  experimental  quantity  1  +  1/pkRT  is  shown  as 
a  function  of  p  for  several  simple  liquids;  the  plot  for  n-decane  is 
along  an  isotherm  (358°K) ,  i.e.,  a  compressed  liquid,  while  the 
others  are  along_the_U_quid_vapor  saturation  line.   The  curves  are  all 
similar  in  shape  and  one  concludes  that  there  exists  a  general  correla- 
tion of  the  form 


1-f  1 


PkRT 


=  F(p)  =  F(pv  )  (4^2) 


Experimental  compressibility  data  for  Ar,  0.,  N  ,  CH, ,  CCl, ,  nC  H 

z.       I  4     4'   10  22* 

"^12"26'  "^IaSo  ""^-^   ''^16"34  ^"^^  "^^'^  *=°  derive  an  analytical 
expression  for  this  general  relationship.  The  first  five  of  these 
substances  are  chosen  among  the  primary  set  because  they  are  simple 
molecules.   Further,  their  characteristic  volumes  can  be  taken  as 

the  critical  volumes  to  yield  a  set  of  curves  for  fl  +  ^   1  that 

\  p<RT  J  ^"" 

are  coincident  within  experimental  error.   The  function  fl  +  ^^ 

t    PkRT 
then  has  a  limiting  value  of  unity  as  p~  -  pv*  ->  1;  however,  detailed 

experimental  data  are  not  available  in  this  lower  reduced  density  range, 

The  reduced  density  range  for  this  set  of  substances  extends  over 

1.5  £  p  <  2.9.   The  long  chain  hydrocarbons  are  included  in  the  primary 

set  to  extend  the  reduced  density  range  of  the  correlation  to  larger 
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Fig.  4-1.   Isothermal  compressibilities  of  pure  liquids. 
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values;  their  compressibilities  are  reported  to  be  of  a  high  degree 
of  accuracy  (1  1%).  The  characteristic  volumes  for  these  hydrocarbons 
are  obtained  by  matching  the  experimental  data  to  the  curve  for  simple 
substanceslin  the  region  of  overlap.  The  algebraic  expression  obtained 
to  represent  the  experir.eatal  data  over  the  range  1.5<_p<_3.7  is 

2  -  C°2  =  1  +  —^  =  e>:p[-. 42704 (P  -  1)  +  2.089(0  -  1)^  -  .42367(P  -  1)^] 

(A. 3) 

Fig.  4-2  is  the  ^i^sphical  representation  of  equation  (4.3)  and 
shows  some  representative  experimental  data  to  indicate  the  average 
error  of  prediction.   For  all  substances  other  than  those  listed  above, 
the  characteristic  voluz:es  were  obtained  by  fitting  experimental 
compressibility  data  to  equation  (4.3).   For  nonpolar  substances,  the 
fitted  characteristic  volumes  are  estimated  to  be  within  the  errors 
of  experimental  and  predictive  methods  for  the  critical  molar  volume. 
Thus,  In  the  absence  of  experimental  compressibility  data,  the 
critical  volume  of  a  nonpolar  substance  or  a  weakly  polar  substance 
may  be  used  as  the  characteristic  volume.   For  polar  substances, 
however,  e.g.,  water  and  methanol,  v  as  fitted  from  equation  (4.3) 
is  found  to  be  less  than  v  .   (On  the  molecular  scale,  v  is  related 
to  the  distance  parameter  in  the  angle  averaged  or  effective 

spherical  potential.)   There  is  no  clear  and  lucid  explanation  of  this 

* 

relationship,  however,  that  would  permit  the  determination  of  v 

from  other  macroscopic  data.   Therefore,  this  parameter  must  be 
obtained  from  at  least  one  experimental  point. 
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Fig.  4-2.  Generalized  correlation  for  isothermal 
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The  results  of  compressibility  calculations  with  equation  (A. 3) 

are  shown  in  Table  4-1  for  nine  representative  substances.   5,  the 

fractional  root-mean-square  deviation  in  1  +  — —   ,  is  always  less 

than  0.08;  the  largest  single  error  in  1  +  — r^r-  encountered  over  an 

p  <Ki 

aggregate  of  70  state  points  is  11%.   Characteristic  volumes  for  all 
the  substances  included  in  the  correlation  are  presented  in  Table  4-2; 
the  list  contains  liquids  of  all  types,  polar,  nonpolar,  long  chain 
and  branched  hydrocarbons,  alcohols,  esters,  inorganic  liquids,  etc. 
In  general,  for  substances  with  known  characteristic  parameters,  both 
polar  and  nonpolar,  equation  (4.3)  provides  a  value  of  <  with  an 
average  deviation  of  6-8%  at  p  greater  than  1.5. 

In  a  binary  solution,  C  „  is  related  to  the  partial  molar 
volume  of  component  1  at  infinite  dilution  in  component  2.   (When 
component  1  is  a  gas ,  whose  critical  temperature  is  appreciably  lower 
than  the  temperature  of  the  solution,  the  infinitely  dilute  state  is  a 
thermodynamically  conceptualized  state.)   The  partial  molar  volume  is 
determined  experimentally  by  observing  the  change  in  volume  of  the 
mixture  as  the  component  composition  is  changed.   It  is  then  extracted 
from  the  experimental  volume  by  the  general  procedure  for  partial  molar 
properties  and  extrapolated  to  zero  composition.   If  the  solution  is 
assumed  to  be  ideal,  the  partial  molar  volume  at  infinite  dilution 
(pmv  aid)  can  be  obtained  from  the  change  in  solute  solubility  with 
pressure.   Eqn.  (2.45)  is 

pv°  =  (1  -  C°2)/(l  -  C°2)  (^.^) 

.-.    v°  =  v°(l  -  C°2)/(l  -  C°2)  (^.5) 

C^2  ^^s  calculated  from  equations  (4.5)  and  (4.3)  with  experimental 
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TABLE  4-2- 

-a 

REDUCING  VOLUMES 

FROM  ISOTHERMAL  COMPRESSIBILITIES^ 

RMS 

* 

Temp.  : 

Range 

#  Data 

Dev. 

V 

Data 

Substance 

°K 

Points 
11 

% 

cc/g  mole 
276^ 

Reference 

CCl^ 

250  - 

343 

4.7 

1 

"S2«26 

298  - 

358 

5 

6.0 

730 

2 

^^14^30 

318  - 

358 

4 

5.0 

845 

2 

-S6»34 

298  - 

358 

4 

5.8 

970 

200^ 
255^ 

2 

S«8 

310  - 

327 

6 

5.5 

6 

^^4^10 

310  - 

344 

6 

8.6 

7 

-S^2 

273  - 

344 

5 

6.2 

309 
308^ 
369^ 
425^ 

3 

iso  C3H^2 

223  - 

298 

4 

3.5 

11 

-^6^14 

273  - 

444 

6 

6.1 

3 

-S«16 

273  - 

473 

7 

3.1 

1 

°S«18 

403  - 

553 

6 

3.2 

489 
541^ 

3 

^S»20 

303  - 

473 

5 

2.9 

3 

"^11^24 

298  - 

473 

3 

5.75 

670 
783^ 

3 

°^13«28 

323  - 

473 

3 

2.9 

3 

°S5»32 

333  - 

408 

3 

2.5 

915 
1035^ 
llOO*' 

4 

-^17^6 

323  - 

523 

3 

3.6 

3 

^^18^38 

333  - 

408 

3 

2.58 

4 

°S0»42 

323  - 

573 

6 

3.4 

1225 

3 

°So^62 

373  - 

573 

3 

2.82 

1880 

3 

^6«5^3 

223  - 

298 

4 

3.7 

312 

11 

S«5^°2 

298  - 

358 

4 

3.3 

321 

8 

^6«5^2 

298  - 

358 

4 

2.78 

285 

8 

C^H3C1 

298  - 

358 

4 

1.0 

306 

8 

C^H3Br 

298  - 

358 

4 

2.68 

321 

8 

0- xylene 

298 

2 

3.3 

363.5 

9 

m-xylene 

298 

2 

2.6 

362.5 

9 

Cumene 

298 

1 

2 

2.1 

298 

9 

Also  for  use  in  partial  molar  volume  correlation 
Within  .5%  of  experimental  or  predicted  critical  volume. 
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TABLE  4-2-a  (Continued) 


Substance 

Temp.  Range 
°K 

Ir   Data 
Points 

2 

EMS 

Dev. 

% 

1.9 

* 

V 

cc/g  mole 
423 

Data 
Reference 

Mesitylene 

298 

9 

Tetrahydronaph- 
thalene      298 

2 

3.0 

430.2 

9 

CH.Cl^ 

303  - 

323 

2 

.06 

168.8 

10 

C2H5OOCH3 

303  - 

323 

2 

7.5 

283 

10 

CHCI3 

298 

2 

.55 

219.6 

9 

Cyclohexane 

298  - 

348 

4 

1.21 

311 

12 

Isopropyl  el 

ther  323 

1 

- 

415.8 

10 

Tetrahydro- 
furan 

303  - 

323 

2 

6.9 

237 

10 

Ethyl  Ether 

273  - 

413 

8 

6.5 

281.2 

15 

(CHCl2)2 

298 

-■ 

286.7 

13 

(CCl2)2 

298 

- 

303.6 

13 

(CH^Cl)^ 

298 

- 

226 

13 

CHCl  CCI2 

298 

- 

258.9 

13 

(CHCl) 2 

298 

- 

208.6 

13 

CH2CHCI2 

298 

233 

13 

CH  COCH 

298 

• 

200.2 

13 

PSU  25^ 

310  - 

408 

4 

3.6 

1550 

14 

PSU  19^^ 

310  - 

408 

4 

2.7 

1465 

14 

PSU  174^^ 

352  - 

408 

4 

3.7 

1332 

14 

PSU  528'^ 

333  - 

408 

4 

8.3 

785 

14 

PSU  532^^ 

333  - 

408 

5 

5.9 

933 

14 

PSU  537^^ 

333  - 

408 

5 

5.1 

1116 

14 

PSU  175^^ 

333  - 

408 

3 

1.1 

1452 

14 

PSU  12^ 

333  - 

408 

3 

1.8 

1368 

14 

PSU  574^^ 

333  - 

408 

3 

3.7 

800 

14 

PSU  575^^ 

333  - 

408 

3 

5.4 

873 

14 

PSU  122^ 

373  - 

408 

3 

.75 

1708 

14 

CS2 

298  - 

358 

4 

2.0 

165 

14 

'Complex  hydrocarbon.   See  reference  14  for  formula. 


TABLE  4-2-a  (Continued) 


RMS 

a 

TsiriD .  Range 

# 

Data 

Dev. 

V 

Data 

Substance 

°K 

Pj 

?ints 
3 

% 
2.7 

cc/g  mole 
101.5 

Reference 

Methanol 

273  -  323 

9 

Propanol 

298 

4 

2.5 

160.3 

9 

NH. 

253  -  313 

16 

S.8 

65.18 

17 
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TABLE  4-2-b 
CHARACTERISTIC  VOLU>tES  DETERMINED  FROM  PARTIAL 


MOLAR  VOLUMES  AT  lOTINITE  DILUTION 


.T^.b 


Substance 

//  Data  Points 

#  Solvents 

V  cc/g  mole 

CO 

10 

3 

93.1^ 

co^ 

6 

2 

80 

«2 

19 

5 

51.5^ 

^2^6 

4 

4 

158 

CF, 

2 

2 

139 

S«4     • 

3 

3 

127 

^2«2 

3 

3 

112.6^ 

^^^3^2° 

3 

3 

169.7 

CH^Cl 

3 

3 

136.5 

SO, 

3 

3 

115 

Also  for  use  in  compressibility  correlation. 
Sources  of  experimental  data  quoted  in  Table  4-3. 
"Assigned  to  be  critical  volume. 


Effective  critical  volume.   See  reference  28. 
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data  on  a  wide  variety  of  solute  (1)  and  solvent  (2)  pairs.   In  cases 
where  the  total  experimental  pressure  had  not  been  reported,  the 
volume  of  the  pure  solvent  was  taken  as  that  of  the  saturated  liquid. 
Hypothesizing  a  general  dependence  of  the  DCF  integrals  on  the  reduced 
solvent  density,  the  variation  of  C^  was  found  to  be  similar,  but  not 
identical,  for  all  solute-solvent  pairs.   A  factor  accounting  for  the 
relative  molecular  sizes,  as  represented  by  the  ratio  of  the 
characteristic  volumes,  was  introduced.   The  empirical  relationship 
obtained  for  C  „  is  shown  in  Fig.  4-3,  with  some  representative  systems. 
The  analytic  expression  for  this  relationship  is 
.62  \ 


_* 

V, 


=  -2.4467  +  2.12074  p°  ;2.0  1  P2  1  2.785 

=  3.02214  -  1.87085  p°  +  .71995  P2^;2.785  1  P^  1  3-2 


(4.6) 


Calculation  of  v  for  a  solute-solvent  pair  of  known  characteristic 
parameters  thus  requires  only  a  knowledge  of  v„ ,  the  pure  solvent 
volume.   C._  is  obtained  from  (4.6),  C-o  from  (4.3)  and  v^^  from  (4.5). 
p_  is  the  reduced  solvent  density  and  is  used  in  both  correlations. 

The  numerical  coefficients  in  equation  (4.6)  were  generated 
using  V-  for  seven  gases  in  different  solvents.   For  Ar,  N-,  0^  and 
CH,  the  characteristic  parsimeters  had  been  determined  earlier  from 
compressibility  measurements.  For  CO  and  C„H  ,  the  choice  of  the 
critical  volune  as  the  characteristic  parameter  proved  satisfactory; 
similarly  the  effective  critical  volume  was  chosen  for  H^. 

Equation  (4.6)  provides  a  satisfactory  representation  of 
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Fig.  4-3.   Generalized  correlation  for  partial  molar 
volumes  of  gases  at  infinite  dilution  in 
liquids. 
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experimental  partial  molar  volumes.   Since  the  experimental  variable 
in  the  correlation  is  the  solvent  density,  which  is  a  function  of  both 
temperature  and  pressure,  the  variation  of  the  volumetric  properties  with 
pressure  is  partially  accounted  for  by  equation  (4.6).   Partial  molar 
volumes  for  varied  gas-solvent  systems  as  calculated  from  equation  (4.6) 

are  shown  in  Table  4-3,  and  compared  with  experimental  data.   Also  shown 

29 
are  volumes  calculated  from  the  generalized  correlation  of  Lyckman, 

which  is  based  on  a  cohesive  energy  density  for  the  solvent.   The 

accuracy  of  the  two  correlations  is  similar  for  saturated  systems  at 

low  temperatures.   However,  since  the  correlation  of  equation  (4.6) 

is  based  on  the  reduced  solvent  density  which  is  a  function  of  both 

temperature  and  pressure,  it  can  account  for  the  pressure  effect  on 

volumetric  properties  to  some  extent,  and  for  the  rapid  change  in  solvent 

volume  at  temperatures  approaching  the  critical.  This  is  clearly  seen 

in  the  hydrogen-benzene  system  where  a  30%  change  in  partial  molar 

volume  with  pressure  is  well  predicted,  and  the  high  temperature  value  is 

correlated  closely. 

The  characteristic  parameters  for  several  substances  other  than 

the  seven  previously  listed,  were  obtained  by  fitting  experimental 

partial  molar  volumes  to  equation  (4.6).   In  each  case,  v^  for  the  gas 

in  several  different  solvents  was  used.  .  They  are  shown  in  Table  4-2-b. 

4.2  Generalized  Isothermal  Equation  of  State  for  Liquids 

Macroscopic  experimental  evidence,  equation  (4.3),  shows  that 

the  function  1  +  — r^z-  is  a  universal  function  of  reduced  density  for  all 
picRT 

liquids.   Substituting  for  pic  from  equation  (4.3)  into  equation  (1.15) 


TABLE  A- 3 

PARTIAL  MOLAR  VOLLI'lES  OF  GASES  AT  INFINITE  DILUTION 

o 


V  cc/gm  mole 


Data 


Solute   Solvent   T^^K     r.  Calculated   Experiment  Lyckman  Reference 

CH,      nC,       300    2.86       56.6  55.4       58.2      17 
4        7 

CF,      CCl,      300    2.845      69.1  79.7       66.3 
4        4 

CF,      nC,       300    2.86       88.7  86.4       73.7 
4        7 

H-       C,H       298    2.865      31.1  35.0       37.1      18 

■^  °  °      403    2.496      57.4  58.3       52.7 

2.521      53.9  53.7 

2.55       51.4  49.7 

433    2.375      66.9  72.1       60.6 

2.412      63.0  64.0 

2.443      56.7  57.5 

473    2.197      92.9  104.4       64.7 

2.251      81.6  85.6 

2.296      69.0  72.5 

H       nC        403    2.646      68.4  63.6       59.8 

^        "  2.683      62.8  58.0 

2.716      57.55  53.3 

473    2.355     104.6  110.3       67.9 

2.422      86.2  87.6 

2.474      77.1  72.7 

298    2.776      59.1  60.2       62.5      19 

76.8  67.9       72. 

88.8  93        93.2 

102.5  104.9  107.3 

50.4  51.5       53.4 

55.8  62.8       59.1 

49.7  55.7       52.7 
298    2.838      33.6  38         39.2      20 

41.8  44        43.5 
298    2.838      49.2  52.4       51.2      21 

64.7  65.9       67.1 

43.1  53.0       48.3 

57.1  60.8       61.7      22 

53.0  53.6       57.3 

67.9  66.7       63.2 


CH, 

nCg 

S«6 

S«8 

t^c, 

Ar 

«2 

°2 

h 

CCl, 
4 

Ar 

CH4 

CCl, 
4 

S«6 

^2 

^2% 

^2«2 

(CH3)20 

59 


TABLE  4-3 

(Continued) 

. 

Solvent 

T,°] 

cC     p 

0 

cc/gn  cole 

Solute 

Calculated 

Experinent 
54 

Lvcknan 

Data 

CH^Cl 

59.4 

58.8 

^t-^cttLllCe 

SO2 

53.6 

53.2 

61.4 

S"6 

C^H^Cl 

298 

3.033 

57.4 

63.4 

63.3 

22 

S"a 

54.6 

57.8 

58.0 

S"2 

47.0 

49.3 

53.8 

(CH3)20 

60.1 

65.3 

73.0 

CH^Cl 

52.9 

52.8 

63.2 

so^ 

47.7 

47 

57.6 

s^ 

^6«6 

298 

2.865 

60.1 

66.1 

64.7 

22 

^2»A 

56.9 

60.4 

59.3 

S«2 

49.4 

50.3 

54.9 

(CH3)20 

62.8 

67.1 

74.3 

CH^Cl 

55.2 

54.1 

64.6 

so^ 

49.8 

46.8 

58.9 

CH, 

43.2 

52.0 

48.8 

21 

^2 

• 

43.1 

52.6 

45.9 

CH, 

"So 

323 

2.998 

58.4 

53.0 

59.1 

23 

273 

3.162 

52.6 

49.0 

_ 

• 

248 

3.242 

46.1 

46.0 

_ 

CH, 

"^8 

223 

3.257 

40.8 

43.0 

52.6 

23 

2A8 

3.172 

48.4 

47.0 

58.0 

273 

2.726 

75.2 

59.0 

87.6 

CH, 

nC^ 

233 

3.156 

46.2 

44.0 

55.2 

23 

CH, 

nC^ 

223 

3.094 

44.5 

45.5 

54.1 

23 

248 

3.00 

47.9 

48.0 

59.8 

298 

2.812 

56.2 

53.5 

73.3 

«2 

CH^OH 

273 

2.566 

35.1 

35. 

28.2 

24 

^2 

• 

48.6 

47. 

25.3 

CO2 

• 

45.4 

40. 

39.5 

CO 

49.6 

47. 

37.2 

°2 

43.2 

42. 

31.7 

6C 


TABLE  4-3 

(Continued) 

Solvent 

T,°K 

P 

^ 

cc/gm  luole 

Data 

Solute 

Calculated 
51.6 

Experiment 
47. 

Lvckman 
39.4 

Reference 

CH4 

h 

CH  OH 

298 

2.500 

39.1 

35. 

30.4 

24 

^2 

54.0 

52. 

26.8 

CO2 

50.5 

43. 

41.8 

CO 

55.1 

51. 

39.5 

°2 

48.0 

45. 

34.0 

CH4 

57.3 

52. 

41.7 

»2 

CH^OH 

323 

2.428 

44.2 

42. 

32.9 

24 

^2 

60.7 

54. 

28.7 

CO2 

56.7 

46. 

44.5 

CO 

61.9 

56. 

42.2 

02 

54.0 

50. 

36.6 

CH^ 

64.4 

55. 

44.4 

h 

H^O 

273 

2.577 

24.9 

24. 

21.2 

24 

^2 

34.8 

41. 

29.4 

CO2 

32.4 

32. 

32.2 

CO 

35.5 

37. 

30.0 

^2 

30.8 

31. 

24.6 

CH^. 

37.0 

36. 

32.2 

«2 

H^O 

298 

2.569 

25.2 

26. 

21.8 

24 

^2 

35.2 

40. 

30.1 

CO2 

32.8 

33. 

32.9 

CO 

36.0 

36. 

30.7 

°2 

31.2 

31. 

25.3 

CH4 

37.4 

37. 

32.9 

«2 

h° 

323 

2.547 

26.1 

34. 

22.7 

24 

^2 

36.4 

38. 

30.9 

CO2 

34.0 

33. 

33.8 

CO 

37.2 

32. 

31.5 

O2 

32.3 

32. 

26.1 

CH, 

38.7 

38. 

33.7 
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TABLE  4-3    (Continued) 


o 

^1 

cc/gm  mole 

Data 

Solute 

Solvent 
""^S^IS 

T,°K 
426 

P 

Calculated 
109.3 

Experiment 
99.6 

Lyckman 
80.5 

Reference 

CO 

2.501 

27 

455 

2.374 

138.4 

128.8 

88.2 

483 

2.218 

189.9 

175.0 

94.9 

512 

2.029 

292.0 

267.3 

99.6 

N? 

Ar 

83.78 

2.671 

40.7 

35.4 

35.6 

25 

N? 

°2 

77 

2.765 

30.1 

35.6 

31.1 

25 

Ar 

°2 

83.78 

2.689 

32.0 

28.7 

31.0 

25 

90 

2.617 

35.7 

29.5 

34.2 

CH, 

Cyclo- 

298 

2.858 

37.0 

53.0 

53.3 

26 

hexane 

310.9 

2.819 

52.3 

57.0 

55.5 

344 

2.697 

65.2 

65.0 

62.7 

^2 

NH^ 

273 

2.448 

48.8 

43.6 

28 
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^2         ^2 
AP  --  I   ^  =  p*  I   dpRT[f(p)  -  1]  (4.7) 

P-,  Pi 

X  1 

^2 


RT 


V   = 


[f(o)  -  l]dp  (4.8) 


J 


'l 


where 


f(p)  =  exp[-.4270A(p  -  1)  +  2.019(p  -  1)^  -  .42367(p  -  1)^]      (4.9) 


Equation  (4.3)  relates  the  isothermal  pressure  change  corresponding  to 
a  density  change,  and  holds  over  the  same  reduced  density  range  as 
(4.3).   For  calculational  convenience,  the  definite  integral  is 
re-expressed  as 


P2  P 


<^2  -  ^>  W  = 


[f(p)  -  l]dp  - 
P^  P 


[f(p)  -  l]dp  (4.10) 

o 


where  p   is  an  arbitrary  base  density.   The  integral  from  p   to  p 
has  been  numerically  evaluated,  and  is  shown  in  Table  4-4.   (The 
integrals  are  tabulated  at  integrals  of  Ap  =  0.1;  they  were  evaluated 
for  Ap  =  0.002.) 

Equation  (4.8)  can  be  used  to  determine  unknown  densities  or 
pressures.   Given  P  ,  p   and  p_,  P„  is  found  easily  from  the  difference 
of  the  integrals  of  equation  (4.10).   Given  P^ ,  P-  and  p^ ,  one  looks 
for  a  reduced  density  p„  such  that  the  corresponding  integral  of  Table  4-4 
satisfies  equation  (4.10). 

The  characteristic  parameters  used  in  the  following  calculations 
were  those  obtained  from  the  experimental  compressibility  data. 
Considerable  manipulation  of  the  experimentally  measured  volumes  is 
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TABLE  4-4 
REDUCED  INTEGRALS  FOR  ISOTHER^IAL  EQUATION  OF  STATE 

P  1.4       1.5       1.6       1.7       1.8       1.9 

P 

J  [f(p-l)]dp      0         .0214     .0603     .1237     .2209     .3641 


P  2.0 

P 


[f(p-l)]dp       .5705 
P^ 


PQ 


2.1      2.2      2.3      2.4      2.5 
.8631    1.2739    1.8469    2.6423    3.7427 


^  2.6      2.7      2.8      2.9      3.0      3.1 

P 

[f(p-l)]dp      5.2596    7.3427   10.1895   14.0564   19.2680   26.2246 


P  ■  3-2      3.3      3.4      3.5      3.6 

P 

[f(p-l)]dp     35.4028   47.3454   62.6364   81.8564  105.5160 


Po 
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necessary  before  compressibilities  are  extracted,  e.g.,  fitting  to  an 
equation  of  state  such  as  the  Tait  equation.  Thus  the  use  of  these 
characteristic  parameters  in  the  integrated  relation  to  predict 
volumes  constitutes  a  severe  test  of  both  the  accuracy  of  the 
parameter  and  the  validity  of  the  generalized  equation  of  state.   Also, 
the  integral  in  Table  4-4  increases  very  rapidly  at  high  reduced 
densities. 

The  results  of  compression  calculation  on  liquids  with  equation 
(4.8)  are  reported  in  Tables  4-5  and  4-6.   In  most  cases,  _the_inij:ial 
state  corresponds  to  the  saturated  liquid.   Then  given  either  the  volume 
or  pressure  of  the  compressed  state,  the  other  is  calculated.   Where 
Tait  equation  parameters  were  available,  calculations  were  made  using 
this  equation  too.   The  Tait  equation  is  a  two-parameter  equation 
usually  written  in  the  form 

'  B  +  P 


V   «  V 

^2    1 


1  -  C  £n 


2 


B  +  P 


(4.11) 


The  constant  C  generally  has  a  value  of  .11  for  all  liquids,  whereas 
B  varies  with  temperature. 

Volumes  calculated  from  equation  (4.8)  agree  to  within  1%  for 
all  the  substances  considered,  over  the  entire  pressure  and  temperature 
range.   The  Tait  equation  with  one  fixed  constant  C  (as  for  IJH-)  gives 
results  which  are  slightly  improved  over  those  of  equation  (4.8). 
However,  when  both  Tait  contants  B  and  C  are  varied  with  temperature 
and  fitted  to  experimental  data  as  for  nC^-H„„,  the  Tait  equation  gives 
decidedly  better  results  than  the  generalized  one-parameter  equation 
.(4.8). 
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Pressure  calculations  are  shown  in  Table  A-6.   The  errors  of 
calculation  of  equation  (4.8)  and  the  Tait  equation  are  of  the  same  order 
of  magnitude,  but  generally  the  Tait  equation  gives  better  agreement 
with  the  experimental  pressure. 

Where  extensive  volumetric  data  for  a  specific  liquid  are 
available,  it  is  preferable  to  extract  the  Tait  constants  and  use 
equation  (4.11)  for  compression  calculations.   Equation  (4.8),  however, 
has  the  appeal  both  of  generality  and  a  firm  basis  in  molecular  theory. 
Its  dependence  on  a  single  characteristic  parameter  reinforces  the 
simplicity  of  use.   Where  there  is  limited  experimental  data,  the  use 
of  equation  (4.8)  will  give  results  of  adequate  accuracy  for  compression 
calculations.   Further,  it  is  equally  possible,  although  it  remains 
uninvestigated,  that  a  characteristic  parameter  extracted  directly 
from  volumetric  data,  i.e.,  equation  (4.8)  rather  than  from 
compressibilities,  will  yield  results  of  accuracy  equivalent  to  those 
of  a  particular  equation  of  state  such  as  the  Tait  or  Huddleston 
equations. 

4.3  Temperature  and  Pressure  Dependence  of  DCF  Integrals 

The  corresponding  states  correlations  for  C  „  and  C  „  describe 
their  variation  with  the  reduced  solvent  density.   The  solvent  density, 
however,  is  a  unique  function  of  the  system  temperature  and  pressure. 
Knowledge  of  the  temperature  and  pressure  behavior  of  the  DCF  integrals 
directly  yields  the  state  dependence  of  the  related  thermodynamic 
properties.   We  have,  in  the  1-2  binary. 
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Taking  a  derivative  w.r.t  temperature  at  constant  pressure 


3C22 

3T 

P 

3T 

I             J 

(4.13) 


3F, 


3p: 


3p 


o  -1 


3P 


P  ' 


3P_ 
3T 


(4.14) 


Now 


dp, 
W 


P2<2  =if/(F^  -  1) 


(4.15) 


1      "     J 

P 

f  ''^1  f 

i  ^4  \ 

hi^\ 


- 1] 


3P_ 
3T 


(4.16) 


Now,  writing  the  isothermal  equation  of  state,  equation  (4. 8), with 

s 
P„,  the  saturation  pressure  of  the  solvent  as  the  base  pressure, 

Po 


P  -  P^ 
* 

RTp^ 


'2 

P, 


[F^(52)  -  l]d52 


[F^(32)  -  l]dp  - 


-s 
^2 


'~o 


[F^(P2)  -  l]d3 


(4.17) 


(4.18) 


=  F^(P2)  -  F^(P2) 


(4.19) 


then 


RTp, 


3P_ 
3T 


dP 
dT 


P  -  P'" 

2  * 
RT  P2 


3T 


,-s. 


P2  2 


(4.20) 
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RTp, 


3P 


dF 

dT 


P   -   P" 

2   * 
RT  P2 


dT 


(4.21) 


3P_ 
3T 


=   RTp, 


P   -   P" 

2   * 

RT  P2 


dT 


dP- 

dT 


(A. 22) 


Substituting    (4.34)    in    (4.29) 


f  '<^  1 

[dFj 

8T 

P 

[dP°J 

^/[F^(P2)    -   1] 


RTp, 


P  -  P^ 

2   * 
RT  p2 


dT 


dT 
dT 


(4.23) 
Then,  to  determine  the  temperature  derivative  of  C„„,  we  need 
(1)    the  compressibility  correlation,  equation  (4.3) 
(li)    an  equation  (relation)  for  the  saturation  pressure  for 
the  solvent  as  a  function  of  temperature;  we  shall 
use  one  of  the  generalized  equations  suggested  by 
Reid  and  Sherwood,  e.g.,  Antoine  equation,  or  Miller- 
Thodos  equation. 

(iii)   an  equation  for  the  density  of  saturated  liquid;  here 

32 
we  use  the  generalized  Rackett  equation, 

(iv)    the  generalized  isothermal  equation  of  state,  equation 

(4.8). 

Differentiating  (4.12)  w.r.t.  pressure 


3C 


22 


3P 


dF^(p2) 
dp° 


9p, 

w 


(4.24) 


(4.25) 
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Similar  expressions  can  be  obtained  for  the  state  dependence  of  C 


12 


-  C 


12 


*   N 

■k 


.62 


f^Cp^) 


(4.26) 


and 


ac^2 

3T 
I               J 

P 

['"^l 

.'"2] 

(          *     "1 
^1 
* 

I   V2  J 

.62 


^/[F^(P2)  -  1] 


9P_ 

3T 


8P 

•                              • 

T 

f   dF2| 

r     *  -( 
^1 
* 
I   V2  J 

,62 


^/[F^(P2)  -  1] 


(A. 27) 


(4.28) 


We  can  now  write  the  pressure  and  temperature  variation  of  the  isothermal 
compressibility  and  partial  molar  volumes  using  equations  (4.23)  through 
(4.28). 


3T 


(1  -  C^2^ 


3p. 
3T 


(4.29) 


3k, 
3F 


(1  -  ^2) 


3C 


22 


3P 


3p, 


(4.30) 


5 

3P 


-o 
^1 


3C 


22 


3T 


(1  -  C22)  I 


-o 

^1 


P   (1  -  C22) 


3P 


T   p°(l  -  C°2) 


3C 


22 


3T 


P   P2(l  -  C22) 


3C 


12 


3P 


-o 
^1 


T   p, 


3C 


12 


3T 


-o 

^1 


P   P, 


The  quantity 


3p, 


dF^(pp 

is  obtained  from  — 7;;; —  in  equation  (4.23). 

ai 


The  details  of  this  derivation  are  in  Appendix  D. 
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Calculations  for  general  gas-solvent  systems  are  shown  in 
Table  4-7.   These  are  order-cf-inagnitude  estimates  of  the  pressure  and 
temperature  derivatives  of  the  compressibility  and  partial  molar 
volume.   Ultimate  accuracy  of  prediction  depends  on  close  correlation 
of  C°2  and  c"^. 


72 


I 

M 

< 


CO 
M 
M 

§ 

u 


o 
o 


M 

PC 
H 

O 

w 

> 

l-t 
H 
< 

> 
M 

w 

Q 

•s 

:=> 

CO 
CO 


i 

PL4 


Ui 

PU 

0 

<u 

vO 

00 

o\ 

in 

iH 

v£> 

iH 

0    r- 

i-H 

<r 

CN 

n 

CN 

-* 

.H 

CM 

i-H 

cs 

1  > 

H 

0 

• 

• 

• 

• 

• 

• 

• 

• 

• 

fO 

<D 

a 

^— ' 

v^ 

H 

'  0 
0 

B 

-■—^ 

m 

/->. 

00 

in 

^'**^ 

CO 

cr» 

p~ 

iH 

a^ 

tH 

(N 

iH 

cs 

01 

0 

0 

0 

t-l 

0 

0 

0 

0 

0 

0    .- 

iH 

• 

• 

• 

• 

• 

• 

m 

• 

• 

1   > 

Oi 

0 

1 

1 

1 

1 

1 

1 

1 

1 

1 

<T> 

fO 

fc 

— ' 

s.^ 

'  0 

u 

0) 

1-1 

CO 

^s 

y'-N 

1-1 

/-^ 

<r 

^— \ 

^-N 

•-^ 

0 

iH 

v£) 

\o 

r^ 

r- 

\0 

rH 

r-l 

vO 

04 

0 

m 

0 

0    r-H 

e 

• 

• 

• 

• 

• 

• 

• 

• 

• 

* 

• 

1   > 

CO 

iH 

f-4 

o\ 

m 

<N 

in 

in 

CN 

vr 

00 

r^ 

cn 

0 

vO 

m 

%o 

ir> 

00 

CO 

n 

^ 

in 

n 

o^ 

~3- 

m 

u 

N-/ 

>^ 

^-^ 

' — ' 

^-^ 

N_X 

VO 


u 


&4 


c 
cr 

01 


CM 

0 
CO 

00 
0 

1 

\0 

1 

00 
1 

1 

00 

1 

m 
00 


1 


CT\ 

-* 

^ 

CO 

00 

00 

00 

CM 

CM 

CM 

0 

0 

0 

0 

0 

0 

• 

0 

0 

0 

• 

0 

• 

0 

• 

>3- 
0 

i-t 

H 

4-1 

CM 
1 
^-\ 

6 

4-1 

(d 

-.006 
(-.007) 

-.006 

(-.011) 

-.006 

ON 

• 

1 

008) 
009 

a. 

C 

3 

1    1 

O        S 


CO 

CM  -* 
O  <N  ^-' 

U 

cr 
« 


1^  .H 


.H  Cl   1-1 

m  r^  in 


r-.  00  r~- 
0\  C%  0^ 


\o  sr  vo 


000 
000 


I     I     I 


00  <n  CO 
000 


CM    tH 


CM   CM    CM 


CO 

CO 

0 

CM 

1 

0 
CM 

1 

1 

vO 

CO 

1 

vO 
vO 

vO 

m 

en 
1 

m 

1 

CU     4J 
(0 


M  o    I 


O 


CO 

-3- 


O 


en 


o 


.-I 


CO 

ON 
CM 


00 
CM 


00 

00 

en 

a\ 

CM 

CM 

n 

01 
CD 

« 

4J 

(3 
01 
U 
Qt 

a. 

c 

•H 

CO 
(U 

3 
iH 
CD 

> 


s 

0)   CM 

4J 
(0 

CO 


0 


CM 

X 


CO 

•H 
X 

CO 
U 

c 


00 

rH 
X 
00 

o 

G 


CM 


\0 

X 

vO 
U 


CM 

X 


vO 

X 


CM 

z 


u 


CM 

X 


u 


Cvl 

z 


0) 

6 

-H 
0) 

cx 

X 

w 


73 


REFERENCES  FOR  CHAPTER  4 

1.  J.  S.  Rowlinson,  Liquids  and  Liquid  Mixtures,  Plenum  Press,  New  York, 
1969. 

2.  P.  S.  Snyder  and  J.  Winnick,  "The  Pressure,  Volume  and  Temperature 
Properties  of  Liquid  n-Alkanes  at  Elevated  Pressures,"  5th  Symposium 
on  Thermophysical  Properties,  ASME  (1970). 

3.  A.  K.  Doolittle  and  B.  Doolittle,  AIChE.J.,  6_,  157  (1960). 

4.  W.  G.  Cutler  et  al.  .  J.  Chem.  Phys.,  29..  727  (1959). 

5.  R.  C.  Reid  and  T.  K.  Sherwood,  The  Properties  of  Gases  and  Liquids, 
McGraw-Hill,  New  York  (1966). 

6.  B.  H.  Sage  et  al. ,  I. E.G.,  26,  1218  (1934). 

7.  W.  B.  Kay,  I. E.G.,  32.,  358  (1940). 

8.  R.  E.  Gibson  and  0.  H.  Loeffler,  J.  Am.  Chem.  Soc.  ,  61,  2877  (1939). 

9.  P.  W.  Bridgman,  Proc.  Am.  Acad.  Arts  &  Sci,  ,  69_,  389  (1934). 

10.  L.  G.  Schomack  and  C.  A.  Eckert,  J.  Phys.  Chem.,  7J*_,    3014  (1970). 

11.  F.  I.  Mopsik,  J.  Chem.  Phys.,  5£.  2559  (1969). 

12.  G.  A.  Holder  and  E.  Whalley,  Trans.  Far.  Soc,  58.,  2095  (1962). 

13.  D.  W.  Newitt  and  K.  E.  Weale,  J.  Chem.  Soc,  London,  3092  (1951). 

14.  D.  A.  Lowitz  et  al.  ,  J.  Chem.  Phys.,  74^,  3014  (1970). 

15.  J.  S.  Rowlinson,  Liquids  and  Liquid  Mixtures,  Butterworths, 
London  (1958). 

16.  K.  Kumagai  and  T.  Toriumi,  J.  Chem.  Engr.  Data,  16^,  293  (1971). 

17.  R.  H.  Schum  and  0.  L.  I.  Brown,  J.  Am.  Chem.  Soc,  75^,  2520 
(1953). 

18.  J.  F.  Connolly  and  G.  A.  Kandalic,  C.  E.  P.  Symposium  Series,  59, 
8  (1963). 

19.  Y.  Ng  and  J.  Walkley,  J.  Phys.  Chem.,  72,  2274  (1969). 

20.  J.  E.  Jolly  and  J.  H.  Hildebrand,  J.  Am.  Chem.  Soc,  SjO,  1050 
(1958). 

21.  J.  C.  Gjalbaek  and  J.  H.  Hildebrand,  J.  Am.  Chem.  Soc,  21,  1072 
(1950). 


74 


22.  J.  Horiuti,  Scientific  Papers  of  Institute  of  Physical  and  Chemical 
Research,  17,  125  (1931). 

23.  J.  P.  Kohn,  "Volumetric  and  Phase  Equilibria  of  Methane-Hydrocarbon 
Binary  Systems,"  A.I.Ch.E.  Meeting,  Houston  (1967). 

24.  I.  Kritchevsky  and  A.  Illinskaya,  Acta  Physicochemica  U.R.S.S., 
20,  327  (1945). 

25.  S.  D.  Chang  and  B.  C.  Lu,  Can.  J.  Chem.  Engr. ,  48,  261  (1970). 

26.  V.  Berry  and  B.  H.  Sage,  J.  Chem.  Engr.  Data,  4^,  204  (1959). 

27.  M.  Orentlicher  and  J.  M.  Prausnitz,  Can.  J.  of  Chem.,  45_,  595  (1967) 

28.  J.  M.  Prausnitz,  Molecular  Thermodynamics  of  Fluid  Phase  Equilibria, 
Prentice-Hall,  Englewood  Cliffs  (1969). 

29.  E.  C.  Lyckman  et  al. ,  Chem.  Engr.  Sci.,  _20,  685,  703  (1965). 

30.  A.  V.  Itterbeek  et  al. ,  Physica,  29,  742  (1963). 

31.  E.  H.  Amagat,  Ann.  de  Chim  et  Phys.,  29,  68-136,  505-574  (1893). 

32.  C.  F.  Spencer  and  R.  P.  Danner,  J.  Chem.  Engr.  Data,  Y]_,    236 
(1972). 


CHAPTER  5 
SOLUTION  THEORY  FOR  SUBCRITICAL  SYSTEMS 

In  previous  sections  we  discussed  the  pressure  and  temperature 
dependence  of  thermodynamic  properties,  and  of  the  corresponding  DCF 
integrals  for  systems  of  constant  composition,  either  pure  solvent 
or  a  solute  at  infinite  dilution.   Here,  we  consider  the  prediction 
of  theinnodynamic  properties  in  subcritical  systems  of  variable 
composition.   The  temperature  or  pressure  of  the  systems  may  be  fixed 
or  variable  in  accordance  with  the  Gibbs  phase  rule. 

The  corresponding  states  correlations  for  the  DCF  integrals 
derived  earlier  for  limiting  compositions  can  be  extended  to  mixtures 
by  the  concept  of  a  pseudo  pure  fluid,  analogous  to  a  one-fluid 
macroscopic  theory.   Here  the  mixture  properties  are  represented  by 
those  of  a  pure  fluid  whose  critical  parameters  depend  on  the 

composition  of  the  mixture.   This  rule  is  similar  in  concept  to  Kay's 

1  2 

rule  of  pseudo  criticals,  and  Prausnitz  and  Chueh's  mixing  rule. 

The  mixing  rules  for  the  critical  parameters  are  necessarily  empirical, 
unless  they  are  based  on  a  rigorous  molecular  theory.   This  approach 
will  be  used  to  study  the  compressibilities  of  binary  and  ternary 
liquid  mixtures. 

Finally,  when  faced  with  a  lack  of  experimental  data  in  the  mid- 
composition  range  of  a  solution,  one  can  postulate  (on  the  basis  of 
empirical  knowledge)  rules  for  the  composition  dependence  of  the  DCF 
integrals.   Two  of  the  simple  rules  are  the  linear  arithmetic  mean  and 
the  linear  harmonic  mean.   Expressions  for  the  thermodynamic  properties 
with  the  linear  rule  have  been  presented  earlier. 
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5.1  One-Fluid  Theory 

Leland  and  Chappelear  discuss  several  methods  of  applying 
corresponding  states  theory  to  mixtures.   The  microscopic  principle 
rests  on  the  specification  of  the  parameters  of  a  pure  fluid  pair 
potential  function  which  generates  the  properties  of  the  mixture  at 
the  same  temperature,  pressure  and  density.   On  the  macroscopic  level, 
the  problem  becomes  one  of  defining  a  relation  between  the  critical 
properties  of  the  pure  components  and  the  pseudocritical  properties 
of  the  hypothetical  pure  fluid.   One  of  the  simplest  examples  of  these 
is  Kay's  rule,  and  an  equally  simple  rule  is  used  here. 

It  is  postulated  that  the  properties  of  the  fluid  mixture 
correspond  to  those  of  a  pseudo  pure  fluid,  whose  characteristic 
parameters  depend  on  the  composition  of  the  mixture.   The  pseudo  pure 
fluid  follows  macroscopic  states  and  its  thermodynamic  properties  are 
described  by  the  correlations  derived  earlier. 

The  macroscopic  correlation  for  compressibility,  equation  (4.3), 
contains  only  one  characteristic  parameter  for  each  substance,  the 
characteristic  volume.   For  the  mixture,  this  is  defined  to  be  the 
simple  mole  fraction  average  of  the  pure  component  volumes 

*         * 
v  ^   =  Z  x.v.  (5.1) 

mix      1  1 

Then,  for  the  mixture 

1 


1  + 


PkRT 


.  Hi^,^)  (5.2) 

mixture 


where  F(p  .  )  is  the  same  generalized  function  of  p  as  equation  (A. 3). 
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Equation  (5.2),  with  the  rule  of  (5.1),  has  been  used  to 
determine  the  isothermal  compressibilities  of  the  binary  mixtures  of 
benzene-cyclohexane,  methyl  acetate-water,  nitrobenzene-aniline,  and 
of  the  ternary  mixture  n-hexane,  n-heptane,  n-octane.   The  results  are 
shown  in  Tables  5-1  and  5-2.   In  all  systems,  there  is  good  agreement 
between  the  calculated  and  experimental  compressibilities  over  the 
entire  composition  range.   This  clearly  demonstrates  the  efficacy  of 
the  approximations  and  the  composition  mixing  rules  of  (5.2)  and  (5.1). 
In  the  binaries  of  the  hydrocarbons,  and  of  nitrobenzene-aniline,  there 
is  a  small,  but  consistent  deviation  between  the  experimental  and 
calculated  compressibilities.   However,  this  probably  results  from  an 
incorrect  prediction  of  the  pure  compressibilities  by  equation  (4.3) 
and  should  not  be  immediately  attributed  to  a  failure  of  the  mixing 
rule.   A  possible  remedy  is  to  adjust  the  characteristic  parameters  of 
the  pure  components  such  that  equation  (A. 3)  predicts  their  experimental 
compressibilities  correctly.   The  results  of  such  calculations  are 
shown  in  Table  5-3,  and  reveal  marked  improvements  in  the  benzene- 
cyclohexane  and  aniline-nitrobenzene  systems.   This  stratagefiv  has  a 
mixed  effect  in  the  water-methyl  acetate  system,  which  is  a  possible 
reflection  of  the  complex  interactions  in  this  mixture. 

The  volume  of  the  binary  mixture  is  required  in  all  of  the 
above  calculations  and  has  been  assumed  an  experimentally  accessible 
quantity.   This  may  not  always  be  the  case  and  the  simplest  approximation 
under  these  conditions  is  one  of  ideal  mixing.   The  effects  of  this 
assumption  on  the  predictions  of  the  one-fluid  theory  are  shown  in 
Figs.  5-4  and  5-5. 


78 


Q 

W 

H 

< 

.J 

:=> 

CJ 

I-) 

<; 

u 

en 

w 

pa 

;=> 

H 

X  i^ 

h-t    o 

:s  M 

H 

3    M 

c-  ;< 

M    O 

^    K^ 

l-M 

>->   A< 

Pi    < 

<--i 

< 

1 

K   a 

lO 

M    t-l 

PQ  ;r> 

w 

J 

>-l 

(^  f^ 

^ 

o 

H 

W    Pi 

w  ;d 

H    Ph 

H 

M    O 

>J   Q 

n  ;j 

«   W 

M    (/) 

M    p^ 

00 

■    W      >H 

c;  m 

CM 

S 

o 

o 

u 
c 

0) 

«-( 

OS 


CO 
Q 


P. 


00 


CO 


NO 


cri  CM  vo  00  CM  r^ 

w  CO  r~.  vo  vD  in 


CM  i-t  o  (yi 

CM  tH  O  00 


CM  ON  m  r-. 

CM  CO  00  m 


o  <yi  o>  00  r-»  r^ 
lA  ^  <i-  ^  -*  <T 


f-\    r-i    ,-i    t-i    r-i    r-i 


CM  CM  (M  t-( 


6 

-* 

■M 

o 

nj 

O 

ro  LT)  CTi  CM  r^  CM 

cr»  LO  n  n 

<T>    OX    1-1    CM 

iH 

^^^ 

r-( 

(Tv  00  r^  r--  \o  vo 

est    r-i    a    (^i 

CM    <T>    00   \0 

00  00  r^  \D  lO  vr 
-*  ^*  «*  sr  p^  -J- 


>  B 


. 

(U 

r-l 

t-) 

>< 

JJ 

O 

•H 

n 

6 

y 

a. 

> 

x: 

o 

Q) 

o 

CT»  »— t  0\  f-^  vO  t— I 

Lo  CO  a\  CM  m  m 

i^  r^  r^  00  00  00 


CM  CM  CM  f-l 


o  rH  i-H  a\ 
r^  O  m  m 
\o  p~  r^  r^ 


O  O  O  -H 
r-  -<r  (7>  tn 

v£>  r^  r~  00 


P--  -*  iH  CM  <■  00 

CN  i-t  •-(  i-H  .-I  i-H 


Csl  CM  CM  CM  CM  C>l 


O  r-l 


n  ~3-  in 
^  00  o^  O 


CM  CM  CM  CM 


CO  .H 
O  00  CO  O 


CM  CM  CM  CM 


00  00 

O  lO  f---  o 


CO  CO  CO  n  en  CO 


o  00  r^  v^  m  o 
o  -<r  CO  CM  a\  o 


iri  CTi  1— I  \o  r^  CTv 
CM  CO  lo  vo  r^  00 
<r  <r  -a-  <»•  <f  -<f 


o  ^  in  vo  CO  <j- 
CO  CO  in  \o  o  cvi 
o  00  CM  vj  m  in 


o^  \o  in  m 

vO  00  O  CM 
CO  CO  -*  ~3- 


rH  CM  0\  O 

iH  \0  O  CO 
CM  CM  v£)  O 


ON  ON  CM  ON 

vo  O  <r  00 

CO  <•  -3-  sj- 


r^  CM  -* 

r^  CTi  CM 
<r  ^  m 


I-H  tn  vo  CM  00  in 

Csl  O  CTl  ON  CO  00 
CO  CO  CM  CM  CM  CM 


I-H  t— I  iH  O  r-~  v£) 

in  o  m  .H  00  iH 

vO  iH  CN  CTv  r~~  -a- 


<r  r^  rH  m  00  iH 
in  in  \o  vo  vo  r^ 


00  CO  00  <• 
CO  ^3-  ^T  m 


00  ON  00  rH 

CO  -d-  in  r~- 


CM  CO  in  CO  CM  i-H 

O  0\  CTv  ON  ON  ON 


O  CO  <!•  r^ 

CJN  r-^ 

00  m 

t-H  -^  00   .H 

cvi  iH  m  CM 

<-i  in 

CO     r-l 

CO    00    ON    ON 

CN  -vT  vo   00   O 

O  CO  vO   O 
•       •       •       • 

CO  vO  o 

•        •        • 

^  \0  r~~  00  O 

CO 

CO 

CO 


CO 


CO 

CO 
CO 


00 
ON 
CM 


vf> 


CM 


e 

(U 

■u 


00 
rH 

00 

o 
c 


o 

(3 


l-» 
X 

c 


NO 

o 


CO 

r-t 
X 
CO 

o 

c 


CM 

o 
in 

u 


CM 

z 

m 

X 

M3 
U 


79 


d) 

u 

c 

9) 

U 

0) 

<4H 

-* 

OJ 

m 

OS 

cn 

« 

4J 

d 

o 

-3- 


V  I 


u 

a 

w 


\o  CS  00 
O  O  <T> 


o\  in  ON 

<M  CM  t-l 


O  O  O  r-~ 

<r  -*  >3-  00 


a 

r 

u 

• 

o 

<o 

o 

fn  cy. 

<r  ON  o 

r--  O  -*  CTi  vo 

\0    lA   IT)  O 

rH 

>-• 

i-H 

•-i    O 

o  (TV  a» 

vr  n  CM  1-1  tH 

CO    CO    n  CTv 

tH   .H 

r-l 

iH   .-(    .-1    iH  rH 

01 

c 

c 
o 
u 


00  o\  o  I-H  in 
in  in  in  m  vo 

00  CO  CO  00  00 


in  in  r~-  o  -* 
r^  r~-  o  (~~  vo 

r^  r^  r~  r--  r~~ 


CM  CM  Ol  CM  CM 


CM  CM  CN  CM  CM 


00  vo  iH  m 

CM  r^  O  ro 

r^  I — ^  00  00 

CM  CM  CM  CM 


I 

m 

H 


■H 

*  e 
> 

•  0) 

I-H  1-1 

KU  O 

>    X    O 
0)  u 


o  o  o  o  o 


o  o  o  o  o 


i-i  r-»  r^  ON  in 
1-1  C7\  00  v£>  in 
fO  CM  CM  CM  CM 


O  <T»  O  vO  O 
00  00  CO  CO  o 


1— I  r^  n  CTv  m 
1— I  c^  00  ^  in 

m  (N  CM  CM  CM 


m  ro  vo  m  m 

O  O  CM  iH  CM 


00  ro  <Ti  -3-  0> 

o  CD  a\  o\  oa 


CM  P^  rvl  r^  CM 
i-l  O  O  o>  <n 


<■  o  o\ 

in  o\  i-( 

*  •  •  • 

-a-  CM  00  in 

\0  r^  r^  CM 

CM 


vO  ON  r» 

vo  m  o  in 
vo  CM  cTi  en 

•  •  •  • 

CO  >sO  r^  ON 
CM  CM  CM  r^ 


in  O  in 

CM  in  r^ 


in  o  m 
O  CM  m  r-~  o 


vo  -<r  o 

M  00  00 

O  <r  r-. 

O  iH  iH  iH  O 


CO 

o\ 

CM 


CO 
CM 

CO 


Ok 
CM 


0 
u 

0) 


CM 

t-i 

X 
vO 

o 


as 
u 


CM 

1-1 

U 


\0 


M 
4) 

« 

I 


0) 


iJ    0) 


80 


+ 


X 


I>1 

--3- 


o 

\0 


CO 


00 


00 


a 

cc 
o 

en 

in 

o 

in 
o 

O 

en 

o 
in 

00 

I-l 

r-l 

iH 

iH 

iH 

>-i 

r-l 

.-1 

CO 

w 
Pi 

H 


M 


I 

in 

w 


H 

O 

J" 
H 


<N 

O 

CO 

m 

O 

00 

.-( 

<r 

X 

cr\ 

O 

<r 

vO 

iH 

iH 

\£> 

CO 

•H 

a\ 

cr> 

00 

r^ 

O 

Oi 

00 

r~. 

E 

, 

• 

• 

• 

• 

• 

• 

• 

(Q. 

CM 

£M 

CM 

CM 

PO 

C>l 

CM 

CM 

(Tv 

CT\ 

CTv 

CT> 

CO 

ro 

CO 

CO 

m 

ro 

CO 

n 

r^ 

r^ 

r«. 

r^ 

vO       \0       ^0       vO 
CM       CN       CM       CM 

vr     -<r     -*     <r 


o\  CT»  CT>  o^ 
CO  CO  CO  CO 
-*     -*     ^    '^ 


Q) 

<■ 

in 

00 

r» 

vO 

CO 

r-« 

• 

r-t 

CO 

CO 

cr> 

CM 

vO 

r» 

o 

in 

X   *-» 

g 

in 

o 

a\ 

CM 

o 

r- 

r-» 

CT> 

■H    a. 

^ 

« 

• 

• 

• 

• 

• 

• 

« 

E    X 

CM 

r^ 

<j^ 

St 

vO 

o 

CO 

r~ 

0) 

a 

■<r 

>a- 

vr 

in 

-* 

m 

m 

in 

o 

<-» 

r-t 

1-H 

rH 

iH 

t-t 

tH 

1-4 

H  i^ 


ro 

00 

CO 

CO 

CO 

00 

CO 

CO 

r-~ 

CTl 

r-l 

CO 

r^ 

a\ 

rH 

CO 

CM 

CM 

CO 

CO 

CM 

CM 

CO 

CO 

6 

0) 
4J 
U) 
>^ 
CO 


\0       00 


vo     CO 


K 

K 

X 

X 

33 

X 

vO 

p- 

00 

vO 

r-~ 

00 

u 

o 

o 

U 

U 

CJ 

^^ 

^« 

s>e 

S^! 

^e 

s< 

vO 

CM 

CM 

.-• 

vO 

CO 

CM 

r^ 

O 

.-1 

iH 

p~ 

• 

» 

• 

• 

• 

• 

in 

r-4 

CO 

a\ 

i-i 

a\ 

CO 

CO 

CO 

iH 

sr 

CO 

00 

at 
a 

0 

0) 

M-l 


81 


Pi 

Q. 


a 
Pi 

H 


(0 

o 


X  -H 
(U     £ 


o 

o 

CM 

in 

o 
rsi 

o 

CO 

CM 

en 

-3- 
O 

00 
o 

CM 

04 

CM 

io 

CO 

i-l 

CM 
CO 

CO 

CM 

o 
en 

en 
en 

CN 

ON 

00 

1-H 

— 1 

CO 

in 

o 

CO 

CO 

in 

CM 

cn 

CO 

00 

in 

00 

CO 

•■o 

rn 

m 
n 

en 

-3^ 
CM 

tM 

\o  in 
r-  o 

CO 

m 
o 

-4" 

in 

o 

\0 

cn 

o 

in 

00 

cn 

VD 

vO 

00  CO 
m  m 

in 

00 

m 

VO 

cn 

in 
rn 

VO 

cn 

CM 

cn 

VO 

cn 

CO 
en 

vO 
C^4 

cn 

CO 

cn 

a\   00 

CM  00 

00 

00 

CO 

CM 

m 

CO 

in 

CSl 
CO 

00 
CM 
CO 

CJ\ 

CO 

CM 
CO 

CO 

i-H 

cn 

o 

00 

VO 

o 

CO 

CS  <M 

<M 

CM 

CM 

CM 

CM 

CM 

CJ 

CM 

CM 

CM 

CM 

CM 

CM 

CO 
r-t 

CO 

CO 

O 

o 
in 
in 

00 
rH 

vO 

in 

vD 

CO 
00 

in 
in 

c-l 

v£) 
vO 

CO 

o 
m 

O 
vO 

VO 

LO 

CM 

ON 

o 

^  CO 

o  ov 

ON 

cn 

O 

CO 

a. 

in 

in 

VO 

00 

vO 

CO 
CM 

vO 
CM 

CM 

C3 

»-] 

b 

b 

O 

cn 

M 

1 

U 

m 

M 

H 

w 

M 

►4 

»J 

^ 

P3 

H 

M 

M 

CM   -*    CSl    0\ 

CO  in  iH  m 

VO  -^   CO   <j> 


rH  .-H  CO 
i-<  CO  vD 
CM    cn    CO 


sj-  in  o  vo 
CO  <r  CO  CM 
Ov  CM  vO  r^ 


CO  -cr  M 

iH  O   >H 

VO  o  cn 


• 

CO 

• 

• 

o 

• 

vO 

00 

CO 

• 

o 

CO 

00 

• 

-3- 

• 

CO 

• 

00 

o 

Ov 

o\ 

00 

cjv 

r^ 

VO 

^ 

in 

VO 

r^ 

vO 

r~ 

r-- 

CO 

Cvl 

CM 

CM 

CM 

CM 

CM 

<r 

<r 

cn 

<r 

CO 

00 

vO 

<r 

r^ 

CO 

vO 

•<r 

O 

rH 

-^ 

00 

tH 

C7V 

CTv 

'-n 

c 

CO 

CO 

r-^ 

.-H 

CO 

00 

CO 

00 

a\ 

Ov 

CO 

-J- 

CO 

CM 

■-I 

m 

eg 

o 

<3- 

r~. 

<T 

VO 

• 

• 

00 

• 

CM 

• 

m 

• 

r^ 

CM 

• 

-3- 

• 

vO 

CO 

• 

.-1 

• 

rH 

• 

.H 

O 

*     CM    6 

>    ^ 

U 
O 


•K     .H 


S 
(U 

to 


CM 


cr> 
t-i 

CO 


o 
o 

cn 

00 
CM 


CO 

cn 


(U 

c 
<u 

C 

<u 

O 


c 


c 
< 


»>. 

o 
cn 


VO 

m 


00 
CI 


O 

vO 
CM 

-a- 


in 

VO 

O 
a\ 


CO 
CO 
CO 


00 

cn 
«o 


o 

c 

<d 

OJ 

X 

a 

01 

ca 

^ 

u 

o 

(X 

1-1 

0) 

u 

K 

^ 

C 

0) 

o 

c 

a 

ed 

01 

4-> 

N 

U 

c 

o 

(U 

BQ 

c 

(D 

W 

0) 

Oi 

<n 

M 

U 

00 

01 

<U 

vO 

4-) 

g 

• 

dJ 

w 

in 

B 

U 

CM 

q) 

ca 

CM 

o. 

o. 

4-1 

<n 

4J 

0) 

a\ 

a 

c 

e>» 

0) 

o 

c 

a. 

o 

E 

o. 

o 

B 

u 

o 

o 

0) 

0) 

3 

u 

o. 

3 

u 

o. 

T3 

4) 

V 

4J 

"O 

4-1 

Jd 

0) 

CO 

» 

4-1 

3 

CO 

1-1 

3 

T3 

1-1 

<a 

0) 

•n 

c 

i-» 

4J 

< 

:=) 

>s 

HI 

J3 

4J 

*J 

(U 

x^ 

y'-N 

0) 

u 

■H 

CN 

S  < 

>-• 

^^ 

82 


5.2  Rules  for  Composition  Dependence  of  DCr  Integrals 

The  thermodynamic  properties  of  a  binary  system  are  completely 
defined  by  the  three  integrals  C^^,  C^  _  ^""^  ^jl   ^^'^^S^  ^^^  expressions 

12       2 

=  1  -  x,C,,  -  x„C„„  -  2x.x„C-_  (5.3) 


p<RT        1  11    2  22     1  2  12 

^;^  -  (1  -  x^C^^  -  X2C^2>  <^-^) 

and 


°  X  2       2 

T,P    ^         1  -  x^C^^  -  X2C22  -  2x^X2C^2 

(5.5) 


The  quantities  C^ - ,  C^„  and  C-^  vary  with  the  composition  of  the 
mixture.   At  any  composition,  they  can  be  determined  only  through  a 
simultaneous  knowledge  of  all  three  thermodynamic  properties. 

However,  experimental  data  are  often  available  only  at  the 
composition  extremes,  and  the  variation  of  the  C   over  the  composition 
range  must  be  hypothesized  to  calculate  mixture  properties  at  other 
compositions.   Equations  describing  this  composition  dependence  may  be 
proposed  and  the  calculated  thermodynamic  properties  compared  with 
experimental  data.   Good  agreement  for  a  wide  range  of  systems  is 
evidence  for  the  accuracy  of  the  composition  rules  which  can  then  be 
(confidently)  applied  to  other  systems.   The  two  composition  rules 
studied  are  the  arithmetic  mean  (5.6)  and  the  reciprocal  mean  (5.7) 
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The  superscripts  *  refer  to  quantities  at  x^  =  0.  while  °   refers  to 
quantities  at  x^  =  1.   The  above  relations  are  then  substitution 
into  equations  (5.3).  (5.4)  and  (5.5)  to  yield  expressions  for  thenno- 
dynamic  properties  of  the  mixture  in  terms  of  DCF  integrals  at  the 
limiting  compositions,  which  in  turn  are  calculated  as  follows 


(5.8) 


1 

0    0 

v°        1  -   C° 
2        ^        ^12 

O           ,               0 

■  ailnY2 

Sx^ 

o 

-   (1 

(5.9) 


c°)   ^'-^12)' 

T,P  (1  -  Cj^i) 


where  the  symbols  have  their  usual  meaning.   Similar  expressions 

1     *    *       * 
relate  C^^,  C^^  ^nd  C^^   to  properties  at  the  compositions  x  =  1. 

The  benzene-cyclohexane  and  aniline-nitrobenzene  systems  were 
studied  initially.   Complete  thermodynamic  data  are  available  for  these 

systems  ^'-^'^.S.e.lO  _.        .  ^ 
^  The  composition  dependent  C^^,  C^^  and  C   in 

the  benzene-cyclohexane  system  were  determined  by  simultaneous  solution 
of  equations  (5.3).  (5.4)  and  (5.5)  and  are  shown  in  Fig.  5-1.   Although 
this  binary  system  is  not  very  nonideal.  the  experimental  C.^'s  show  a 
quadratic  dependence  on  composition.  This  is  consistent  with  the 
quadratic  and  higher  order  composition  dependence  of  the  thermodynamic 
properties. 
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Fig.  5-1.   Experimental  DCF  integrals  in  benzene  (1) 
cyclohexane  (2)  system  at  298°K. 
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The  third -order  correlation  function  integrals  C^^,  C^^^* 
C  „„  and  C_„„  of  the  system  can  be  obtained  from  the  complete  composition 
dependence  of  C^,  C,  -  and  C-„.   The  triplet  and  pair  correlation 
function  integrals  are  related  by  the  recurrence  relations  (Appendix  B) 


T,P,N, 


T.P.N, 


*       _*     *      _*  * 


l,j=l,2 


=C°.(l-pv°).C°.,-pv°C°., 


(5.11) 


(5.12) 


i,j=l,2 


These  provide  three  relations  between  the  four  quantities  C..,  at  each 
composition  extreme.   The  fourth  relation  is  obtained  from  the 
density  derivative  of  the  pure  liquid  compressibility  (Appendix  B) 


^iii  "  '^i 


3Pi 


-  C 


ii 


(5.13) 


T,P,N 


J 


The  corresponding  states  correlation,  equation  (4.3),  was  used  to 
determine  the  density  derivative 


ii 


3p. 


The  composition 


derivatives  were  evaluated  analytically  by  fitting  the  experimental 
C.'s  to  quadratics.   The  derived  C...  's  are  shown  in  Table  5-4. 
They  all  have  large  negative  values,  much  more  so  than  the  C   's. 
This  is  in  total  contradiction  to  the  HNC  assumption  for  mixtures  which 
requires  all  C.    =  0  (Appendix  F).   Interestingly  enough,  there  is  a 
consistent  relation  of  the  form 
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which  is  analogous  to  the  relation  often  used  for  the  cross  third 
virial  coefficients  in  gas  mixtures. 

The  essentially  quadratic  composition  dependence  of  the  C 
indicates  that  the  linear  rule  can  only  be  a  first  order  approximation. 
Thermodynamic  properties  calculated  with  equations  (5.6)  and  (5.7) 
are  shown  in  Table  5-5.  Although  the  hypothetical  and  experimental 
C   's  in  the  benzene-cyclohexane  system  differ  by  more  than  10%  at 
some  compositions,  the  resulting  thermodynanic  properties  are  in  much 
closer  agreement  with  experimental  data.   Figs.  5-2  and  5-3  show 
comparisons  of  the  isothermal  compressibility  and  partial  molar  volumes 
with  experiment.   The  calculated  partial  molar  volumes  agree  very 
well  with  the  experiment  over  the  complete  composition  range.   The 
isothermal  compressibility  is  not  as  well  predicted  and  the  calculated 
curve  is  of  a  different  shape  than  the  experimental.  The  maximum 
error  occurs  near  the  mid-composition  range  and  is  about  3%.   Fig.  5-2 
also  shows  the  isothermal  compressibility  calculated  from  the  one- 
fluid  theory  with  nonadjusted  and  adjusted  pure  component  characteristic 
parameters.   Isothermal  compressibilities  calculated  from  the  one- 
fluid  theory  show  the  same  composition  dependence  as  the  experimental 
data. 

For  calculations  of  the  aniline-nitrobenzene  system,  the 
partial  molar  volumes  at  infinite  dilution  were  assumed  equal  to 
the  pure  component  molar  volumes,  i.e.,  ideal  mixing.  The  composition 
dependences  of  the  activity  coefficient,  calculated  with  equations 
(5.6)  and  (5.7),  follow  the  experimental  result  quite  closely  and  are 
shown  in  Fig.  5-4.   The  calculated  isothermal  compressibility  is  of  the 


TABLE  5-5 
RULES  FOR  COMPOSITION  DEPENDENCE  OF  DCF  INTEGRALS.-  CALCULATIONS 


System 

X,      C.^       C,^       C, 


(1)    (2)       °K      1       22        12        11 


Benzene  Cyclohexane  298  0.0     (-37.46)  (-31.45)  (-26.35) 

.2S93   -43.12  -35.50    -29.17 

(-40.64)  (-32.43)  (-29.41) 

.5498   -'.8.21  -39.16    -31.75 

(-47.35)  (-33,63)  (-31.31) 

.7856   -52.82  -42.46    -34.08 
(-76.55) 

1.0     (-57.03)  (-45.47)  (-36.20) 


Aniline  Nitrobenzene  338  0.0  (-60.36)  (-53.70)  (-47.37) 

.25     -63.31    -56.39  -49.87 

.50     -66.27    -59.09  -52.37 

.75     -69.22    -61.78  -54.88 

1.0  (-72.17)  (-64.48)  (-57.38) 


/-' 


Experimental  values  in  parentheses. 
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OF   THERM0DYNA:IIC  properties    IX   SUBCRITICA.L   SYSTEMS" 


p<RT 


39.82 

(37.32) 

39.75 
(36.71) 

38.76 


pv. 


(38.47)         (    .84) 


.871 
(    .872) 

.908 
(    .909) 

.951 


pv. 


(37.20)         (1.0) 


(1.0) 

1.053 
(1.052) 

1.113 
(1.112) 

1.179 


(1.249) 


dx. 


T.P,N, 


(-   .0247) 

-  .0212 
(-    .0175) 

-  .0116 
(-   .0111) 

-  .0072 


(   0.0) 


Si^ny, 


<?x. 


T,P,N, 


(  0.0) 

-  .0086 

-  .0142 

-  .0265 
(-    .0247) 


(61.36) 

( 

.892) 

(1.0) 

60.88 

.916 

1.028 

60.21 

.942 

1.058 

59.36 

.971 

1.089 

(58.38) 

(1 

.0) 

(1.122) 

(-    .394) 

-  .277 

-  .170 

-  .077 
(  0.0) 


(  0.0) 

-  .092 

-  .170 

-  .231 
(-   .271) 
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Fig.  5-2.   Isothermal  compressibility  function  of  benzene 
(1)  -  cyclohexane  (2)  mixtures. 
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Fig.  5-3.   Partial  molar  volumes  in  benzene  (1)  -  cyclohexane 
(2)  system  at  298°K. 
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Fig.  5-4.   Activity  coefficients  in  aniline  (1) 
nitrobenzene  (2)  system  at  338°K. 
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same  form  as  in  the  benzene  cyclohexane  system  and  deviates  from  the 
experimental  value  at  the  mid-composition  range  (Fig.  5-5).   The 
predictions  of  the  one-fluid  theory  with  adjusted  parameters,  on 
the  other  hand,  agrees  well  with  the  experimental  data.   From  the 
limited  study  of  the  two  binary  systems,  one  can  only  speculate  on  the 
efficacy  of  the  linear  composition  rule  for  the  DCF  integrals  in 
calculating  thermodynamic  properties.   However,  the  results  for  these 
systems  indicate  that  with  a  knowledge  of  the  six  limiting  values  of 
the  DCF  integrals,  the  partial  molar  volume  can  be  estimated  to  within 
5%  at  any  composition,  and  the  logarithm  of  the  activity  coefficient 
at  infinite  dilution  to  about  10%.   The  isothermal  compressibility  is 
best  determined  from  the  one-fluid  theory  with  adjusted  parameters, 
using  an  ideal  mixing  rule  for  the  volume  of  the  mixture  if  this  is 
unkno;*'n . 

In  some  situations,  one  may  not  have  all  the  thermodynamic 
information  necessary  to  calculate  the  six  DCF  integrals  at  the  limiting 
compositions.   Under  these  circumstances,  the  only  alternative  to  total 
surrender  is  a  reasonable  and  realistic  approximation  to  relate  the 
C, ,'s  at  each  composition.   The  experimental  C   's  in  the  benzene- 
cyclohexane  system  (Fig.  5-1)  provide  a  clue  in  this  direction.   Over 
the  entire  composition  range,  C^ _  is  approximately  the  arithmetic  mean 
of  C^ ^  and  C„„,  and  this  assumption  is  made  in  subsequent  calculations. 
We  now  propose  to  determine  activity  coefficients  in  the  binary  mixture 
from  the  molar  volumes  and  isothermal  compressibilities  of  the  pure 

components  and  their  partial  molar  volumes  at  infinite  dilution  in  the 

o    o    *       * 
other.   C^^,  C^^,   C„-  and  C  _  are  calculated  as  in  equations  (5.8)  and 
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Fig.  5-5.   Compressibility  function  of  aniline  (1)  - 
nitrobenzene  (2)  mixtures. 
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(5.9).   C^-  and  C--  are  evaluated  from 


*      *     * 

^11  "   ^^12  "  ^2 


C°   =  2C°   -  C° 
^22   ^^12        ^11 


(5.14) 
(5.15) 


The  linear  composition  rules  are  used  to  determine  C   at  intermediate 
compositions.   The  calculated  composition  derivatives  of  the  activity 
coefficients  are  integrated  to  yield  activity  coefficients  at  infinite 
dilution. 


iny. 


x^=0 


£nY°  + 


ax. 


dx^ 
T.P.N, 


=  0  - 


Siny. 


I 


3x 


dx. 


(5.16) 


"""   ''T.P.N, 


The  infinite  dilution  activity  coefficients  in  some  binary 
mixtures  are  shown  in  Table  5-6.   The  agreement  with  experimental 
data  is  satisfactory  only  in  the  acetone-benzene  system.   In  other 
systems,  the  calculated  slope  of  the  activity  coefficient  has  too  large 
a  negative  value.   This  slope  is  calculated  as  a  small  difference  of 
large  numbers,  and  is  accurate  only  when  the  C.'s  at  the  composition 

extremes  are  precisely  known.   Evidently,  equations  (5. 14)  and  (5.15) 

o       *  ' 

do  not  yield  values  of  C  _  and  C   which  satisfy  this  criterion.   The 

calculation  for  aniline-nitrobenzene  at  SSS^K  illustrates  this  point. 
The  deviation  of  C  „  from  the  arithmetic  mean  is  designated  5  „,  i.e., 
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o       * 

The  experimental  values  of  5  _  and  i  „  as  determined  from 

Table  5-5  are  -  .0246  and  -  .0030,  respectively.   These  values  of 
5^«  in  conjunct:ion  with  the  linear  rules  yield  slopes  of  the  activity 

coefficients  which  are  in  close  agreement  with  experiment  (Fig.  5-5). 

*     o 
However,  when  5^„  =  6,^  =  0,  equations  (5.14)  and  (5.15),  the  resulting 

activity  coefficients  differ  by  more  than  20%  from  experiment. 

A  possible  correction  for  this  shortcoming  is  to  include 
information  on  the  limiting  activity  coefficient  of  one  component 
in  the  other,  e.g.,  a  polar  solute  at  infinite  dilution  in  a  non- 
polar  solvent.   This  would  provide  the  precise  knowledge  of  the  C   's, 
and  hence  5  „,  at  the  pure  solvent  composition.   Then,  assuming  a 
constant  value  for  i5^„,  the  C,.'s  at  the  pure  solute  composition  can 
be  calculated.   In  all  likelihood,  the  calculated  activity  coefficient 
for  the  nonpolar  solvent  at  infinite  dilution  in  the  polar  solute 
will  be  in  closer  agreement  with  experiment  than  that  calculated  by 
the  assumptions  of  equations  (5.14)  and  (5.15). 
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CHAPTER  6 

DETERl-IINATION  OF  EXPKRDIENTAL  THERMODYNAMIC  PROPERTIES 
OF  GAS-LIQUID  SYSTEM'S 

The  experimental  method  of  Kay  and  coworkers,  and  Connolly 

2 

and  Kandalic,  vas  used  to  determine  the  P-V-T  relationships  of 

constant  composition  fluid  mixtures.   A  small  precisely  measured 
amount  of  the  constant  composition  fluid  sample  is  contained  over 
mercury  in  an  inverted  capillary  tube,  sealed  at  the  top.   The  tube 
is  mounted  on  a  steel  block,  within  a.   constant  temperature  bath. 
Pressure  is  exerted  on  the  fluid  tlirough  the  mercury.   The  volume  of 
the  sample  at  a  fixed  temperature  and  pressure  is  determined  by 
measuring  the  length  of  the  liquid  column  with  a  cathetometer;  the 
cross-sectional  area  of  the  tube  is  accurately  calibrated. 

6.1  Description  of  Apparatus 

The  apparatus  is  divided  into  pressure  and  vacuum  sections. 

A  symbolic  drawing  of  the  vacuum  section  is  shown  in  Fig.  6-1. 
The  glass  manifold  M,  approximately  100  cm  long,  A  cm  in  diameter,  is 
connected  to  a  mechanical  vacuum  pump  through  the  stopcock  H 
(A  mm  bore),  and  a  cold  trap  with  liquid  nitrogen.   Stopcocks  D,  E, 

F,  G  and  I  are  2  mm  bore.   A  tilting  Mcleod  gauge,  reading  to  1 
micron  Hg,  is  connected  at  D  and  indicates  the  ultimate  vacuum  obtained 
in  the  manifold.   A  Texas  Instruments  pressure  gauge  is  connected  at 

G.  (The  Texas  gauge  has  interchangeable  pressure  measuring  Bourdon 
capsules  for  different  pressure  ranges.)   J,  B  and  C  are  A  mm  bore 
vacuum  stopcocks,  and  K  is  a  ball  and  socket  joint.   The  entire 
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apparatus,  to  the  right  of  K,  can  be  swivelled  through  180°  for 
use  during  calibration.   A  is  a  Cajon  Ultra  Torr  stainless  steel  fitting 
that  provides  a  vacuum  connection  between  the  glass  tubing  and  the 
stainless  steel  compression  block.  (CB)  which  holds  the  experimental 
capillary  tube  ECB.   L  is  a  three-way  glass  stopcock  through  which  mercury 
may  flow  into  the  system.   The  length  of  the  glass  column  from  L  to  A 
is  approximately  80  cms. 

The  pressure  section  of  the  apparatus  is  shown  in  Fig.  6-2. 
The  double  walled  glass  jacket  (VJ)  is  supported  on  the  base  of  the 
compression  block.   The  annular  space  in  the  jacket  is  permanently 
evacuated.   T  and  T„  are  temperature  probes  for  a  Hewlett  Packard 
quartz  thermometer  (HPT).   The  connecting  lines  in  the  pressure 
section  are  1/8"  CD.  stainless  steel  tubing.   Valves  1  through  10 
are  high  pressure  stainless  steel  needle  valves.   SG  is  a  thick 
walled  plexiglass  tube.   It  is  used  as  a  sight  glass  up  to  8  atm 
pressure  to  indicate  the  level  of  mercury  in  the  system,  and  is  shut 
off  during  operation  at  higher  pressures.  MLI  is  a  high  pressure 
mercury  level  indicator,  made  out  of  stainless  steel.   It  is  constructed 
by  connecting  three  pieces  of  insulated  copper  wire  through  pressure 
seals  1-1/2"  apart.   The  wires  are  led  to  lamps  through  a  6- Volt 
battery.   When  the  mercury  in  the  system  rises  to  the  level  of  each 
wire,  the  corresponding  lamp  is  lit.   HG  and  BG  are  Bourdon-type 
pressure  gauges.   The  former  is  a  Heise  gauge  reading  up  to  600  psi, 
and  was  calibrated  with  a  dead  weight  tester.   The  latter  gauge  reads 
up  to  6000  psi  and  is  used  as  an  indicator.   DWG  is  a  Ruska  Instrument 
dead  weight  pressure  tester  and  gauge;  its  pressure  range  is  0.4  to 
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2 

170  kgm/cm  .   PI  is  a  piston-type  pressure  intensifier  manufactured 

by  High  Pressure  Co. ;  it  can  produce  pressures  up  to  10,000  psi. 


6.2   Experimental  Procedure 

The  materials  used  were  spectrograde  quality  benzene,  n-octane, 
and  carbon  tetrachloride  of  certified  purity  +99  mole  %.   H„  and  N„ 
gases  were  prepurified  grade  99.95%.   The  densities  of  the  liquid 
samples  were  determined  by  measuring  the  volumes  in  a  standard  5  ml 
pycnometer  at  25%,  and  then  measuring  the  mass  on  a  Sartorius  balance. 

The  liquid  sample  was  transferred  into  a  clean  dry  glass  tube 
and  attached  by  clear  rubber  tubing  at  stopcock  E.   The  clean 
experimental  capillary  tube  (ECT) ,  with  a  small  stainless  steel  ball, 
2  mm  diameter,  dropped  into  it, was  mounted  on  the  compression  block, 
and  connected  at  A.   C  was  closed,  with  clean  mercury  above  it.   The 
system,  with  E  closed,  was  evacuated  for  several  hours,  the  ultimate 
vacuum  obtained  was  about  50  u  Hg.   At  this  point,  E  was  opened  and 
air  in  contact  with  the  liquid,  together  with  some  liquid,  exhausted. 
E  was  closed  and  the  system  evacuated  again.   On  attaining  the 
ultimate  vacuum,  H  was  closed  off  and  a  liquid  N„  flask  held  around  the 
ECT.   On  opening  E  some  of  the  liquid  volatilized  into  the  ECT  which 
was  at  liquid  nitrogen  temperature.   E  was  closed  and  the  system 
partially  evacuated  with  the  liquid  nitrogen  flask  in  place.   D  was 
then  closed,  and  C  opened,  dropping  some  mercury  onto  the  frozen 
liquid  sample.   The  flask  of  liquid  nitrogen  was  quickly  removed.   With 
the  sample  still  frozen,  the  compression  block  was  disconnected  and 
mounted  upright  on  its  base  on  the  pressure  section.   Valves  9,  11  and 


2  were  closed  and  the  pressure  system  evacuated  with  the  vacuum  puiap. 
When  the  lower  section  of  the  frozen  nercury  column  began  to  thaw, 
pumping  was  stopped,  valves  1  and  12  shut  and  2  opened.   Mercury  from 
the  reservoir  (MR)  was  forced  up  into  the  system  by  pressurii:lng 
with  compressed  air,  usually  the  level  was  raised  up  to  the  middle 
copper  probe  of  the  MLI.   Valve  2  was  then  shut.   The  system  was  now 
loaded  and  ready  for  measurement.  The  vapor  jacket  (VJ)  war?  mounted 
on  its  base,  around  the  ECT,  and  connected  to  the  cartesian  manostat 
(CM)  through  a  separation  tube.   The  manostat  controlled  the  pressure 
in  the  vapor  jacket  to  within  i  .1  mm  Hg.   The  bath  heater  was  turned 
on,  and  the  temperature  stabilized  in  the  vapor  jacket  in  about  30-40 
minutes.   Ethylene  glycol  was  used  as  the  fluid  in  the  vapor  bath. 
The  bath  temperature  was  measured  at  the  upper  and  lower  ends  of  the 
vapor  jacket. 

For  system  pressures  up  to  8  atmospheres,  the  mercury  level 
was  read  in  the  sight  glass  SG.   For  higher  pressure  operation,  it  was 
shut  off  through  valves  4  and  7.   Pressures  up  to  600  psi  gauge  were 
read  on  the  Heise  gauge.   Higher  pressures  were  measured  with  the  Ruska 
dead  weight  tester.   The  system  pressure  was  adjusted  by  introducing 
compressed  nitrogen  gas  from  the  cylinder  at  desired  delivery  pressures 
up  to  1500  psi.   By  manipulating  the  pressure  intensifier  (PI) ,  the 
level  of  mercury  in  the  MLI  was  maintained  quite  constant  as  indicated 
by  the  lamps  of  the  MLI.   The  level  of  mercury  in  the  ECT  was  read  on 
the  Wild  Herbregg  cathetometer . 

After  completing  measurements  on  the  liquid  sample,  the 
pressure  was  released  through  valve  9,  and  the  vapor  jacket  removed. 


Two  small  metal  cups  were  mounted  on  the  ECT.   Each  had  a  hole  in 
the  bottom  somewhat  larger  than  the  diameter  of  the  ECT  with  retaining 
0-rings.   Both  cups  CI  and  C2  were  placed  lower  than  the  mercury 
level  in  the  ECT.   Liquid  N  was  poured  into  C2,  the  lower  cup.   When 
the  mercury  froze,  valve  1  was  opened,  11  shut  and  the  system  evacuated 
through  valve  12.   The  frozen  mercury  at  CI  kept  the  liquid  sample 
isolated  from  the  vacuum.   After  evacuating  to  500  y  Hg,  valves  12, 
4,  3  and  6  were  shut,  and  valve  11  opened  to  introduce  gas  into  the 
system  up  to  C2.   Valve  3  was  shut,  4,  5  and  12  opened  and  the  system 
evacuated  once  more.   Finally,  after  evacuating  at  500  p  Hg  for  10 
minutes,  12  and  1  were  shut,  2  opened,  and  mercury  from  the  reservoir 
raised  up  into  the  system  and  the  ECT,  sealing  off  the  gas  between  it 
and  the  mercury  frozen  at  C2.   Valve  2  was  then  closed  and  liquid 
nitrogen  poured  into  CI  freezing  the  mercury  in  its  vicinity;  this 
mercury  was  kept  frozen  throughout  the  following  operation.   Nitrogen 
gas  from  the  cylinder  forced  the  mercury  further 
up  into  the  ECT  below  the  trapped  gas.   There  was  also  a  visible 
mercury  head  in  SG.   The  mercury  at  C2  was  then  allowed  to  thaw;  as 
it  thawed  and  warmed  to  room  temperature,  the  gas  expanded  to  about 
twice  its  original  volume  and  pushed  the  mercury  downward.   The  system 
pressure  was  adjusted  such  that  the  mercury  level  was  just  visible 
above  the  compression  block  CB.   The  length  of  the  gas  column,  the 
mercury  levels  in  the  ECT  and  SG,  and  the  pressure  at  HG  and  TIG  were 
noted.   This  gave  a  measurement  of  the  gas  volume  at  room  temperature. 

Valves  4  and  5  were  then  shut  and  the  mercury  at  CI  allowed 
to  thaw.   Cups  CI  and  C2  were  removed,  and  the  system  pressurized  up 
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to  100  psig,  as  read  on  the  Heise  gauge.   The  mercury  column 
separating  the  gas  and  liquid  vas  broken  by  coving  the  steel  ball 
(placed  inside  initially)  through  it  with  a  poverful  magnet.   The 
two  components  were  mixed,  and  stirred  with  the  steel  ball  for  about 
10  minutes.   The  system  was  again  ready  for  P-V-T  measurements. 

The  system  was  pressurized  in  regular  increments  of  20  psi, 
with  vigorous  stirring  between  pressure  changes.   The  bubble  pressure 
of  the  system  was  reached  when  the  last  bubble  of  gas  disappeared  on 
stirring.   By  reducing  and  increasing  the  pressure  around  the  bubble 
pressure,  it  could  be  reproduced  to  within  ^  1  psi  for  all  systems 
except  the  CCl,  system  for  which  this  range  vas  ^  2  psi.   When 
reducing  the  pressure,  the  bubble  pressure  was  indicated  by  the 
appearance  of  the  first  gas  bubble.   For  pressures  above  the  bubble 
pressure,  P-V-T  measurements  were  conducted  as  described  earlier  for 
the. liquid  system. 

The  mixed  solvent  could  not  be  volatilized  into  the  vacuum 
manifold,  due  to  the  differences  in  volatility  of  the  components. 
It  was  made  up  in  the  desired  composition  by  mixing  measured  amounts 
of  the  pure  components  in  a  clean  measuring  bottle,  and  then  injected 
into  the  ECT  through  a  syringe.   The  ECT  was  mounted  on  the  compression 
block,  and  then  connected  at  A.   The  liquid  nitrogen  flask  was 
attached,  and  the  system  evacuated  while  the  liquid  mixture  was 
frozen.   B  was  then  closed,  the  nitrogen  flask  removed,  and  the 
mixture  thaxi?ed.   The  freeze  and  thaw  procedure  was  repeated  until  no 
gas  bubbles  appeared  on  thawing.   No  losses  of  material  could  be 
detected  by  the  cathetometer  in  this  sequence  of  operations.   After 
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the  last  thawing  operation,  in  all  cases,  the  level  of  liquid  in  the 
ECT  was  the  same  (as  detected  by  the  cathetotneter)  as  when  first 
introduced  into  the  ECT.   The  procedure  from  this  point  onwards  was 
the  same  as  for  the  pure  liquid  solvents. 

The  ECT  was  calibrated  on  the  vacuum  section.   The  limb  on 
the  right  of  joint  K  was  swivelled  through  180°,  so  that  the  ECT 
was  mounted  upright.   The  system  was  evacuated  to  the  ultimate  vacuum. 
The  mercury  level  was  raised  up  in  the  limb  CL  until  it  stood  a 
little  below  C.   The  system  was  pumped  out  for  another  2  hours. 
Stopcock  J  was  then  closed  and  the  mercury  raised  up  all  the  way 
into  the  ECT  with  the  levelling  bulb.   A  weighing  bottle  was  held 
under  L;  L^  and  L.  were  connected  and  a  small  amount  of  mercury 
collected  in  the  weighing  bottle.   The  amount  of  mercury  was  weighed 
and  its  value  determined  from  the  known  density.   The  mercury  level 
in  the  ECT  was  noted.   The  procedure  was  repeated  in  steps  until  the 
ECT  was  completely  drained.   A  calibration  curve  was  made  of  the  volume 
of  mercury  collected  corresponding  to  a  given  level  in  the  ECT. 

6.3  Treatment  and  Analysis  of  Experimental  Data 

The  densities  of  the  liquid  samples  were  measured  at  25°C  in 
a  pycnometer.   These  are  reported  in  Table  6-1,  and  compared  with 
values  from  the  literature.   The  measured  density  of  pure  octane  is 
A%  larger  than  the  literature  value;  all  the  other  densities  agree  to 
within  .05%. 

The  mass  of  liquid  in  the  experimental  capillary  was  determined 
from  the  measured  volume  (in  situ)  and  the  known  density.   The  virial 
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TABLE  6-1 
DENSITIES  OF  PURE  AND  MIXED  SOLVENTS  AT  25 °C 


Solvent 

Compi 

onents 
2 

^1 

1.0 

Densi 

ty> 

gms 

/cc 

1 

Measured 
.8736 

Literature 

Benzene 

.8737-^ 

Benzene 

CCl^ 

0.5 

1.2443 

1.2420^ 

CCl, 

4 

1.0 

1.5850 

1.5844-^ 

Benzene 

Octane 

0.5 

.7662 

- 

Octane 

1.0 

.7047 

.7019^ 

109 


equation  of  state,  up  to  the  second  virial  coefficient,  was  used  with 
the  measured  volume  to  determine  the  moles  of  gas  added  to  the  system. 

Table  6-2  contains  a  key  indicating  the  compositions  of 
the  mixture  studied.   The  bubble  pressures  measured  for  each  gas 
composition  are  also  shown  in  this  table.   These  were  reproducible 
in  both  increasing  pressure  and  decreasing  pressure  to  ±  1  psia  for 
the  benzene-octane  systems  and  to  ±  2  psia  for  the  benzene-carbon 
tetrachloride  systems.   The  bubble  pressures  are  used  later  to 
derive  Henry's  constants  and  activity  coefficients. 

Experimental  volumes  were  measured  to  a  precision  of  ±  .25% 
and  pressures  to  +  .6%.  The  volumes  of  constant  composition  liquid 
mixtures  were  fitted  to  equations  of  the  form 

v(T,P)  =  a  +  I  +  cP  +  dP^  (6.1) 

at  each  isotherm.   A  least  squares  fitting  program  was  used  for  this 
purpose.   Wherever  possible,  the  three-parameter  equation,  with  d  =  0, 
was  used.   The  fitted  constants  a,  b,  c,  d  are  reported  in  Table  6-3. 
6  is  the  sum  of  squares  of  the  relative  error  in  representing  the 
experimental  volume  by  equation  (6.1). 

The  isothermal  compressibility  from  equation  (6.1)  is 


<  = 


b/P^  -  c  -  2dP 

~ 2  (6.2) 

a  +  b/P  +  cP  +  dP 


and 


.J: ^  [a  +  b/P  +  cP  +  dP^I^ 

P<RT  ~~~2~  (6.3) 

RT[b/P  -  c  -  2dP] 

The  uncertainty  in  determining  <   from  the  experimental  volumes  is  ±  10%. 
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TABLE  6-2 

MEASURED  BUBBLE  PRESSURES   OF  GAS-SOLVENT   SYSTEMS 


Bubble  Pressure 

Run 

atm 

System 

# 

X 

gas 

298°K 

378°K 

418°K 

N--benzene 

A2 

.0063 

14.98 

11.98 

13.89 

A3 

.0129 

31.45 

22.60 

23.28 

A4 

.01971 

49.56 

33.96 

33.08 

A5 

.0272 

— 

46.79 

44.20 

A6 

.0330 

- 

57.03 

52.96 

N--benzene 
octane 

B2 

.0052 

4.580 

5.94 

8.63 

B3 

.0112 

9.92 

11.13 

13.84 

B4 

.0191 

17.20 

18.09 

20.80 

B5 

.0267 

24.52 

24.93 

27.62 

B6 

.0379 

35.87 

35.29 

37.90 

N.-octane 

C2 

.01065 

7.08 

7.08 

8.03 

C3 

.02113 

14.16 

13.61 

14.30 

C4 

.0436 

30.22 

27.64 

27.64 

C5 

.06335 

45.20 

40.16 

39.48 

N^-CCl^ 

U2 

.0073 

11.71 

11.44 

14.02 

U3 

.0132 

21.38 

19.06 

21.24 

. 

U4 

.02078 

34.45 

29.14 

30.50 

U5 

.0278 

47.11 

37.97 

39.21 

U6 

.0331 

57.18 

45.99 

45.88 

■• 

U7 

.03906 

54.46 

53.51 

N^-benzene 
-octane 

V2 

.0073 

16.13 

12.75 

14.77 

V3 

.0138 

31.20 

22.45 

23.61 

V4 

.0209 

48.60 

33.30 

32.91 

V5 

.0280 

67.31 

44.47 

42.65 

V6 

.0384 

- 

51.07 

48.37 

H_-benzene 

X2 

.0059 

14.57 

13.89 

15.11 

X3 

.0143 

35.40 

30.91 

29.27 

X4 

.02104 

52.55 

44.79 

40.84 

X5 

.0292 

- 

61.67 

54.87 

X6 

.0323 

- 

- 

60.18 

H„-benzene 
-octane 

Y2 

.003 

7.38 

7.53 

9.81 

Y3 

.012 

29.70 

25.79 

27.05 

Y4 

.01868 

46.70 

39.57 

40.19 

Y5 

.0232 

58.43 

49.02 

49.27 

Ill 


TABLE  6-2    (Continued) 


Bubble  Pressure 

Run 

X 

gas 

atm 

System 

298'K 

378°K 

418°K 

H  -octane 

Z2 

.0113 

16.75 

12.53 

12.25 

Z3 

.02102 

31.59 

22.74 

21.51 

Z4 

.04105 

45.20 

44.25 

40.30 

Z5 

.06078 

63.31 

65.62 

58.95 
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TABLE  6-3 


VOLUMETRIC  PtlO?ERTIES  OF  CONSTANT  COMPOSITION  MIXTURES- 
COEFFICIENTS  FOR  EQUATION  (6.1) 


Run 


Al 


A2 


A3 


Bl 


B2 


B3 


B4 


CI 


C2 


C3 


Ul 


T 

°C 

10^ 
.000 

a 

b 

2 

10  c 

-  .87 

lO^d 

25 

89.420 

.001 

0.0 

105 

.004 

99.183 

.007 

-2.03 

0.0 

145 

.060 

105.780 

.153 

-3.45 

0.0 

25 

.001 

89.204 

-  .614 

-  .94 

0.0 

105 

.070 

99.159 

-1.356 

-2.33 

0.0 

145 

.002 

105.779 

-  .803 

-3.60 

0.0 

105 

.002 

98.869 

-  .311 

-2.07 

0.0 

145 

.006 

105.450 

4.379 

-3.37 

0.0 

25 

.003 

125.482 

.232 

-1.35 

0.0 

105 

.010 

139.216 

.517 

-3.26 

0.0 

145 

.050 

148.485 

.788 

-6.02 

0.0 

25 

.005 

125.085 

.249 

-1.35 

0.0 

105 

.010 

138.858 

.844 

-3.25 

0.0 

145 

.007 

148.147 

1.643 

-5.96 

0.0 

25 

.002 

124.592 

.788 

-1.29 

0.0 

105 

.010 

138.364 

2.841 

-  .03 

0.0 

145 

.004 

147.778 

2.601 

-  .06 

0.0 

25 

.004 

123.981 

1.228 

-1.30 

0.0 

105 

.010 

137.947 

1.174 

-3.27 

0.0 

145 

.000 

147.374 

2.647 

-6.06 

0.0 

25 

.020 

162.067 

.199 

-2.06 

0.0 

105 

.010 

180.132 

-  .431 

-6.25 

0.0 

145 

.000 

192.096 

.315 

-10.39 

0.0 

25 

.010' 

160.852 

.685 

-1.93 

0.0 

105 

.030 

178.697 

2.402 

-5.05 

0.0 

145 

.050 

190.739 

3.798 

-9.23 

0.0 

105 

.001 

177.679 

1.730 

-5.21 

0.0 

145 

.010 

189.600 

6.748 

-9.14 

0.0 

25 

.000 

97.070 

-  .012 

-1.07 

0.0 

105 

.050 

107.895 

-  .039 

-2.36 

.17 

145 

.000 

115.174 

.017 

-2.96 

-.01 
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TABLE  6-3    (Continued) 


Run 

T 
°C 

23 

10^ 
.000 

a 

b 

2 

10  c 

-1.08 

10"d 

U2 

96.732 

.006 

0.0 

105 

.060 

107.661 

-  .322 

-2.26 

0.0 

145 

.060 

115.546 

-9.361 

-3.59 

0.0 

U3 

105 

.005 

107.374 

1.465 

-2.13 

0.0 

VI 

25 

.030 

93.189 

.179 

-8.70 

0.0 

105 

.070 

103.406 

.393 

-1.92 

0.0 

U5 

.000 

110.436 

.042 

-3.08 

.18 

V2 

25 

.040 

92.780 

1.774 

-  .76 

0.0 

105 

.001 

103.223 

.390 

-1.97 

0.0 

145 

.004 

110.317 

.155 

-2.99 

0.0 

V3 

105 

.001 

103.032 

.669 

-1.98 

0.0 

145 

.000 

110.098 

2.818 

-2.87 

0.0 

XI 

25 

.000 

89.428 

-  .011 

-  .89 

0.0 

105 

.003 

99.187 

-  .055 

-2.04 

0.0 

145 

.003 

105.801 

-  .015 

-3.49 

0.0 

X2 

25 

.000 

89.103 

-  .007 

-  .89 

0.0 

105 

.001 

98.913 

.142 

-2.04 

0.0 

145 

.002 

105.585 

.   .109 

-3.51 

0.0 

Yl 

25 

.005 

125.536 

-  .116 

-1.56 

-.22 

105 

.004 

139.315 

-  .114 

-3.63 

-.38 

145 

.040 

148.591 

.765 

-6.02 

0.0 

Y2 

25 

.100 

125.298 

-  .497 

-1.53 

0.0 

105 

.02 

138.964 

.469 

-3.26 

0.0 

145 

.03 

148.259 

.765 

-6.04 

0.0 

Y3 

105 

.070 

138.004 

3.184 

-3.03 

0.0 

145 

.010 

147.789 

1.708 

-5.98 

0.0 

Zl 

25 

.012 

162.057 

.304 

-2.05 

0.0 

105 

.007 

180.143 

-  .495 

-6.33 

-1.36 

145 

.000 

192.128 

.086 

-10.51 

.14 

Z2 

25 

.001 

160.616 

.816 

-1.97 

0.0 

105 

.005 

178.437 

2.003 

-5.13 

0.0 

145 

.030 

190.428 

3.299 

-9.26 

0.0 

Z3 

105 

.006 

177.033 

7.173 

-4.93 

0.0 

145 

.000 

189.102 

6.612 

-9.07 

0.0 

11: 


The  pressure  range  of  compression  dictates  the  accuracy  with  V7hich  < 
can  be  determined  experimentally.   The  presenL  experiment  was  originally 
designed  for  compressions  up  to  200  atmospheres  but  due  to  repeated 
fracture  of  the  glass  capillary  tubing,  the  maximum  operating  pressure 
was  limited  to  65  atmospheres. 

Volumes  and  compressibilities  calculated  from  (6.1)  and  (6.2) 
are  reported  in  Table  6-4,  and  compared  with  corresponding  values 
taken  from  the  literature.   At  room  temperature,  the  volumes  in  the 
benzene  and  carbon  tetrachloride  systems  agree  to  _  .5%  with  the 
literature  value.   The  compressibilities  agree  to  within  t   10%.   In 
the  benzene-carbon  tetrachloride  systems,  however,  the  measured 
compressibility  of  the  mixture  is  less  than  that  of  either  pure 
component.   At  the  higher  temperature,  the  experimental  compressibility 
of  octane  is  40%  larger  than  that  reported  by  Connolly,  and  the  volume 
1.5Z  smaller.   The  literature  values  will  be  used  in  further  analysis 
of  the  benzene-octane  system  at  this  temperature. 

Partial  molar  volumes  of  the  components  at  constant  pressure 
and  temperature  were  determined  from  the  composition  dependence  of  the 
volumes.   The  experimental  volumes  were  fitted  to  equations  of  the  form 

V  -  V  +  bx  +  ex  (6.4) 

where  x  is  the  solute  mole  fraction,  and  v  is  the  volume  of  the 
solute-free  solvent.   This  equation  was  used  only  at  conditions  x-zhere 
the  volumes  of  at  least  4  compositions  including  the  pure  solvent 
were  known.  At  low  pressures,  hence  fewer  compositions,  the  coefficient 
b  was  estimated  from  the  simplified  relation  setting  c  equal  to  zero. 
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Expression  for  the  partial  molar  volume  at  infinite  dilution, 
and  the  excess  volume  at  composition  x  are  obtained  from  equation 
(6.4).   Those  are 

v°  =  v°  +  b  (6.5) 

V  =  cx^  (6.6) 

The  uncertainty  in  the  fitted  value  of  v^  is  of  the  same  order  as 
that  of  V,  i.e.,  i  .4  cc/mole.   The  corresponding  uncertainty  in  b, 
however,  is  much  less,  "t   .11  cc/mole.   b  is  the  unit  change  in  volume 
of  the  solvent  on  addition  of  the  first  incremental  amount  of  solute, 
and  the  uncertainties  in  the  initial  and  final  volumes  are  correlated. 
The  experimental  partial  molar  volumes  should  be  accurate  to  i  .5 
cc/mole.   The  accuracy  of  estimating  the  excess  volume  parameter,  c, 
is  dependent  on  the  highest  solute  composition  of  the  liquid  mixture. 
A  typical  composition  of  0.2  mole  fraction  of  solute  will  yield  an 
excess  volume  parameter  accurate  to  1  5%.   In  the  experiments  reported 
here,  the  maximum  composition  attained  was  .06  mole  fraction  for  one 
system.   At  best,  only  an  estimate  of  c  can  be  obtained  from  this  data. 

The  calculated  partial  molar  volumes  are  shown  in  Table  6-5 
and  compared  with  values  obtained  from  the  literature.   At  each 
temperature,  the  partial  molar  volumes  are  shown  at  the  lowest  and 
highest  pressures  at  which  they  were  obtained. _  In  all  cases,  the 
experimental  v^  agrees  with  the  literature  value  to  within  the  reported 
error. 

The  last  column  in  this  table  shows  partial  molar  volumes 
determined  from  the  generalized  correlations  (4.3)  and  (4.6).   The 
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TABLE 

6-5 

PARTIAL 

MOLAR  VOLLTEES  OF 

SOLUTES 

cc 

^1 

,/mole 

System 

T 

P 
atm 

Expt. 

Literature 

Correlations 

N„-benzene 

298 

21.4 

48.1 

52.6' 

45.3 

2 

55.5 

46.7 

- 

43.2 

378 

14.6 

73.8 

- 

68.8 

55.5 

69.0 

- 

65.1 

418 

21.4 

92.7 

- 

90.1 

55.5 

88.3 

- 

87.1 

N„-benzene 

298 

21.4 

49.4 

- 

46.2 

-octane 

55.5 

48.6 

- 

44.5 

378 

14.6 

73.2 

- 

66.2 

55.5 

72.7 

- 

63.7 

418 

21.4 

92.7 

- 

86.8 

55.5 

89.1 

- 

82.7 

N„-octane 

298 

21.4 

53.7 

— 

48.2 

2 

55.5 

53.2 

- 

47.4 

378 

21.4 

77.2 

- 

67.1 

55.5 

75.7 

- 

65.7 

418 

21.4 

95.4 

- 

87.3 

55.5 

91. 

- 

84.0 

H^-benzene 

298 

41.8 

35.3 

35.0 

+ 

2^ 

33.0 

2 

55.5 

35.1 

- 

31.2 

378 

35.0 

52.7 

52.2 

+ 

39 

45.1 

55.5 

51.3 

49.3 

+ 

39 

47.7 

418 

35.0 

65.9 

66.6 

+ 

39 

67.9 

55.5 

64.3 

63.4 

+ 

39 

64.4 

H„-benzene 
-octane 

298 

35.0 

34.1 

- 

33.2 

55.5 

32.6 

- 

32.1 

378 

35.0 

52.7 

- 

50.6 

55.5 

49.3 

48.2 

418 

28.2 

67.0 

- 

64.3 

55.5 

65.0 

- 

61.2 

H„-octane 

298 

35.0 

38.0 

35.9 

-t- 

3^ 

35.1 

2 

55.5 

37.6 

35.9 

+ 

3^ 

34.2 

378 

28.2 

57.2 

56.0 

+ 

50.3 

55.5 

53.6 

53.8 

+ 

39 

48.2 

418 

28.2 

69.3 

75.5 

+ 

$ 

66.3 

55.5 

66.2 

69.3 

+ 

62.3 
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TABLE  6-5    (Continued) 


-o 

^1 

T 

P 

atm 

cc/mole 

System 

Expt. 

Literature 

Correlations 

N^-CCl^ 

298 

35.0 

52.0 

53.0 

2^° 

46.6 

55.5 

51.5 

- 

45.0 

378 

35.0 

78.0 

- 

71.1 

55.5 

75.1 

- 

69.0 

418 

21.4 

99.8 

— 

98.3 

55.5 

97.6 

- 

95.1 

N„-benzene 
4 

298 

21.4 

49.6 

_ 

46.2 

55.5 

48.7 

- 

44.1 

378 

14.6 

75.8 

— 

69.8 

55.5 

73.4 

- 

67.1 

418 

21.4 

95.6 

- 

94.7 

55.5 

94.3 

— 

91.5 
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mixed  solvents  were  represented  by  a  one-fluid  model  and  the 
experimental  volumes  used  to  calculate  reduced  densities.   The  resulting 
v^  is  within  15%  of  experiment,  indicating  that  the  one-fluid  model 
used  in  conjunction  with  equations  (4.3)  and  (4.6)   yields  reliable 
estimates  of  the  solute  partial  molar  volume. 

Activity  coefficient  parameters  and  Henry's  constants  for 
solute-solvent  systems  were  determined  from  the  experimental  bubble 
pressures.   The  solute  activity  coefficient  at  low  concentrations 
is  represented  by  a  one-term  Margules  expression 

A  two-term  Margules  expression  represents  the  excess  free  energy  of 
the  mixed  solvent.   The  activity  coefficients  for  the  solvent  com- 
ponents in  the  mixture  are 

Iny^^   "  ^2^  ■*"  ^^^^  "  ^^^3  "^  ^(A  -  B)X3]  (6.8) 

Any 3^  "  ^l""!  '^   ^^^^  '   ^^^2  ^  ^^^   "  ^^^^2^  ^^'^^ 

X-  and  X_  are  the  solute  free  mole  fractions  in  the  solvent  mixture. 
A  and  B  are  the  coefficients  of  the  excess  free  energy  expression  and 
are  listed  in  Table  6-6.   The  temperature  variations  in  the  excess 
free  energies  were  estimated  by  the  Gibbs-Kelmholtz  equation. 

The  vapor  phase  properties  are  characterized  by  the  Redlich- 
Kwong  equation.   The  computational  method  and  algorithm  for  the 
determination  of  Henry's  constant  and  activity  coefficient  parameters 
from  bubble  pressures  is  described  in  Appendix  E.  The  results  of  these 
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TABLE  6-6 
EXCESS  FREE  ENERGIES  OF  EQUIMOLAR  SOLVENT  MIXTURES 


T 

i 

E 

System 

298 

A 

B 
.0331 

calS; 
82 

/mole 

11 
Benzene 

.0331 

.00 

carbon  tetrachloride 

378 

.0233 

.0233 

73, 

.14 

418 

.0196 

.0196 

68, 

.16 

12 
Benzene 

343 

.85 

.11 

70, 

.84 

n-octane 

298 

1.1445 

.1481 

95, 

.38 

378 

.6208 

.0803 

65, 

.62 

418 

.3613 

.0468 

42. 

,23 
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TABLE  6-7 


HENRY'S  CONSTANTS  AND  ACTIVITY  COEFFICIENT  PARAMETERS 
FROM  EXPERIMENTAL  DATA 


l,m 

T 
°K 

298 
378 
418 

h 

atm 

System 

Experiment 

1438 

1045 

929 

Literature 

H„-octane 

-   .36 
.054 
.205 

1464-'--^ 

N„-octane 

298 
378 
418 

-  .28 
+  .142 
.346 

648 
606 
588 

673" 

N^-CCl^ 

298 
378 
418 

.452 

,911 

1.076 

1544 
1238 

1144 

1590" 

N2-^6»6 

298 
378 
418 

.502 

.874 

1.129 

2290 
1541 
1338 

2280^^ 

^2-^6«6 

298 
378 
418 

.510 

.962 

1.232 

2429 
1994 
1646 

3882^^ 

2121 

1644 

^-^6^6 

CCl. 
4 

298 
378 
418 

-1.127 

-  .330 

-  .075 

2077 
1406 
1254 

^2-^6»6 
octane 

298 

378 
418 

-  .761 

.381 

2.189 

850 
869 
959 

«2-^6»6 

298 
378 
418 

-  .684 

.519 

2.270 

2408 
2043 
2072 
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calculations  are  shown  in  Table  6-7.   There  is  an  experimental 
uncertainty  of  i  6%  in  the  Henry's  constants.   The  major  contribution 
to  this  is  from  the  uncertainty  ir  determining  the  mole  fraction  of 
the  solute.   For  systems  of  gases  dissolved  in  single  solvents,  the 
measured  Henry's  constant  agrees  with  literature  values  to  within 
these  limits.   The  sole  exception  is  hydrogen  in  benzene  at  298°K. 
The  disagr-'fc-ent  here  is  35%;  there  is  no  apparent  experimental 
factor  to  account  for  a  discrepancy  of  this  magnitude.   Both  values 
will  be  used  later  to  analyze  the  mixed  solvents  data  with  the 
expressions  of  the  distribution  function  solution  theory. 

The  experimental  f„  increases  with  increasing  temperature. 
This  is  in  accordance  with  the  corresponding  states  correlations 
discussed  in  the  next  chapter.   Some  negative  values  of  f_  are  encountered 
In  the  mixed  solvents  systems,  indicating  a  positive  slope  for  the 
solute  activity-coefficient  at  zero  composition.   However,  this  is  not 
inconsister.t  with  the  direct  correlation  function  solution  theory. 
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CHAPTER  7 
SOLUTION  THEORY  FOR  GAS-LIQUID  SYSTEMS 

The  expressions  for  thermodynamic  properties  obtained  from  the 
solution  theory  are  most  generally  applicable  to  a  system  of  one  or 
more  gases  in  a  multicomponent  solvent.   This  system  is  characterized 
by  the  complexities  of  liquid-liquid  and  gas-liquid  types  of  interactions. 
Those  methods  previously  found  successful  in  describing  these  interactions 
individually  are  combined  to  analyze  the  complex  multicomponent  system. 

Generalized  correlations  of  the  nonideality  parameters  for 
gases  in  single  solvents  are  investigated  first.   These  provide  a 
quantitative  description  of  the  gas-solvent  interaction.   The  hydrogen- 
benzene  system,  for  which  complete  thermodynamic  data  is  available,  is 
studied  to  examine  the  validity  of  expressions  for  the  composition 
dependent  thermodynamic  properties,  and  the  characteristics  of  triplet 
correlation  functions  in  gas-liquid  systems  as  opposed  to  those  in 
liquid-liquid  systems.   The  experimental  DCF  integrals  in  ternary 
systems  are  next  studied  for  their  composition  and  macroscopic  state 
dependence.   In  addition  to  the  usual  thermodynamic  properties,  these 
quantities  determine  the  relations  between  the  Henry's  constants  of 
the  solute  in  the  mixed  and  pure  solvents,  respectively.   Finally,  the 
Henry's  constants  and  activity  coefficients  in  the  mixed  solvent  are 
determined  from  individual  gas-solvent  parameters  and  used  to  calculate 
vapor-liquid  equilibrium  in  the  ternary  mixture. 

7.1  Correlations  for  Activity  Coefficient  Parameters 

A  simple  representation  of  the  unsymmetric  convention  activity 
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coefficients  in  a  binary  mixture  is 


ilnv*  =  f°(xj  -  2x^)  +  f°(xj  -|xj)  (7.1) 


and 


iny^   =  f°x2  +  f°xj  (7.2) 


of  the  pure  solvent  at  system  conditions.   C. , ,  which  represents  the 


In  terms  of  DCF  integrals,  f^  is  given  by 

-2f°=  (1-C°,) f-  (7.3) 

^^    (1  -  C°2) 

Orentlicher  has  used  expressions  (7.1)  and  (7.2),  with  f^  =  0,  to 
represent  the  activity  coefficients  for  hydrogen  in  croygenic  solvents. 
He  finds  the  coefficient  f°  to  vary  inversely  with  absolute  temperature, 
Equations  (A. 3)  and  (4.6)  relate  C^^   and  0°^  to  the  density 
pure  solvent  at  system  conaitions.   ^iii 
self-interaction  of  solute  molecules  in  the  solvent  medium,  is  also 
likely  to  depend  on  the  solvent  density.   At  this  point  one  would 
attempt  to  correlate  C^ ^  as  a  function  of  reduced  solvent  density. 
However,  a  close  examination  of  equation  (7.3)  shows  that  the  accuracy 
required  in  correlating  C^ ^  is  of  the  same  order  as  the  required 
accuracy  of  f°.   For  the  experimental  systems  studied  by  Orentlicher, 
f°  varies  from  0.0  to  1.0.   At  low  solute  compositions,  an  error  of 
10%  in  the  activity  coefficient  corresponds  to  an  equivalent  error  in 
f°,  i.e.,  about  1  0.1  absolute  error.  This,  in  turn,  requires  an 
accuracy  of  i  .2  in  the  absolute  value  of  C^,  which  has  a  large 
negative  value. 

Isothermal  vapor-liquid  equilibrium  data  were  used  to  determine 
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experimental  values  of  f„  and  the  solute  Henry's  constant.   The 
activity  coefficients  for  the  solute  and  solvent  were  represented  by 
equations  (7.1)  and  (7.2),  respectively.   Vapor  phase  nonidealities 
were  characterized  by  the  Redlich-Kuong  equation.   The  details  of  the 
calculations  and  the  relevant  equations  are  presented  in  Appendix  E. 

The  experimental  data  for  hydrogen  systems  were  those  studied  by 

2 
Orentlicher,  and  the  hydrogen-benzene  data  of  Connolly.    Data  on  the 

3  A 

methane  systems  were  taken  from  Sage   and  from  Shim.    The  fitted 

parameters  for  these  systems  are  shown  in  Table  7-1.   The  calculated 

f„  shows  a  significant  variation  with  temperature,  or  reduced  solvent 

density  for  all  systems.   In  the  methane  systems   the  contribution 

of  f.  to  the  activity  coefficients  is  significant.   Under  these 

conditions  a  proper  representation  of  the  activity  coefficient  requires 

correlations  for  f„  and  f„. 

The  next  step  was  to  search  for  an  effective  value  of  f- 

which  would  be  an  average  of  the  actual  f„  and  f..   This  was  achieved 

(JlnY*)^''^  xj  -  3/2  x^ 

by  plotting  — ^ versus  — and  selecting  the  best 

x^  -  2x  x^  -  2x, 

horizontal  line  to  fit  the  data  over  the  composition  range.   A 

horizontal  line  implies  that  f^  =  0.   The  zero  intercept  of  the  function 

eff  eff 

gave  f-   .   The  extracted  values  of  f„   for  the  hydrogen-solvents 

system  and  methane-solvents  systems  are  shown  in  Figs.  7-1  and  7-2, 

'         .       eff 
respectively.   There  is  an  approximate  linear  relationship  for  f- 

with  solvent  density.   For  the  methane  systems   it  is 

£2"  =  4.355  -  1.570  p°  (7.4) 

and  for  hydrogen  systems 
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Fig.  7-1.   Effective  f  for  H  -solvents  systems. 
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Fig-  7-2.   Effective  f.  for  CH. -solvents  systems. 
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f^^^  =  5.182  -  1.679  ^2  (7-5) 

The  validity  of  these  approximations  was  tested  by  comparing  calculated 
bubble  pressures  and  vapor  compositions  for  some  systems  with  experimental 
data.   The  results  are  shown  in  Table  7-2  with  calculations 
for  some  systems  not  included  in  the  formulation  of  correlations 
(7.4)  and  (7.5).   The  calculation  of  bubble  pressures  is  similar  to 
that  used  in  the  fitting  of  parameters.   There  is  good  agreement 
between  calculated  and  experimental  vapor  compositions.   This  agreement 
does  not  depend  on  a  good  calculation  of  the  bubble  pressures   and  may 
indicate  an  insensitivity  to  the  equations  describing  the  nonideality. 
The  calculated  bubble  pressures  agree  well  with  experimental  data  in 
systems  of  moderate  composition.   At  larger  solute  concentrations   the 
pressures  may  differ  by  20%  from  experiment.   The  calculations  for 
methane  in  hexane  and  hydrogen  in  butane  indicate  that  for  systems  in 
which  f^   from  equation  (7.4)  differs  considerably  from  the  experimental 
value,  the  maximum  error  in  the  bubble  pressure  is  15%. 

The  sensitivity  of  calculated  component  mole  fractions  to  the 
activity  coefficient  relation  is  an  important  consideration  in  design 
calculations.  (A  method  of  operation  of  separation  equipment  is  to 
control  the  pressure  and  heat  input  to  a  unit  and  allow  the  stream 
compositions  to  adjust  themselves  to  equilibrium  values.)  The 
calculated  compositions  for  some  representative  systems  are  reported 

in  Table  7-3.   Comparisons  are  shown  with  experimental  results  and  the 

10 
predictions  of  the  activity  coefficient  relation  of  Chueh  and  Prausnitz. 

This  is  a  widely  used  representation  for  liquid  phase  activity  coefficients. 
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They  recommend  a  modified  Van  Laar  expression.   To  correct  for  the 
size  differences  between  solute  and  solvent  molecules,  they  use  the 
solute  volume  fraction  rather  than  the  simple  mole  fraction  as  the 
composition  variable.   For  many  systems   this  equation  predicts  activity 
coefficients  in  good  agreement  with  experimental  data  up  to  about  0.3 
solute  mole  fraction.   However,  the  model  parameters  are  characteristic 
for  each  solute-solvent  interaction  and  must  be  previously  determined 
from  experimental  data.   A  comparison  of  equation  (7.1)  with  Chueh's 
equation  is  meaningful  in  light  of  the  loss  of  accuracy,  if  any,  that 
is  incurred  in  going  from  a  generalized  parameter  equation  such  as 
equations  (7.4)  and  (7.5),  to  a  nongeneralized  equation  to  calculate 
vapor  liquid  equilibrium.   The  calculated  vapor  compositions,  in 
Table  7-3   are  uniformly  good  and  agree  consistently  with  experimental 
data.   Both  methods  yield  an  equally  accurate  liquid  composition  for 
the  hydrogen  system.   In  the  methane  systems   the  liquid  mole  fractions 
calculated  by  equation  (7.4)  do  not  agree  as  well  with  experimental  data 
as  do  those  from  the  Chueh  equation.   The  deviations  at  the  highest 
pressures  are  up  to  20%  for  equation  (7.4)  in  the  methane-propane  system 
and  10%  in  the  methane-hexane  system.   Therefore,  the  use  of  Chueh's 
equation  for  systems  of  known  parameters  yield  uniformly  good 
predictions  of  equilibrium  composition  at  fixed  temperature  and  pressure. 
At  moderate  compositions   equations  (7.4)  and  (7.5)  which  need  only 
one  generalized  parameter  for  each  component  will  yield  comparable  results. 

The  DCF  integral  C^ ^  is  more  significant  than  the  quantity  f„ 
in  understanding  the  molecular  correlations  in  the  fluid  mixture.   The 
qualified  success  in  correlating  activity  coefficients  with  f   suggests 
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that  correlations  of  C^  for  individual  gas-solvents  systenPare  possible. 
Experimental  values  of  C^   for  hydrogen-solvents  and  methane-solvents 
were  obtained  through  equation  (7.3).   C„„  and  C^^  were  calculated  from 
equations  (4.3)  and  (4.6),  respectively.   The  experimental  values  of 
fj   were  those  obtained  from  fitting  bubble  pressures  and  composition  and 
are  reported  in  Table  7-1. 

C^ ^  for  the  methane  systemswas  correlated  by  the  equation 


£n 


-  C 


11 


*  ^ 

* 


=  .36856  +  .9458  p, 


(7.6) 


The  corresponding  equation  for  hydrogen  systems  and  one  nitrogen-decane 
system  is 


In 


-  C 


11 


*  ^ 

* 


=  -  3.049  +  1.92606  p, 


(7.7) 


The  accuracies  with  which  the  equations  represent  experimental  activity 
coefficients  are  shown  in  Table  7-4.   f„(exp)  is  the  actual  representation 
of  experimental  data.   It  is  compared  with  the  correlations  of  f_, 
equations  (7.4)  and  (7.5),  and  the  values  of  f„  from  the  correlation 
of  C^ ^ ,  equations  (7.6)  and  (7.7).   The  correlations  for  C^  generally 
provide  a  value  of  f„  in  closer  agreement  with  the  experimental,  and  the 
effective  f.,  than  the  correlation  of  f„  itself.   Under  these  circumstances, 
calculations  of  binary  vapor-liquid  equilibrium,  using  f„  from  a  C 
correlation,  will  agree  even  more  closely  with  experimental  data  than  from 
an  effective  f„.   In  the  minority  of  systems,  where  the  f„  correlations 
are  better,  e.g.,  methane-butane  at  310°K,  f„(exp),and  f„ (effective) , 
is  itself  small,  and  moderate  errors  in  its  prediction  will  have  little 
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effect  on  the  calculated  activity  coefficient.  The  correlations  for 
C^,  equa;:ions  (7.6)  and  (7.7),  are  preferred  over  equations  (7.4) 
and  (7.5)  for  representing  activity  coefficients.   Their  importance 
lies  in  that  they  provide  a  state  dependence  of  the  DCF,  C^^,  in 
binary  systems. 


7.2  Detailed  Analysis  of  Hydrogen-Benzene  Systen 

There  is  a  twofold  interest  in  studying  this  system.   On  the 
one  hand,  complete  experimental  data  of  high  accuracy  is  available  to 
derive  the  composition  dependence  of  the  DCF  integrals.   Secondly,  the 
predictions  of  the  macroscopic  correlations  derived  earlier  can  be  tested 

with  experimental  data. 

2  12  13 

Connolly  and  coworkers  '   '   have  made  a  systematic  investiga- 
tion of  the  thermodynamic  properties  of  this  system  over  a  wide  range 
of  temperature  and  pressure.  They  provide  analytical  representation 
of  the  following  experimental  relations: 

(i)    Volume  of  pure  benzene  with  temperature  and  pressure, 

(ii)   Volume  of  hydrogen-benzene  mixtures  with 
temperature,  pressure  and  composition, 

(iii)   Activity  coefficient  with  temperature  and  composition, 

(iv)    Henry's  constant  with  temperature. 
These  expressions  are  sufficient  to  extract  the  composition  dependent 
DCF  integrals  from  the  following  expressions 


PkRT 


=  l-4^U-x^C22-  2x^X2C^2  ^^'^^ 


mix 
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pw^   =  (1  -  x^C^^  -  ^^C^^)/a/D<-BI)^^^  (7.9) 


*  ^  2 


3x^ 


x„  (l/p<RT)  . 

(7.10) 


where  hydrogen  is  component  1.  The  integrals  ware  determined  from  the 
experimental  data  at  selected  temperatures  and  pressures.   The  set  for 
SIS'K  are  shown  plotted  in  Fig.  7-3.   The  thermodynamic  properties  vary 
rapidly  with  composition  at  the  higher  temperatures,  and  the  DCF 
integrals  also  vary  considerably  with  composition.  An  assumed  linear 
composition  dependence  will  not  reproduce  the  experimental  quantities. 
However,  such  an  assumption  would  not  be  useful  since  the  experimental 
values  of  C„„,  C  ,  and  C^  are  not  usually  available  in  the  mid- 
composition  range. 

In  systems  having  only  one  limiting  composition  which  is 
physically  realizable,  series  expansions  about  this  composition  provide 
a  method  for  describing  thermod>'namic  properties  at  other  compositions. 
The  solution  theory  expressions  for  the  thermodynamic  properties  are 


1  -  C°  +  x,[2C°  -  2C°,  -  C?,,  -  C°  (1  -  pv°)  +  C°        pv°] 


p<RT        22    1'   22     12    122    22'    "^    V  222  ^  1" 

+  ...  (7.11) 

,0   o     o— O  ,0    O   »   0—0  o  ,,   0— o«   o  ,    o— o, 2, 

-     o-o  ,     ^-SrSl2^P2^l(^S2^-S22^-P2-^lS2(^^2"l^-S22(P2-^l^  ^ 
pv.  -  p,V.  +  X 

^   ^  ^   '  (1  -  c^p 

+  ...  (7.12) 

Any*  =  f2(2x^  -  xj)  +  f°(xj  -  |  x^)  (7.13) 
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The  coefficients  in  these  expressions  are  DCF  integrals  at  zero 
composition  and  can  be  obtained  from  the  correlations  derived  earlier; 
the  third  order  functions  are  obtained  froci  the  density  derivatives 
of  the  pair  functions. 

The  expressions  for  compressibility  and  partial  molar  volume 
are  linear  in  the  hydrogen  mole  fraction.   Consequently,  these 
equations  cannot  predict  any  quadratic  or  higher  order  composition 
dependence.   However,  it  is  possible  to  include  higher  order  composition 
dependence  by  assuming  an  expression  for  the  excess  volume  in  the 
mixture.   The  simplest  such  expression  in  the  asymmetric  convention  is 


v^  »  Ax^  (7.14) 


This  leads  to 


or 


-o  ^  -E 
V        V   +  V 

pv,  -  — =•  (7.15) 

^  1   V      -o  ,    -o  ,   E 

X- V   +  X_V   +.  V 

P^l S^TT ,  o,  2  ^^  •  ^^^ 

^l'^2"^l  "^  ^2   ^2  1 

which  is  a  complete  composition  dependent  expression  for  pv^ .   The 
coefficient  A  is  related  to  the  slope  of  the  partial  molar  volume  at 
zero  composition.   From  (7.16) 

r 


9pv^ 
3x^ 


=  2Ap°  +  p°v°  -  (p°v°)2  (7.17) 


>T,P,x^-0 

A  comparison  of  equations  (7.12)  and  (7.17)  now  yields  A  in  terms  of  the 
DCF  integrals.   The  simplified  expression  is 
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P^^ 


[-(C°^  -.  C°^2>  ^   2p°v°(C°2  -.  C°,p 


(P^V°)^C°,,.C°,)1 


2(1  -  q^) 


(7.18) 


An  expression  for  the  isothermal  compressibility  is 


p<RT 


-o     -o  ,  .  2v2 
(x^v  +  x-v  +  Ax  ; 

RT(x^v°'  +  X2V°'  +  A'xJ) 


(7.19) 


where  the  primes  denote  partial  derivatives  vj.r.t.    pressure  at 
constant  temperature.   This  expression  when  rearranged  and  truncated 
to  second  order  in  x^ ,  yields 


1 

p<RT 


,    ,    o-o  ^   o.  2,2 
(x^  +  ^iP2'^i  +  P2^'^i^ 

s    '  '  '     ■—-■——-■■■--  ■■■.■■     .  ■  ■■    .1^  — 

-Oi 

o  o„„,    ,    ^1  _^  A'    2, 

p  K-Rj  (x  +  x  — —  +  — -  X  ) 

2.   Z  I  1  -ot    -oi   1 

V       V 

2     2 


(7.20) 


where 


-Ot 

-Oi 

^2 


^122  +  ^  "  P2^1^S22  ''  ^2^ 

(1  -  c°p 


(7.21) 


A'  denotes  the  pressure  derivative  of  the  excess  volume  and  is  usually 
small. 

Equations  (7.20)  and  (7.16)  were  used  in  calculation  of  the 
isothermal  compressibility  and  the  partial  molar  volumes.   The  pair 
DCF  integrals  were  obtained  from  the  correlations  (4.3),  (4.6)   and 
(7.7).   The  triplet  function:  C°„,  C°  _>  ^^^^  ^\Vl  ^^^^   obtained  from 
the  density  recurrence  relations 


pO    _   o 

^ij2  -  P2 


9P° 


-  c 


ij 


i,j-l,2 


(7.22) 


T,P 
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Lastly,  p  A  was  evaluated  from  equation  (7.18)  and  A'  was  set  equal 
to  zero.   The  calculated  compressibilities  and  partial  molar  volumes 
are  shown  in  Figs.  7-4  and  7-5.   Each  figure  shows  two  sets  of 
calculations.   The  solid  lines  represent  calculations  wherein  both 
terms  of  equation  (7.22)  were  evaluated  from  the  correlations.   They 
reflect  the  abilities  of  the  correlations  to  predict  these  thermo- 
dynamic properties  from  a  minimum  of  experimental  data.   The  only 
experimental  number  used  in  this  set  of  calculations  is  the  molar 
volume  of  pure  benzene  at  the  system  temperature  and  pressure.   There 
is  good  qualitative  agreement  at  all  state  conditions,  and  the 
composition  dependences  show  the  same  trends  as  the  experimental  data. 
But  the  errors  in  predicting  the  infinite  dilution  values  for  both 
properties  is  a  limiting  factor  in  this  calculation. 

To  overcome  this  restriction,  a  set  of  calculations  was  made 
using  the  experimental  values  of  C   in  equation  (7.22),  the  other  term 
again  being  obtained  from  the  correlations.   These  calculations  are  repre- 
sented by  the  broken  lines  in  the  figure.   As  required,  they  coincide  with 
experimental  values  at  infinite  dilution  but  the  slopes  at  zero  composi- 
tion are  too  small  (in  absolute  value).   This  error  can  arise  from 
erroneous  values  of  C„„„,  C^„„  or  C^^^*   T°  examine  the  sources  of  the 
error,  experimental  and  calculated  values  of  these  quantities  are  shown  in 
Table  7-5.   The  calculated  pair  DCF  integrals  are  in  good  agreement  with 
experimental  data,  as  is  Cj„2»  which  is  the  density  (or  pressure)  derivative 
of  the  isothermal  compressibility.   But  the  calculated  C  „„  and  C   „  have 
the  opposite  sign  as  the  experimental  values   and  differ  appreciably  in 
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Fig.  7-4.   Isothermal  compressibility  and  partial  molar 
volume  in  Hj-benzene  system. 
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Fig.  7-5.   Isothermal  compressibility  and  partial  molar 
volumes  in  H  -benzene  system. 
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T^^LE  7-5 

DCF  INTEGRALS  IN  HYDROGEN-BENZDJE  SYSTEM 
AT  413° K,  50  ATM 


Experinient 
-10.69 

Calculated 

b 

a 

,o 
22 

-10.69 

-  9.10 

,o 
12 

-   6.01 

-  6.01 

-  5.84 

,o 
11 

-      .94 

-     .94 

-     .48 

,o 
'222 

-79.14 

-77.23 

-78.82 

,o 
122 

+  1.59 

-24.75 

-24.58 

o 

■1  1  o 

10.00 

-  5.00 

-  4.62 

fl  o 

C..  from  correlations  used  in  equation  (7.22). 

experimental  C.  used  In  equation  (7.22). 
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magnitude.   The  correlations  do  not  yield  a  satisfactory  value  of  the 
pressure  derivatives  of  the  partial  raolar  volume  and  the  activity 
coefficients,  although  the  properties  themselves  are  well  predicted. 
Since  C- „2  ^nd  C^ ^ „  appear  in  the  expression  for  p^A,  the  excess 
volume  in  this  system  cannot  be  satisfactorily  determined  from  the 
present  correlation  alone. 

It  is  interesting  to  note  that  neither  the  Percus-Yevick  nor 
the  MC  theories  (Appendix  F)  provides  relations  for  the  triplet 
functions  that  are  consistent  with  the  experimental  data.   The  HNC 
approximation  requires  that  all  C. .,  be  zero,  which  is  clearly  in 
contradiction  to  the  data.   For  the  Percus-Yevick  approximation,  one 

has 

C222  -  -  (^22)^  (7.23) 

C°22  =  -  (C°2)'    .  (7.24) 

The  first  relation  is  approximately  consistent  with  the  data,  but 

the  second  is  not.   Further,  no  relation  is  available  for  C^ 

112 

Refinements   °f  the  correlations  for  C  „  and  C^^  should  improve 

the  accuracy  of  the  predicted  triplet  functions.   However,  this  stage 

of  development  calls  for  examining  and  correlating  experimental  data 
which  Is  scarcely  available. 

7.3  Thermodynamic  Properties  of  Gas-Mixed  Solvents  Systems 

The  study  of  thermodynamic  properties  of  subcritical  systems 
suggests  that  liquid  mixtures  are  satisfactorily  represented  by  the  one- 
fluid  model  with  a  composition  dependent  characteristic  parameter. 
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IChen  the  mixed  solvent  in  the  multicomponent  gas-solvent  system  is 
so  represented,  the  subsequent  analysis  is  that  for  a  gas  in  a  single 
solvent. 

With  the  assumption  of  the  one-fluid  model  for  the  solvent 
mixture,  the  expressions  for  thermodynamic  properties  of  the  1-2-3 
ternary  become 

■'•       =  1  -  x?C,,    -  x^C       -  2x,x  C,  (7.25) 


PkRT  1   11  m  mm  1  m  Im 


-^  =  1  -  x.C.,    -  X  C,  (7.26) 

kRT  1  11  m  Im 

(1  -  c° )   (1  -  c°y 

f^  = 5-^^  +  ^^-  (7.27) 

2(1  -  C°   ) 
mzu 

where  the  subscripts  m  denote  mixture  quantities.   It  is  important  to 

note  that  C   and  C.   are  not  mole  fraction  averages  of  the 
mm      Im 

respective  pure  component  integrals  C^j*  ^-jo*  etc.   For  example,  in 

general 

C   ^  x_C„_  +  x.C__  (7.28) 

mm    2  22    3  33 

for  the  one-fluid  model,  although  the  equality  may  hold  in  certain 

cases.   However,  C_-  and  C--  are  the  limiting  values  of  C at  the 

22      33  inm 

extreme  compositions.   The  variation  of  C   with  composition  is 

tutti 

determined  by  the  reduced  density  of  the  mixture.   Similar  remarks 

hold  for  C,   and  C. , . 
Im      11 

Equations  (7.25),  (7.26)  and  (7.27)  were  used  to  analyze  the 
experimental  data  described  in  Chapter  6.   The  generalized  correlations 
(4.3)  and  (4.6)  were  used  to  screen  the  experimental  data  in  the  process 
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of  obtaining  DCF  integrals.   The  results  are  presented  in  Table  7-6. 

C  's  as  predicted  by  the  generalized  correlations  are  listed  as 

calculated  values.   For  most  systems,  the  quantities  C   for  the 

mm 

solvent  mixture  lie  between  those  of  either  pure  component.   This 
behavior  is  anticipated  but  not  required  for  consistency  with  the  one- 
fluid  model.   The  trend  is  not  as  marked  in  the  integrals  C,   and  is 

Im 

totally  absent  for  the  integrals  C   . 

More  significant  is  the  consistent  variation  of  all  these 
quantities  with  the  reduced  density  of  the  solvent.   They  all  increase 
positively  with  increasing  temperature  indicating  that  on  the  molecular 
scale  the  negative  contribution  to  c(r)  from  the  repulsive  part  of  the 
intermolecular  potential  decreases  faster  than  the  positive  contribution 
from  the  attractive  part  of  the  potential.   Further,  regardless  of 
the  solvent  composition,  the  integrals  involving  solute  interactions 
are.  always  less  negative  than  those  for  solvent  interactions.   This 
suggests  that  the  size  and  energy  parameters  characterizing  the  solute- 
solute  interaction  are  always  smaller  than  those  of  the  solvent- 
solvent  interaction,  and  the  solute-solvent  parameters  lie  in  between. 

For  comparison  with  experiment  the  isothermal  compressibilities 

and  partial  molar  volumes  for  the  N„-C-H,-C„H- „  system  were  calculated 

zoo   o  lo 

with  equations  (7.20)  and  (7.16).  The  calculated  composition 
dependence  is  shown  in  Figs.  7-6  and  7-7.   The  solid  lines  represent 
calculations  for  which  the  experimental  DCF  integrals  (of  Table  7-6) 
were  used,  with  triplet  functions  from  the  generalized  correlations. 
The  calculations  are  in  good  agreement  with  experimental  data  over  the 
entire  composition  range  of  the  solute.  The  dashed  lines  represent 
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Fig.  7-6.   Isothermal  compressibility  of  N_-benzene- 
octane  mixtures  at  298°K. 
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Fig.  7-7.   Partial  molar  volume  of  N2  in  N2-benzene-n-octane 
mixtures. 


calculations  perfonr.ed  with  all  DCF  integrals  (pair  and  triplet) 
obtained  from  the  generalized  correja'icns.   Tliid  only  experimental 
data  required  for  this  calculation  is  il.e  volume  of  the  solvent  mixture. 
The  experimental  compressibilities  of  che  pure  ;;olvents  are  not 
satisfactorily  represented  by  equation  (4.3),  and  this  causes  a 
considerable  error  in  calculating  the  cc:r.pressibility  of  the  solvent 
mixture,  as  indicated  by  the  difference  between  the  solid  and 
dashed  lines  at  zero  gas  composition.   Much  of  the  overall  error  can 
be  eliminated  by  adjusting  the  pure  component  characteristic  volumes 
to  obtain  agreement  between  their  calculated  and  experimental 
compressibilities.   The  mixture  characteristic  volume  is  then  revised 
to  be  the  mole  fraction  average  of  the  adjusted  pure  solvent 
characteristics.   The  resulting  compressibilities  and  partial  molar 
volumes  of  the  gas-mixed  solvent  system  are  represented  by  the  dashed- 
dotted  lines.   The  improvement  in  both  calculated  properties  is  quite 
marked,  and  brings  them  into  satisfactory  agreement  with  the  experimental 
data.   Thus  the  predictive  method  using  the  corresponding  states 
correlations  with  adjusted  pure  component  characteristics  gives  a 
good  representation  of  the  experimental  volumetric  properties. 

The  experimental  DCF  integrals  also  provide  a  relation  between 
Henry's  constants  of  component  1  in  the  2-3  mixture  and  in  the 
individual  solvents,  respectively.   From  equation  (1.9)  we  have 

M  c  ^    "     C^ 

JlnH,    -   y  X.£nH,  ,  =  f^^^   -     I     X.f^^^  (7.29) 

where  H,   and  H-  .  are  the  Henry's  constants  of  the  solute  in  the 
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mixture  and  in  solvent  j ,  respectively.   f  '  and  f „   are  the 

, .       .  .       .  - .  ,  - (m)  .    .     , 

corresponding  activity  coefficient  parameters,   f^   is  given  by 

equation  (7.27).   A  simplification  of  (7.29)  results  from  assuming 

that  C^ ^  is  a  mole  fraction  average  of  its  limiting  value  at  each  of 

the  pure  solvent  compositions,  i.e., 

V     (J) 

hi "  J  ^J^ll  ^^'^^^ 

Then,  equation  (7.29)  reduces  to 

M  (1  -  C  )^    M  X(l  -  C  )^ 

l,m    .^^     j    l.j    2(1  -  C^)    j^^  2(1  -  C.J 

-   2  -   2 

(pv,  )    M    (pv,  .) 

'   ,        P^  -  I  X.  iJr-  (7.31) 

The  above  equation  relates  the  Henry's  constants  deviation  to  the 
partial  molar  volumes  of  the  solute  and  isothermal  compressibilities 
in  the  mixture  and  the  pure  solvents.   These  quantities  are  more 
readily  obtainable  than  activity  coefficient  parameters;  thus 
equation  (7.31),  if  valid,  is  more  convenient  than  (7.29).   Two  other 
relations  for  the  Henry's  constant  deviation  are  of  interest  here.   The 
first,  used  by  Chueh  and  Prausnitz,    is 


M  E 

£nH^  ^  -  I   X^£nH,  ^  -  -  |:;:(3olvents)  (7.32) 

y 


'i.^-ivh.r-^^' 


14 
O'Connell   presents  a  relation  based  on  the  study  of  several 


experimental  systems 

+ 


M  X. 


£nH,    -  y  X.JLnH-  .  =  -r^ 


3x_ 


A23*2*3  (7.33) 


T.P.N^ 
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where  1050 |v„  ~  v  | 

^23  -   372 ^^-^"^^ 

T    V2V3 

$2  and  $_  are  component  volume  fractions  in  the  solute  free  solvent 
mixture.   An  interesting  and  perhaps  restrictive  feature  of  both 
(7.32)  and  (7.33)  is  that  no  solute  properties  are  involved   and  the  re- 
sulting Henry's  constant  deviation  is  independent  of  the  nature  of  the 
gas.   On  the  other  hand,  equations  (7.29)  and  (7.31)  do  account  for 
the  solute  interactions. 

The  Henry's  constant  deviations  determined  from  equations 
(7.29),  (7.31),  (7.32)  and  (7.33)  for  the  experimental  systems 
described  previously  are  shown  in  Table  7-7  and  compared  with  the 
experimental  results.   The  results  of  equations  (7.33),  (7.32)  and 
(7.31)  (with  experimental  DCF  integrals)  agree  with  each  other  to 
within  the  range  of  uncertainty  normally  associated  with  this  difference. 
They  are  also  in  reasonable  agreement  with  experiment  for  the  N^ 
systems  but  not  for  the  H„  system.   Equation  (7.31)  with  calculated 
DCF  integrals  also  gives  satisfactory  agreement  with  the  N„-Benzene- 
carbon  tetrachloride  system  but  the  trend  with  temperature  is  not 
entirely  consistent  with  experimental  values  in  the  N  -benzene-octane  . 
systems.   Henry's  constant  deviations  from  experimental  activity 
coefficient  parameters,  (7.29),  are  much  larger  than  the  experimental 
in  all  cases  but  one. 

On  the  basis  of  the  current  experimental  data,  expression 
(7.29)  must  be  considered  unsatisfactory  as  a  representation  of  the 
Henry's  constant  deviation.   Equation  (7.31)  may  be  as  accurate  or  more 
so  than  (7.33)  and  (7.32)   but  this  cannot  be  ascertained  on  the  basis 
of  the  present  limited  experimental  data. 
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7. A  Vapor-Liquid  Equilibrium  in  Mulricomponent  Systems 

The  test  of  most  practical  iir.portance  for  an  expression 
for  activity  coefficients  and  the  corresponding  reference  fugacities 
in  a  mixture  is  the  calculation  of  vapor-liquid  equilibrium.   The 
calculated  bubble  pressures  and  equilibrium  compositions  can  be 
compared  directly  with  experimental  values.   The  generalized 
correlations  for  f„  and  C^^,  which  provide  activity  coefficient 
parameters  and  the  Henry's  constant  deviation,  are  tested  here  vith 
the  use  of  the  one-fluid  model  for  the  solvent  mixtures.   Bubble 
pressure  calculations  were  made  for  two  experimental  systems,  H„-CH,- 
C^U^   at  116°K,  and  CH^-C^Hg-C^H^^  ^^  310.9°K.^°  The  representation 
of  vapor  and  liquid  phase  fugacities  followed  that  described  in 
Chapter  6.   The  excess  free  energies  of  the  solvent  mixtures 
CH^-C^H  and  ^2^R~^5^\2  ^^^^   determined  by  Chueh's  method.    However, 
this  quantity  could  be  obtained  from  an  independent  source,  if  possible 
or  desired.   Calculated  bubble  pressures  and  vapor  compositions  are 
shown  in  Table  7-8.   There  is  generally  good  agreement  between  the 
calculated  and  experimental  vapor  compositions,  despite  some  discrepancy 
in  the  calculated  bubble  pressure.   This  insensitivity  of  the  vapor 
mole  fractions  was  also  observed  in  the  calculations  for  binary 
systems. 

The  CH^-C^Hg-nC  H    (1-2-3)  system  is  discussed  first.   The 
solvent  volume  was  assumed  to  be  the  mole  fraction  average  of  the  pure 
component  volumes.   The  characteristic  parameters  for  the  pure 
components  each  required  less  than  1%   adjustment  to  match  the  experimental 
compressibility  with  that  predicted  by  (4.3).   (1  -  C   )  for  the  mixture 
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was  then  obtained  from  the  one-fluid  theory.   Next,  the  quantities 
C^  (v./v  )*   ,  j=2,3,  were  determined  from  the  experimental  v   .,  and 

equation  (A. 3).   These  were  plotted  on  the  graph  of  the  generalized 

o   *  *  62 
correlation  for  ^■lo^'^l^'^l^'      '  ^'1"^*^^°'^  (A. 6),  and  a  straight  line 

o   *  *  .62 

drawn  between  them.   The  mixture  quantity  C,  (v  /v^ )    was  determined 

im  m  i 

*   ,       o 

from  this  linear  relation  at  a  reduced  density  p  =  v  .  /v  .  .   C. ^ 

•^   m    mix  mix    11 

in  the  mixture  was  determined  similarly  from  the  generalized  correlation 

for  methane  in  solvents,  equation  (7.6).   The  derived  DCF  integrals 

C   ,  C,   and  C, ,  are  shown  in  Table  7-6  as  experimental  values, 
mm   Im      11 

Corresponding  values,  obtained  from  the  generalized  correlations  at 
the  reduced  density  of  the  solvent,  are  shown  as  calculated  values. 
The  mixed  solvent  Henry's  constants  and  related  f^   were  each 
determined  from  both  of  the  relations  (7.31)  and  (7.29).   The  results 
with  relation  (7.31)  are  the  first  set  in  Table  7-8;  those  of  (7.29) 
are  the  second  set.   Both  results  are  comparable,  and  the  K- factors 
agree  well  with  experiment,  although  the  calculated  pressures  are 
somewhat  lower  than  experiment.   Chueh's  method  yields  pressures  in 
closer  agreement  with  experiment.  The  third  set  of  calculations  for 
this  system  used  the  calculated  DCF  integrals  (Table  7-6)  in  relation 
(7.31)  to  obtain  the  Henry's  constants  and  activity  coefficient 
parameters.   The  resulting  bubble  pressures  and  component  K-f actors 
compare  very  poorly  with  experimental  values.  This  resulted  from  the 
large  negative  value  obtained  for  the  Henry's  constant  deviation. 

Three  solvent  mixture  compositions  were  studied  in  the  H.-CH,- 
C-H,,  (1-2-3)  system.   The  first  set  of  calculations  are  similar  to 
those  of  the  CH  -C  H  -C^-H   system.   The  only  differences  were  the  use 

H    J  O   J  iZ 


:6o 


of  critical  volumes  for  the  characteristic  parameters  of  the  pure 
solvents,  and  v^  .  .  being  taken  frop.  the  g'jneralizad  correlation, 
e.g.,  (4.6).   The  Henry's  constant  deviation  and  the  relation  f^ 
were  determined  with  experinental  C   and  C,   in  (7.31).   The  calculated 
pressures  differ  by  more  than  10%  from  the  experimental,  but  the  K- 
f actors  are  in  satisfactory  agreement.   The  second  set  of  calculations 
is  quite  different  from  the  first,  in  that  an  effective  f^'   is 
used  for  the  solvent  mixture.   This  is  prompted  by  the  fact  that  the 
experimental  C   data  for  H„  in  CH,  and  H  in  Ar  did  not  correlate 
well  with  those  for  other  H„-solvent  systems  in  equation  (7.7). 
Therefore,  an  effective  value  of  f„   was  determined  from  equation 
(7.5).   The  Henry's  constant  deviation  was  obtained  from  the  simplest 
of  the  solvent  excess  free  energy  expressions,  equation  (7.32). 
These  were  all  the  data  required,  no  other  binary  data  was  necessary. 
The  calculated  bubble  pressures  and  K-factors  show  remarkably  good 
agreement  with  experimental  results,  and  are  somewhat  better  than 
those  of  the  Chueh  method.   This  agreement  is,  indeed,  more  than  a 
fortunate  happenstance,  because  it  is  consistently  good  over  all 
three  solvent  compositions  considered. 

The  success  of  the  above  method  indicates  its  reliability  for  vapor 
liquid  equilibrium  calculations.   Its  usefulness  is  important  for 
physical  systems  in  which  the  detailed  data  required  by  Chueh 's  method 
are  not  available.   The  essential  difference  in  the  two  methods  is  the 
determination  of  the  activity  coefficient  parameter.   In  Chueh 's  method, 
it  is  calculated  from  sets  of  fitted  coefficients  for  each  binary 
pair,  whereas  in  the  present  method  it  is  obtained  from  generalized 
correlations. 
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The  method  can  be  extended  to  systeus  where  more  than  one 
component  is  treated  as  a  solute.   Smith   has  shovm  that  when  several 
solutes  are  present  at  low  concentrations  in  a  liquid  mixture,  the 
concentration  of  each  one  of  these  affects  only  its  oum  chemical 
potential  and  that  of  the  solvent  components.   In  simple  terms,  a 
solute  molecule  "sees"  only  solvent  molecules  around  it.   Then  the 
activity  coefficient  of  each  solute  can  be  represented  by  a  Margules 
(or  other)  type  expression,  involving  only  its  own  mole  fraction. 
Further,  its  Henry's  constant  in  the  liquid  mixture  will  be  independent 
of  the  other  dilute  solutes.   Bubble  pressure  calculations  can  be 
performed  with  the  solute  components  independent  of  each  other  in  the 
liquid  phase.   They  may  be  considered  interacting  in  the  vapor  phase, 
if  necessary. 


162 


REFERENCES  FOR  CHAPTER  7 


1.  M.  Orentlicher  and  J.  M.  Prausnitz,  Chem.  Eng.  Sci.  ,  _19_,  775 
(1964). 

2.  J.  F.  Connolly,  J.  Chem.  Phys. ,  36,  2897  (1962). 

3.  B.  H.  Sage  and  W.  N.  Lacey,  Thermodynamic  Properties  of  Lighter 
Hydrocarbons,  New  York,  API  (1950),  (1955). 

4.  J.  Shim  and  J.  P.  Kohn,  J.  Chem.  Engr.  Data,  7_,  3  (1962). 

5.  H.  Volk  and  G.  D.  Halsey,  J.  Chem.  Phys.,  32.  1132  (1960). 

6.  A.  Azarnoosh  and  J.  J.  McKetta,  J.  Chem.  Engr.  Data,  8_,  321 
(1963). 

7.  H.  H.  Reamer  and  B.  H.  Sage,  J.  Chem.  Engr.  Data,  4_,  98  (1959). 

8.  H.  J.  Aroyan  and  D.  L.  Katz,  I.E.C.,  43,  185  (1951). 

9.  A.  L.  Benham  and  D.  L.  Katz,  A.I.Ch.E.J.,  3_,    33  (1957). 

10.  J.  M.  Prausnitz  and  P.  L.  Chueh,  Computer  Calculations  of  High 
Pressure  Vapor  Liquid  Equilibria,  Prentice-Hall,  Englewood  Cliffs, 
N.J.  (1968). 

11.  H.  H.  Reamer  et  al. .  I.E.C.,  52,  534  (1950). 

12.  J.  F.  Connolly  and  G.  A.  Kandalic,  J.  Chem.  Engr.  Data,  7_,  137 
(1962). 

13.  J.  F.  Connolly  and  G.  A.  Kandalic,  CEP  Symposium  Series,  59 , 
#44  (1963). 

14.  J.  P.  O'Connell,  A.I.Ch.E.J.,  _17»  ^^3  (1971). 

15.  J.  E.  Smith,   Statistical  Mechanical  Theory  of  the  Henry's 
Constants,   Ph.D.  Dissertation,  University  of  Rochester  (1965). 


CHAPTER  8 
CONCLUDING  REMARKS 

The  isothermal  compressibility,  the  partial  molar  volume 
and  the  composition  derivative  of  the  chemical  potential  for  fluid 
mixtures  are  directly  related  through  a  rigorous  statistical 
mechanical  formalism  to  the  volume  integrals  of  molecular  distribution 
functions,  equations  (2.16)  through  (2.32).   Other  expressions  can 
be  obtained  for  the  internal  energy  and  heat  capacity. 

The  molecular  distribution  functions  and  their  volume  integrals 
are  complex  functions  of  the  macroscopic  state  of  the  system.   At 
present  lack  of  accurate  expressions  for  the  distribution  functions 
prohibit  their  use  in  calculating  thermodynamic  properties  directly 
from  molecular  parameters.   When  a  complete  theory  for  composition 
dependent  molecular  distribution  functions  does  become  available,  these 
expressions  should  be  more  generally  useful.   Functional  expansions  of 
the  expressions  for  thermodynamic  properties  are  made  about  limiting 
compositions  to  partially  circumvent  the  problem  of  describing  their 
possibly  complex  composition  dependence,  equations  (2.33)  -  (2.48). 
These  expressions  describe  the  limiting  slope  of  the  thermodynamic 
properties  quite  accurately  but  hold  over  restricted  composition 
ranges. 

Molecular  distribution  functions  of  different  orders  have  been 
related  through  density  and  composition  derivatives.   These  relations 
can  be  used  to  determine  the  higher  order  functions  from  the  density 
dependence  of  the  corresponding  lower  order  distribution  functions,  thus 
adding  to  the  present  understanding  of  molecular  interactions. 
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A  derived  first  order  perturbation  expression,  equation  (3.3), 
for  the  distribution  functions  uses  the  hard  sphere  fluid  as  a  base, 
with  an  attractive  perturbing  potential.   The  calculated  R.D.F.'s  for 
a  Lennard-Jones  6-12  fluid  in  the  liquid  region  agree  vith  nolecular 
dynamics  results.   The  density  dependence  of  the  reduced  DCF  integral, 
equation  (3.8),  is  similar  to  the  macroscopic  density  dependence  of 
the  isothermal  compressibility  of  pure  fluids.   This  is  the  basis  of 
a  qualitatively  accurate,  one-parameter  microscopic  theory  for  the 
reduced  DCF  integrals  for  pure  components  and  mixtures,  with  simple 
mixing  rules  for  mixture  parameters.   Thermodynamic  properties  at 
Infinite  dilution  in  subcritical  systems  as  calculated  from  this 
microscopic  theory  agree  only  qualitatively  vith  experimental  data, 
Table  3-2.   Quantitative  agreement  could  possibly  be  obtained  by 
Introducing  a  temperature  related  parameter  into  the  microscopic, 
and  per  force,  the  macroscopic  theories. 

The  one-parameter  macroscopic  generalized  correlation  for 
isothermal  compressibility  as  a  function  of  reduced  solvent  density, 
equation  (4.3),  is  a  good  representation  of  the  experimental  data. 
The  calculated  pressure  and  temperature  derivatives  of  the  compressi- 
bility, equations  (4.53)  and  (4.54),  also  agree  well  with  the  limited 
experimental  data  available.   In  the  absence  of  compressibility  or 
volumetric  data,  the  generalized  isotherm.al  equation  of  state, 
equation  (4.10),  is  useful  for  com.pression  calculations. 

The  compressibility  correlation  breaks  down  at  low  reduced 
densities  due  to  the  divergence  of  the  pure  liquid  isotherms  as  the 
critical  density  is  approached.   It  may  be  possible  to  correlate  the 
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divergent  behavior  in  this  region  by  introducing  a  temperature  (or 
energy)  related  parameter  into  the  macroscopic  nodel.   On  the  other 
hand,  the  present  correlation  is  quite  satisfactory  in  the  highly 
compressed  liquid  region,  and  might  be  used  to  interpret  the  molecular 
significance  of  the  coefficients  in  the  Tait  equation. 

Experimental  partial  molar  volumes  for  gases  in  liquid  solutions 
are  represented  by  equation  (4.6).  The  calculated  temperature  and 
pressure  derivatives  may  be  useful  for  preliminary  design  calculations, 
but  are  not  accurate  enough  for  an  understanding  of  the  molecular 
interactions  in  the  mixture  and  their  dependence  on  the  macroscopic 
state.   Possible  improvements  of  this  correlation  could  begin  with  the 
use  of  an  energy  parameter  for  the  solvent  obtained  from  data  at 
lower  reduced  densities. 

The  macroscopic  one-fluid  theory,  equation  (5.1),  with  a 
composition  dependent  characteristic  parameter,  equation  (5.2),  is 
found  satisfactory  to  represent  the  isothermal  compressibility  of 
binary  and  ternary  liquid  mixtures  by  equation  (4.3).   Improved 
agreement  is  obtained  by  adjusting  the  characteristic  parameters  of 
the  individual  components.   Little  error  is  introduced  by  assuming 
an  ideal  mixture  to  determine  its  volume.   In  binary  liquid  mixtures . 
the  assumption  that  allof^he  DCF  integrals  vary  linearly  between 
composition  extremes  yields  accurate  partial  molar  volumes  and  activity 
coefficients  at  intermediate  compositions,  though  inaccurate 
compressibilities.   However,  activity  coefficients  determined  from 
experimental  pure  component  compressibilities  and  partial  molar  volumes 
at  infinite  dilution  and  a  simple  mixing  rule  for  the  DCF  integrals  do 
not  represent  the  experimental  data  satisfactorily  (Table  5-6). 
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Improved  agreement  may  result  with  an  experimentally  deterin-.ned 
deviation  parameter  from  activity  coefficients  at  infinite  dilution. 
An  alternate  approach  is  the  development  of  correlations  for  partial 

molar  volumes  and  activity  coefficients  of  liquids  in  liquids,  i.e., 

_    o       o 

for  C^2  ^^^  *^i  1  '  ^^^^^  components  1  and  2  are  both  subcritical.   Such 

correlations  would  certainly  contribute  to  the  understanding  of 
distribution  functions  in  liquid  mixtures. 

A  detailed  analysis  of  the  H^-benzene  system  is  possible  with 
functional  expansions  about  the  pure  solvent  composition.   The  slopes 
of  the  compressibility  and  partial  molar  volume  at  zero  composition 
involve  both  pair  and  triplet  correlation  functions.   The  triplet 
correlation  functions  are  obtained  from  density  derivatives  of  the 
generalized  correlations  for  the  pair  functions.   The  calculated  zero 
slopes  do  not  agree  well  with  the  experimental  data,  the  largest 
source  of  error  being  the  density  derivative  of  C^_,  i.e.,  of  the 
partial  molar  volume.   The  Percus-Yevick  or  HNC  theories  do  not 
yield  triplet  correlation  functions  in  agreement  with  experiment. 

The  infinite  dilution  activity  coefficient  parameter  f„  and 
the  related  DCF  integral  C   ,  for  CH,  and  H  in  single  liquid  solvents 
are  represented  by  generalized  functions  of  the  reduced  solvent  density 
(7.4),  (7.5),  (7.6)  and  (7.7).   The  results  of  vapor-liquid  equilibrium 
calculations  using  equations  (7.4)  and  (7.5)  to  determine  activity 
coefficients  are  quite  encouraging.  Tables  7-2  and  7-3.   Calculated 
vapor  compositions  agree  well  with  experiment,  and  bubble  pressures 
generally  are  predicted  satisfactorily.   The  generalized  activity 
coefficient  parameter  offers  considerable  improvement  over  other 


167 


methods  wlierein  activity  coefficients  are  determined  from  coefficients 
characteristic  of  each  solute-solvent  pair. 

In  systems  of  gases  in  binary  solvents,  the  solvent  mixture  is 
represented  by  the  one-fluid  model.   The  volumetric  properties  of  the 
gas-solvent  mixture  are  described  well  by  the  functional  expansions 
with  the  experir.ental  DCF  integrals.   The  experimental  and  predicted 
DCF  integrals  become  more  negative  as  the  solvent  density  increases. 
The  integrals  of  the  solute-solute  interaction  are  always  less 
negative  than  those  of  the  solvent-solvent  interaction.  Table  7-6.   One 
concludes  that  the  size  and  energy  parameters  characterizing  the 
solute-solute  interaction  are  always  smaller  than  those  of  the 
solvent-solvent  interaction,  and  the  solute-solvent  parameters  lie  in 
between. 

Several  relations  between  the  Henry's  constants  of  a  gas  in  a 
mixed  solvent  and  in  the  corresponding  pure  solvents  are  analyzed 
with  experimental  data,  Table  7-7.   One  of  these  relations,  (7.31), 
Involving  the  compressibilities  of  the  pure  and  mixed  solvents,  and  the 
partial  molar  volumes  of  the  gases  in  these,  offers  advantages  in 
practical  applications,  but  its  validity  is  not  clearly  established  on 
the  basis  of  the  limited  experimental  data.   Further  investigation 
of  this  and  other  expressions  is  clearly  warranted  due  to  the 
increasing  use  of  mixed  solvents  in  industrial  applications. 

Activity  coefficients  of  CH,  and  H_  in  mixed  solvents  can  be 
determined  from  the  generalized  correlations,  equations  (7.6)  and  (7.7), 
when  the  mixed  solvent  is  represented  by  the  one-fluid  model.   The 
results  of  vapor-liquid  equilibrium  calculations  are  satisfactory. 
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Table  7-8.   Although  the  predicted  bubble  pressures  often  disagree  by 
20%  from  the  experimental  pressure,  the  calculated  K-factors  are  in 
close  agreement  with  the  experir.ental  data. 

This  5i.-;nsralized  method  for  the  activity  coefficient  of  the 
solute  can  be  extended  to  systems  of  two  or  more  solutes  present 
at  dilute  concentrations  since  the  activity  coefficient  of  each  solute 
is  independent  of  the  concentration  of  other  solutes. 
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APPENDIX  A 
MULTICOMPONENT  OR-N'STEIN-ZEMIKE  EQUATION 

For  an  m-component  system 


N.x 
where  n  (y)  =  number  density  =  — — 


(A..1) 
(A. 2) 


Multiplying  through  by  d_r'  and  integrating 


Ihr^   ^ij  ""   \     ^  jCSi/^'Z)  -  IJdX  j  c   (y,r')dr'     (A.3) 
•^       ■*    £=l,m   J    -J  J    J 


where 


G,.  -p 


CgijCr.r')  -  l]dr' 


(A. 4) 


and  n  is  taken  independent  of  position  in  a  homogeneous  system.   In 
the  integral  of  equation  (A.3),  2.   ^^d  _r'  are  dummy  indices 


1g 


_  V  r.         V       V 

"  N  ^ii  "^  „  ^    ''Z  N  ^i£  N  ^ 


£=l,m 


^j 


G  =  C  +  G  X  C 


(A.5) 
(A.  6) 


J.  L.  Lebowitz,  Phys.  Review,  133,  A895  (1964) 
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APPENDIX  B 
RECURRENCE  RELATIONS  FOR  DCF  INTEGRALS 


For  the  change  of  variables  from  constant  pressure  to  constant 


volume 
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We  also  have 


3c 


(2) 


3N, 


(B.A) 


T.V.N 


m=^k 


Equations  (B.3)  and  (B.4)  are  used  to  derive  the  recurrence  relation 

between  C.  ,  and  C.  .,  . 
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and  using  Liebnitz's   rule  for  the  partial  derivative  of  the  integral 
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dC 


(2) 
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But  c. .  '(V)  is  the  microscopic  d.c.f.  when  the  separation 
between  molecules  is  of  the  order  of  the  macroscopic  volume  of  the 
systerj;  it  is  zero. 
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where  equation  (B.3)  has  been  substituted  into  (B.6).   Then 
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The  quantities  C.  and  C.,.  can  also  be  related  through  density 
ij      ijk  *        ^ 

derivatives.   We  have,  from  their  definitions. 
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Using  equation    (B.A)    above 
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APPENDIX  C 

FIRST-uyjJHR  PERTURBATION  THEORY  FOR  THE  DISTRIBUTION 
FUNCTIONS  OF  A  DENSE  FLUID 


Percus  shoved  that  molecular  distribution  functions  are 

functional  derivatives  of  the  molecular  density  at  any  point  in  the 

fluid  with  respect  to  a  change  in  potential  on  the  fluid.   This 

provides  a  precise  relation  between  molecular  distribution  functions 

of  different  orders.   The  change  in  potential  may  be  due  to  the 

introduction  of  an  external  molecule  either  the  same  as,  or  different 

from,  the  molecules  of  the  fluid  under  consideration.   The  functional 

derivatives  are 

6n  (r') 
g[_g^(Q)j  n^(r')c(0  -  r')  +  n^Cr' .0)  -  n^(r')n^(0)  (C.l) 

^^^/"^T)^    [n^(r')]~^6(r  -  r')  -  C2(r,r')  (C.2) 

Equation  (C.l)  relates  the  change  in  n^(r'),  the  single-particle 
density  at  r',  due  to  a  change  in  potential  at  0,  U(0),  in  terms  of 
one-  and  two-particle  distribution  functions  n^(0),  n^(r'),  and 
n-CO.r').   The  functional  inverse  shoiim  in  equation  (C.2)  defines 
c_(r,r'),  the  second-order  direct  correlation  function.   Multiplying 
equations  (C.l)  and  (C.2)  and  integrating  over  r'  yields  the  Ornstein- 
Zernike  relation 


[n2(0,r)/n^(r)]  -  n^(0)  =  n^(0)c2(0,r)  + 


C2(r,r')[n2(0,r') 


-  n^(r')n^(0)]dr'  (C.3) 

In  addition,  using  the  functional  relations  (C.l)  and  (C.2),  the 
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quantity  n  (r)  exp[SU(r)],  which  is  used  in  the  Percus-Yevick  theory, 
may  be  expanded  in  terms  of  a  perturbation  in  the  density  to  give 

n*(r)exp[BU*(r)]  =  n°(r)exp[3U°(r) ] [1  +  j  c°(r ,r ' ){n*(r' )  -  n°(r')}dr' 


-i 


(c^Cr.r'.r")  +  c^Cr .r' )c^(r,r") }{n^(r ' )  -  n°(r')} 
{n^(r")  -  n°(r")}dr'dr"  +  ...]  (C.4) 


* 


o 


In  equation  (C.A),  the  superscript   refers  to  the  unperturbed  fluid 

* 
and  the  superscript  asterisk  (  )  to  the  real  fluid.   If  the  change  in 

external  potential  causing  the  density  fluctuation  at  r'  is  due  to 

placing  a  particle  or  changing  the  kind  of  particle  at  the  origin  of 

coordinates,  then  the  first-order  distribution  functions  n  (r') 

and  n  (r')  are  second-order  distribution  functions,  since  they 

represent  the  probability  of  a  particle  being  at  r'  when  it  is  known 

that  there  is  a  particle  at  the  origin,  i.e., 

nj_(r'/0)  =  n2(0,r')/n^(0)  =  n^(r '  )g2(0,r '  )  (C.5) 

where  n^^Cr'/O)  is  the  one-particle  distribution  function  at  r'  in 
the  presence  of  a  particle  at  0.   Then,  to  the  first  order  in  density, 
equation  (C.5)  becomes 
* 


g2(0,r)exp[eu  (r) ]  =  g°(0,r)exp[6U°(r) ] [1  + 


c°(r,r')n^(r') 


{g2(0.r')  -  g°(0,r')}dr']  (C.6) 


We  shall  now  explicitly  use  equation  (C.6)  in  constructing  a 
real  molecular  interaction  from  a  hard-sphere  interaction.  Consider 
the  two  states  shown  in  Fig.  C-1.   In  state  1,  the  fluid  is  composed 
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State  1 


State  2 


Q   Molecule  of  real  fluid 

O   Molecule  of  hard  sphere  fluid 


Fig.  C-1. 
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of  molecules  of  the  "real"  fluid  and  a  molecule  whose  interaction 
with  the  others  can  be  represented  by  a  hard-sphere  potential  is 
situated  at  the  origin.   [A  molecule  at  r  and  one  at  r'  are  shown; 
the  total  potential  at  r  due  to  all  the  other  molecules  in  the  system 
is  V(r).   The  interactions  between  the  molecules,  e.g.,  c^Cr.r'), 
are  those  typical  of  the  "real"  fluid.]   In  state  2,  a  molecule  of 
the  "real"  fluid  replaces  the  "hard-sphere"  molecule  at  the  origin. 
Equation  (C.6)  may  be  applied  to  the  change  in  state  from  1   2: 

g*(0,r)exp[6[V(r)  +  4)*(0,r)]]  =  g°(0,r)exp[6[V(r)  +  <t)°(0,r)]] 

•  {1  +  [  c*(r,r')[n*(r')g*(0,r')  -  n°(r')g°(0.r' ) ]dr' }      (C.7) 

Writing  the  Omstein-Zernike  equation  for  the  states  1  and  2 

c*(r,r')n°(r')[g°(0,r')  -  l]dr'        (C.8) 
c*(r,r')n*(r')[g*(0,r')  -  l]dr'        (C.9) 


g°(0,r)  -  1  =  c°(0,r)  + 


g2(0,r)  -  1  =  C2(0,r)  + 

We  shall  now  assume  that  g°  =  g^^,  the  true  hard-sphere  distribution 
function;  i.e.,  it  is  the  solution  of  equation  (C.8)  when  C2(r ,r ' )n^(r' ) 
is  replaced  by  pc^(r,r')  within  the  integral  of  equation  (C.8), 
and  p  is  the  real  fluid  density.   In  equation  (C.9),  all  quantities 
refer  to  the  real  fluid  (rf)  and  nj^(r')  =  p.   Subtracting  (C.8)  from 
_(C.9)  and  substituting  in  (C.7),  we  get 

g''^r)exp[B/^(r)]  =  g^^r)exp[6^^r)  ]{1  +  g'^^r)  -  g^^r) 

-  c^'^r)  +  c^^r)}  (CIO) 

Equation  (C.IO)  is  valid  at  all  values  of  r,  so  the  intermolecular 
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potential  and  distribution  functions  of  the  real  fluid  are  different 

from  those  of  the  hard-sphere  fluid  over  the  entire  range.  In  the 

range  r  <  a,  v/here  a  corresponds  to  the  hard-core  distance  of  the 
hard-sphere  fluid,  we  have 

g^'^Cr)  =  g^^(r)exp[B:f^^(r)][l  -  c^^(r)  +  c^^r)]/{exp[B<J,''^(r)] 

-  g^^(r)exp[B(^^^(r)]}  (C.ll) 

The  quantities  c   (r)  and  c   (r)  and  the  g   (r)exp{B(j!   (r)}  are 

rf 
all  finite.   Thus,  g   (r)  varies  inversely  as  the  exponential  of  the 

real  fluid  potential  and  quickly  becomes  very  small  as  r  decreases. 

The  distribution  functions  in  the  real  fluid  are  now  completely  defined 

in  terms  of  the  hard-sphere  distribution  functions  and  the  differences 

in  the  intermolecular  potentials. 

The  reference  hard-sphere  fluid  may  be  chosen  to  be  an  ideal 

gas ,  where  - 

g^^(r)  =  1,  (C.12) 

c^^(r)  =  <{,^^(r)  =  0       for  all  r 

Then,    _g^^(r)exp[B<^''i(r)]  =  g^^(r)  -  c^'^lr} (C.13) 

which  is  the  relation  of  the  usual  first-order  Percus-Yevick  theory. 
Thus,  equation  (CIO)  is  a  general  form  of  the  Percus-Yevick  equation 
with  a  hard-sphere  reference. 

:  -    Equation  (CIO)  can  be  used  for  any  intermolecular  potential 
and  reference  distribution  function  with  the  procedure  described 

above.   When  it  is  substituted  into  the  Ornstein-Zernike  equation,  an 

rf        rf 
integral  equation  for  either  g   (r)  or  c   (r)  results. 


APPENDIX  D 
TEMPKRATUPviL  DERIVATIVE  OF  LIQUID  DENSITY  AT  CONSTANT  PRESSURE 

FroiT:  aquation  (4.8) 


P  -  P^ 
KTp^ 


[FCp^)  -  l]dD2 


(D.l) 


Taking  a  partial  derivative  w.r.t.  T  at  constant  P,  and  using  Liebnitz's 
rule 


d_ 
dx 


B(>:) 


B 


F(x^t)dc 


[  i£ff^  .t  .  r(x.BU))  g 


A(x) 


-  F[x,A(x)]  1^   (D.2) 
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0  - 

r  sp^  ^ 

P  _  pS    (Fl(p)  • 

-  1) 

aP2] 

3T 

RTP2 

1  3T   ^ 

P 

2  *         * 
RT  P2        P2 

(f^Cp^)  -  1) 

9T 

* 
P2 

P 

(D.3) 


Now 


[  aP2^ 

3T 

P 

3T  ■ 
3P 

3P_ 
3T 


=  -1 


(D.-4) 


3p2 
3T 

P 

r  SP  ^ 

P2 

3T 

P 

3P  ' 

3T  _ 

P2  RT[F^(p2)  -  1] 


(D.5) 


From  (D.3) 


3p, 


[F^(P2) 


1]  -  [F^(P2)  -  1^  dT"  RT  dT 


p  -  P' 

RT*^ 


(D.6) 


which  gives  the  temperature  dependence  of  the  density  a  compressed  liquid, 


179 


180 


The  pressure  variation  with  temperature  at  constant  density  can  be 
determined  from  (D.5)  snd  (D,6). 


APPENDIX  E 

EXTRACTION  OF  HENRY'S  CONSTANTS  AND  ACTIVITY  COEFFICIENTS 
FROM  VAPOR-LIQUID  EQUILIBRIUM  DATA 


For  binary  vapor-liquid  equilibrium 


^1*1^  =  Xj^Yj^CT,P2,x^)H^  ^iT,?^2,x)ex'p 


f^  V  (T.x) 

R^   dP  (E.l) 

s 


^2 
s 
where  P_,  the  saturation  pressure  of  the  solvent,  is  the  reference 

pressure  for  the  fugacity  of  componentl;  v^(T,x,P)  is  the  partial 

molar  volume  of  the  solute  at  composition  x.   At  low  concentration  x^ , 

V  (T,x,P)  can  be  approximated  by  the  partial  molar  volume  at  infinite 

dilution,  i.e., 

v°(T,P,x)  =.  v°(T,P)  (E.2) 

Then    ^       ^  f^   v°(T,P) 


y^^^?   -  x^Y^(T,x^)H^  2^'^^^^P 
The  solvent  fugacity  is 


RT 
s 


■'•      dP  (E.3) 


v  ?^ 


P  o 

RT 


^  dP  (E.4) 


-  s         -     ^' 

and  f?  =  P^*^  (E.5) 

^     2  2 

where  P„ ,  the  solvent  saturation  pressure,  and  $_,  its  fugacity 
coefficient  in  the  pure  liquid  state  are  calculated  from  generalized 
xForrelations.    The  two  (or  three)  term  Margules  equation  is  chosen  to 


J.  K.  Prausnitz  et  al. ,   Computer  Calculations  for  Multlcomponent 
Vapor  Liquid  Equilibrium,  Prentice-Hall,  Englewood  Cliffs,  N.J.  (1968). 
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describe  the  composition  dependence  of  the  liquid  phase  activity 
coefficients 

£nY*(x^,T)  =  A(x^  -  2x^)  +  B(x^  -  |  x^)  (E.6) 

2     3 
iny^  =  Ax^  +   Bx^  (E.7) 

Vapor  phase  properties  are  represented  by  the  Redlich-Kwong 

2 
equation,  as  modified  by  Chueh  and  Prausnitz.    They  use  the  general 

Redlich-Kwong  equation 

V  -  b   t""'  v(v  +  b) 

but  determine  the  parameters  empirically.   Mixing  rules  are 
empirically  derived  for  the  critical  properties  of  mixtures.   The 
relation  for  the  fugacity  coefficient  is 

jin$   =  An— ^  +  ^-  ^n^^^^  .  2  !  y.a.,/RT^/2^ 
k,      V  -  b   v-b        V       .  V,   1  ik 

1=1 

ab 

k     oV  +  b     b       oPv  /-n  n\ 

The  fugacities  of  the  two  components  in  both  phases  are  now 
completely  defined  in  terms  of  experimental  quantities  and  the 
parameters  to  be  evaluated. 

At  each  state  point,  the  temperature,  bubble  pressure,  liquid 
compositions,  and  partial  molar  volumes  for  each  component  are  supplied 


2 
P.  L.  Chueh  and  J.  M.  Prausnitz,  Computer  Calculations  for 

High  Pressure  Vapor-Liquid  Equilibria,  Prentice-Hall,  Englewood 

Cliffs,  N.J.  (1970). 
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to  the  computer  program.   The  calculational  algorithm  is  schematically 
presented  in  Fig.  E-1.   At  the  given  temperature  and  pressure,  and 
guessed  values  of  the  parameters,  the  liquid  composition  is  adjusted 
until  the  vapor  mole  fraction  sum  to  unity.   The  sum  square  error 
in  liquid  composition  is  calculated  at  this  point.   The  program 
searches  over  different  parameter  sets  until  a  minimum  in  the  sum 
square  error  of  the  liquid  composition  is  obtained.   (The  searching 
algorithm  (RMINSQ)  is  based  on  the  Fletcher  Powell  method^ 

For  a  binary  solvent  system,  the  two  suffix  Margules  equation 
represents  the  excess  free  energy  of  the  solvent  mixture,  e.g., 

RTJlny^  =  (2D  -  C)x^  +  2(C  -  D)xi  (E.IO) 

The  excess  free  energy  in  the  1-2-3  mixture  relative  to  the  2-3 
mixture  is  again  represented  by  a  Margules  equation  which  gives 

RT£nY*jj^  =  A^23<4  '  ^""l^  ^^'-^^^ 

RTinYj^^  =  A^23^i  +  ^^D  -  OX^  +  2  (C  -  D)X^  (E.12) 

RTJlnY2j^  "  '^123^1  "^  ''^^  "  ^^'^l   "*"  ^  (D  "  C)X^  (E.13) 

X^  and  X_  are  mole  fractions  in  the  solute-free  mixture. 

The  fugacity  reference  states  for  the  two  solvent  components 
are  their  respective  saturation  pressures.   For  the  solute  component, 
it  is  the  infinitely  dilute  solution  at  its  equilibrium  pressure, 
i.e.,  the  saturation  pressure  of  the  liquid  mixture.   The  Henry's 
constant  and  solute  activity  coefficients  are  referred  to  this 
pressure  and  not  to  zero  pressure  as  in  Prausnitz  and  Chueh. 
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APPENDIX  F 

PERCuS-YEVICK  AND  HYPERNZTTED  CHAIN  APPROXIMATIONS 
FOR  TRIPLET  CORRELATION  FJNCTIONS 


The  first-order  PY  and  KNC  expression  for  the  pair  distribution 
functions  are  obtained  from  general  density  expansion  by  setting 
terms  of  second  and  higher  order  to  zero."^'   These  conditions,  in 
turn,  provide  relations  for  the  triplet  and  higher  distribution 
functions.   For  an  M-component  mixture,  under  the  Percus-Yevick 
approximation,  we  have 


]=1  L=l  -*    •  -     -J      J 


I^)-c{2)(^.rj)e{2)(^.^)] 


•  UijC^i.lj)  -  lHg.j^(r^.r^)  -  lldr^dr^  =  0  (F.l) 

for  all  i 

Then,  setting  the  integrands  to  zero,  and  assuming  that  [g .  .  (r..  ,£.)  -  1] 
and  ig. .(£.,£,)  -  1]  are  not  zero,  there  results 

^^I^Sll  -"  4l^  ^  2x^x^(C^^.  +  C^^C^^)  +  x2(C^22  ^  ^L^  =  ° 

(F.2) 
and 

4^-112  ■"  ^12>  ^  2x^X2(C^22  *  ^12^2^  ^  ^2^^222  ^  ^22^      ^^'^^ 


J.  K.  Percus  in  Classical  Fluids,  edited  by  H.  L.  Frisch 
and  J.  L.  Lebovitz,  Eenjaaiin,  New  York,  1964. 

2 
J.  L.  Lebowitz,  Physical  Review,  133,  A895  (196A). 
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in  a  binary.   For  infinite  dilution  of  component  1  in  2 


and 


C°   =  -  C°^ 
122     ^12 


'222     ^22 


(F.4) 


(F.5) 


Similarly,  for  the  HInC  approximation 


M   M 


i=l  k=l   -•    J  J    -1    ■'-      2      ^-        IJ   1   J 


l][gi^(r,.r,^) 
for  all  i 


leading  to 


and 


4Sll  +  2^l'^2^112  +  ^^122  =  ° 


^1^112  ^  ^^1^2^22  +  ^2^22  =  ° 


IJdr.dr,  =  0 
(F.6) 


(F.7) 
(F.8) 


Again,  for  infinitely  dilute  1  in  2 


and 


^22  =  ° 


S22  =  ° 


(F.9) 


(F.IO) 
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